Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,364,443 T→C 12.1% L594L (CTT→CTC yhcD → predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,364,4430TC12.1% 221.6 / 19.7 84L594L (CTT→CTCyhcDpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base C (4/6);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACAT  >  W3110S.gb/3364375‑3364512
                                                                    |                                                                     
cACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGc                                                                     >  1:3279006/1‑71 (MQ=255)
cACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGc                                                                     >  1:2938234/1‑71 (MQ=255)
cACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGc                                                                     >  1:2397828/1‑71 (MQ=255)
cACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGc                                                                     >  1:2241391/1‑71 (MQ=255)
   gacgaACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGc                                                                  <  1:1661724/71‑1 (MQ=255)
          gTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTCGCCGCAGGAGCa                                                           <  1:1130731/71‑1 (MQ=255)
          gTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTCGCCGCAGGAGCa                                                           <  1:3298018/71‑1 (MQ=255)
          gTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTCGCCGCAGGAGCa                                                           <  1:2468175/71‑1 (MQ=255)
            gAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAg                                                         >  1:3463460/1‑71 (MQ=255)
               aCGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCgg                                                      <  1:1251776/71‑1 (MQ=255)
                    cTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:2530977/71‑1 (MQ=255)
                    cTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:2487030/71‑1 (MQ=255)
                    cTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:2296336/71‑1 (MQ=255)
                     ttGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:995446/70‑1 (MQ=255)
                       gAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGa                                                <  1:306342/69‑1 (MQ=255)
                         aCGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGc                                            <  1:3042098/71‑1 (MQ=255)
                             cAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTa                                        <  1:2608571/71‑1 (MQ=255)
                             cAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTa                                        <  1:1144436/71‑1 (MQ=255)
                               gTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:743628/60‑1 (MQ=255)
                               gTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:3163410/60‑1 (MQ=255)
                               gTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:1295587/60‑1 (MQ=255)
                               gTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:1099487/60‑1 (MQ=255)
                               gTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:2230309/60‑1 (MQ=255)
                               gTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTg                                                 <  1:36234/60‑1 (MQ=255)
                                 aTCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGAt                                               >  1:2805263/1‑60 (MQ=255)
                                 aTCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACtt                                     >  1:1624916/1‑70 (MQ=255)
                                 aTCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTa                                    >  1:2405692/1‑71 (MQ=255)
                                  tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAg                                                   >  1:1709967/1‑55 (MQ=255)
                                  tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACtt                                     >  1:201436/1‑69 (MQ=255)
                                  tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTa                                    >  1:283778/1‑70 (MQ=255)
                                  tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTa                                    >  1:3154207/1‑70 (MQ=255)
                                  tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTa                                    >  1:3568751/1‑70 (MQ=255)
                                  tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTa                                    >  1:471021/1‑70 (MQ=255)
                                  tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTCGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTa                                    >  1:1299069/1‑70 (MQ=255)
                                  tCAGTTCCAATCTCCATTACACCCAAGTCAGCCTCGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTa                                    >  1:1275670/1‑70 (MQ=255)
                                    aGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACaa                                 >  1:165584/1‑71 (MQ=255)
                                    aGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACaa                                 >  1:3201648/1‑71 (MQ=255)
                                    aGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACaa                                 >  1:1690668/1‑71 (MQ=255)
                                    aGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACaa                                 >  1:260522/1‑71 (MQ=255)
                                    aGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACaa                                 >  1:258003/1‑71 (MQ=255)
                                      ttCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGt                                                  <  1:2955401/52‑1 (MQ=255)
                                       tCCAATCTCCATTACACCCAAGTTAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACgg                              >  1:1663608/1‑71 (MQ=255)
                                       tCCAATATCCATTACACCCAAGTTAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACgg                              >  1:3509/1‑71 (MQ=255)
                                        ccAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGa                                                >  1:2505650/1‑52 (MQ=255)
                                        ccAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACgg                              >  1:3294005/1‑70 (MQ=255)
                                        ccAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACgg                              >  1:2137276/1‑70 (MQ=255)
                                        ccAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACgg                              >  1:3325128/1‑70 (MQ=255)
                                        ccAATCTCCATTACACACAAGTCAGCCTAGCCCCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACg                               >  1:3631840/1‑69 (MQ=255)
                                          aaTCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTAc                           <  1:399987/71‑1 (MQ=255)
                                            tctcCATTACACCCAAGTCAGCCTCGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTAc                                   <  1:2708497/61‑1 (MQ=255)
                                              tcCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATg                       >  1:280824/1‑71 (MQ=255)
                                              tcCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATg                       >  1:2037093/1‑71 (MQ=255)
                                              tcCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATg                       >  1:183796/1‑71 (MQ=255)
                                                cATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGa                         >  1:1282860/1‑67 (MQ=255)
                                                cATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGa                         >  1:2190656/1‑67 (MQ=255)
                                                cATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGa                         >  1:2865908/1‑67 (MQ=255)
                                                cATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGa                         >  1:2138316/1‑67 (MQ=255)
                                                  ttACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAg                   >  1:2838209/1‑71 (MQ=255)
                                                    acacCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGa                         >  1:1054059/1‑63 (MQ=255)
                                                    acacCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGa                         >  1:2768343/1‑63 (MQ=255)
                                                    acacCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGa                         >  1:265686/1‑63 (MQ=255)
                                                    acacCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGa                         >  1:2438764/1‑63 (MQ=255)
                                                    acacCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAgg                  <  1:709279/70‑1 (MQ=255)
                                                    acacCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAgg                  <  1:1496496/70‑1 (MQ=255)
                                                    acacCCAAGTCAGCCTCGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggt                 >  1:2630721/1‑71 (MQ=255)
                                                    acacCCAAGTCAGCCTCGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggt                 >  1:1142396/1‑71 (MQ=255)
                                                    acacCCAACTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACg                          >  1:1559791/1‑62 (MQ=255)
                                                     cacCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggt                 <  1:1320252/70‑1 (MQ=255)
                                                      acCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggtgg               <  1:3561038/71‑1 (MQ=255)
                                                      acCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggtgg               <  1:764899/71‑1 (MQ=255)
                                                      acCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggtgg               <  1:1833594/71‑1 (MQ=255)
                                                      acCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggtgg               <  1:2988728/71‑1 (MQ=255)
                                                      acCCAAGTCAGCCTTCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggtggt              <  1:3238981/71‑1 (MQ=255)
                                                      acCCAAGTCAGCCTCGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggtgg               <  1:2238896/71‑1 (MQ=255)
                                                      acCCAAGTCAGCCTCGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggtgg               <  1:47157/71‑1 (MQ=255)
                                                       cccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGt                      <  1:20724/63‑1 (MQ=255)
                                                       cccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGt                      <  1:3497093/63‑1 (MQ=255)
                                                       cccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGt                      <  1:2498727/63‑1 (MQ=255)
                                                       cccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGt                      <  1:2598216/63‑1 (MQ=255)
                                                       cccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGt                      <  1:1585116/63‑1 (MQ=255)
                                                       cccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGt                      <  1:120987/63‑1 (MQ=255)
                                                       cccAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAggtggt              <  1:970929/71‑1 (MQ=255)
                                                               aGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCg       <  1:1446097/70‑1 (MQ=255)
                                                               aGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCg       <  1:2095659/70‑1 (MQ=255)
                                                                   ttGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACAt  >  1:663806/1‑71 (MQ=255)
                                                                    |                                                                     
CACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTCAGGTGGTATTGCCGTACAT  >  W3110S.gb/3364375‑3364512

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: