Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,967,225 C→G 8.1% intergenic (+52/+67) PP_1762 → / ← yciK hypothetical protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,967,2250CG 8.1% 275.7 / 12.5 86intergenic (+52/+67)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (36/43);  new base G (4/3);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGG  >  NC_002947/1967077‑1967374
                                                                                                                                                    |                                                                                                                                                     
cacgagTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCAcc                                                                                                                                                      <  1:824753/147‑1 (MQ=255)
           acccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                           >  1:1821810‑M1/101‑150 (MQ=255)
                aTGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGAtt                                                                                                                                      <  2:1735775/150‑1 (MQ=255)
                     cctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGcc                                                                                                                                 <  2:2601677/150‑1 (MQ=255)
                      ctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGccc                                                                                                                                <  1:4134521/150‑1 (MQ=255)
                       tgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                                                                                               <  2:1732110/150‑1 (MQ=255)
                       aagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                                                                                               >  1:4367404‑M1/89‑150 (MQ=255)
                        agagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctg                                                                                                                              <  1:3764886‑M1/63‑1 (MQ=255)
                         ttCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGc                                                                                                                             <  1:1834668‑M2/150‑2 (MQ=255)
                           gatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgcct                                                                                                                           >  2:822444‑M1/85‑150 (MQ=255)
                           ccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcct                                                                                                                           >  2:2370566‑M2/1‑147 (MQ=255)
                             tcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                                         >  2:617367‑M1/83‑150 (MQ=255)
                              ggAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtg                                                                                                                        >  1:1623930‑M2/1‑144 (MQ=255)
                              ggAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtg                                                                                                                        >  1:2785016‑M2/1‑144 (MQ=255)
                               gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgg                                                                                                                       >  2:1834668‑M2/1‑143 (MQ=255)
                               gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgg                                                                                                                       >  2:1281101‑M2/1‑143 (MQ=255)
                               atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                       >  2:709298‑M1/81‑150 (MQ=255)
                               atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCCCAGGCAGGCCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                       >  1:5071304‑M1/81‑150 (MQ=255)
                               atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                       >  1:2328536‑M1/81‑150 (MQ=255)
                               atgccggtgtacctgcacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCCAGGATTGCGCCCTGCCTGTgg                                                                                                                       >  2:1637006‑M1/81‑150 (MQ=255)
                                tgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                             >  1:2447975‑M1/80‑143 (MQ=255)
                                tgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                             <  2:2447975‑M1/64‑1 (MQ=255)
                                 gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGAGTGCGCCCTGcctgtggga                                                                                                                     >  1:4169069‑M2/1‑141 (MQ=255)
                                 gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgtccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                     >  2:4021396‑M1/79‑150 (MQ=255)
                                 gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTTGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                     >  1:816259‑M1/79‑150 (MQ=255)
                                 gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                     >  2:4413078‑M1/79‑150 (MQ=255)
                                 gccggtggacctggacctgctggggccggacagcaccggggtgaacgggcaggcgttcaacgcccaagccacctgcagCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCACCCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                     >  2:4070167‑M1/79‑150 (MQ=255)
                                  ccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAg                                                                                                                    >  2:3764886‑M1/78‑150 (MQ=255)
                                   gcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagc                                                                                                                   >  1:1722670‑M2/1‑139 (MQ=255)
                                    ggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCg                                                                                                                  >  2:4863141‑M1/76‑150 (MQ=255)
                                    ggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCg                                                                                                                  >  1:2324250‑M1/76‑150 (MQ=255)
                                    cTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagtg                                                                                                                  >  1:4310178‑M2/1‑138 (MQ=255)
                                    cTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcg                                                                                                                  >  2:2342340‑M2/1‑138 (MQ=255)
                                     tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgg                                                                                                                 >  2:3095719‑M2/1‑137 (MQ=255)
                                     tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgg                                                                                                                 >  2:494488‑M2/1‑137 (MQ=255)
                                     tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                                 >  2:1199323‑M2/1‑137 (MQ=255)
                                     tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                                 >  1:612301‑M2/1‑137 (MQ=255)
                                     tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                                 >  2:2124840‑M2/1‑137 (MQ=255)
                                      gcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcggg                                                                                                                >  2:4145372‑M2/1‑136 (MQ=255)
                                         tgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgggcaa                                                                                                             >  2:757647‑M2/1‑133 (MQ=255)
                                         tgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgggcaa                                                                                                             >  1:175508‑M2/1‑133 (MQ=255)
                                         tgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaa                                                                                                             >  2:4124384‑M2/1‑133 (MQ=255)
                                                    gggCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGGTTGCCCCCTCCCACAGGGTCCCCGCAACCCCCCAGGAGTGCGCCCGGccggtgggagcgggccagcccccccaga                                                                                                  <  2:3936963‑M2/150‑29 (MQ=255)
                                                              aaGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTTCCCGCTCCCACAGGGTCACCGCAACCCTCAAGGATTTCGCCCTGCCAGTGGGAGCGGGCAAGCCCCCGAAGAGGCCGATGCa                                                                                        >  2:3687647/1‑150 (MQ=38)
                                                               aGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGTTCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTAAAGGATTGCGCCCTGCCTGTGGGAGCGGGAAAGCCCGCGAAGAGGCCGATGCAg                                                                                       <  2:360771/150‑1 (MQ=255)
                                                                                                                 GCCTCTGCGCGGGCGTGCCCGCTCCCACCGGATCCCCGCAACCCTCAAGGTTTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGccc                                     >  2:2856058‑M1/1‑150 (MQ=255)
                                                                                                                        CGCGGGCTGGCCCGCTCCCCCCGGATCCCCGCAACCCCCAAGGAGTGCGCCCGGCCTGTGGGCGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcagg                              >  2:4783786‑M1/1‑150 (MQ=255)
                                                                                                                         GCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggt                             <  2:3760621‑M1/150‑1 (MQ=255)
                                                                                                                               TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                       <  2:502165‑M1/150‑1 (MQ=255)
                                                                                                                               TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTGTATTGGGCGTTGAACGCCTGCCCGTTCACCCCGGTGCTGTCCGGCCCCTTcaggtacaggt                       <  2:3196967‑M1/150‑1 (MQ=255)
                                                                                                                                 GCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCaggtacaggtac                     <  2:4266871‑M1/150‑1 (MQ=255)
                                                                                                                                  CCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAggtacaggtaca                    <  1:2288707‑M1/150‑1 (MQ=255)
                                                                                                                                   ccGCTCCCACAGGATCACCGCAACCATCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgca                   <  1:463148‑M2/150‑108 (MQ=255)
                                                                                                                                   CCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                   <  2:4710050‑M1/150‑1 (MQ=255)
                                                                                                                                   CCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                   <  1:3081414‑M1/150‑1 (MQ=255)
                                                                                                                                   CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATc                                 >  1:720505‑M1/1‑136 (MQ=255)
                                                                                                                                   CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATc                                 <  2:720505‑M1/136‑1 (MQ=255)
                                                                                                                                   CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                   <  1:2665218‑M1/150‑1 (MQ=255)
                                                                                                                                    CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc                  <  2:112904‑M1/150‑1 (MQ=255)
                                                                                                                                    CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc                  <  1:4979621‑M1/150‑1 (MQ=255)
                                                                                                                                    CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc                  <  1:4049659‑M1/150‑1 (MQ=255)
                                                                                                                                    CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc                  <  1:3153429‑M1/150‑1 (MQ=255)
                                                                                                                                    CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc                  <  1:217940‑M1/150‑1 (MQ=255)
                                                                                                                                     gCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagc                 <  2:1567509‑M2/150‑110 (MQ=255)
                                                                                                                                     gCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagc                 <  1:4574084‑M2/150‑110 (MQ=255)
                                                                                                                                       tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagtggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgc               <  2:4698679‑M2/150‑112 (MQ=255)
                                                                                                                                       tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgc               <  1:4950923‑M2/150‑112 (MQ=255)
                                                                                                                                       tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgc               <  1:5019474‑M2/150‑112 (MQ=255)
                                                                                                                                         ccACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctc             <  2:3145951‑M2/150‑114 (MQ=255)
                                                                                                                                         ccACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccgatgcaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctc             <  2:914177‑M2/150‑114 (MQ=255)
                                                                                                                                         CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt             <  2:5079722‑M1/150‑1 (MQ=255)
                                                                                                                                         CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt             <  1:5008525‑M1/150‑1 (MQ=255)
                                                                                                                                          cacaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcc            <  2:4543458‑M2/150‑115 (MQ=255)
                                                                                                                                          cacaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcc            <  2:340173‑M2/150‑115 (MQ=255)
                                                                                                                                          CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGGTGTCCGGCCCCATCAGGTACAGGTACACCGGCATg            <  1:760424‑M1/150‑1 (MQ=255)
                                                                                                                                           acaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcca           <  2:4430828‑M2/150‑116 (MQ=255)
                                                                                                                                           acaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcca           <  2:4430877‑M2/150‑116 (MQ=255)
                                                                                                                                           acaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcca           <  2:2949155‑M2/150‑116 (MQ=255)
                                                                                                                                           ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa           <  1:1998524‑M1/150‑1 (MQ=255)
                                                                                                                                           ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa           <  1:1998512‑M1/150‑1 (MQ=255)
                                                                                                                                           ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa           <  1:5103459‑M1/150‑1 (MQ=255)
                                                                                                                                           ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa           <  1:5103522‑M1/150‑1 (MQ=255)
                                                                                                                                           ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa           <  1:161857‑M1/150‑1 (MQ=255)
                                                                                                                                           ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa           <  2:3374847‑M1/150‑1 (MQ=255)
                                                                                                                                              GGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAtct        <  1:2931500‑M1/150‑1 (MQ=255)
                                                                                                                                                 tCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCtt     >  1:4166923/1‑150 (MQ=255)
                                                                                                                                                 tCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCtt     >  2:1441170/1‑150 (MQ=255)
                                                                                                                                                    ccGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGCACACCGGCATGATCTCTTCgg  >  1:2442109/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGG  >  NC_002947/1967077‑1967374

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