Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,967,224 A→G 9.4% intergenic (+51/+68) 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,2240AG 9.4% 239.5 / 11.9 87intergenic (+51/+68)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (33/43);  new base G (5/3);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTT  >  NC_002947/1967075‑1967371
                                                                                                                                                     |                                                                                                                                                   
tGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCa                                                                                                                                                     >  2:4134521/1‑150 (MQ=255)
  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:1281101‑M2/1‑143 (MQ=255)
                                 gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgg                                                                                                                    >  2:1834668‑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                                                                                                                          <  2:2447975‑M1/64‑1 (MQ=255)
                                  tgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                          >  1:2447975‑M1/80‑143 (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                                                                                                              >  1:612301‑M2/1‑137 (MQ=255)
                                       tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                              >  2:1199323‑M2/1‑137 (MQ=255)
                                       tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                              >  2:2124840‑M2/1‑137 (MQ=255)
                                        gcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcggg                                                                                                             >  2:4145372‑M2/1‑136 (MQ=255)
                                           tgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgggcaa                                                                                                          >  1:175508‑M2/1‑133 (MQ=255)
                                           tgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgggcaa                                                                                                          >  2:757647‑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                              <  2:720505‑M1/136‑1 (MQ=255)
                                                                                                                                     CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATc                              >  1:720505‑M1/1‑136 (MQ=255)
                                                                                                                                     CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                <  1:2665218‑M1/150‑1 (MQ=255)
                                                                                                                                      CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc               <  1:4049659‑M1/150‑1 (MQ=255)
                                                                                                                                      CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc               <  1:4979621‑M1/150‑1 (MQ=255)
                                                                                                                                      CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc               <  1:3153429‑M1/150‑1 (MQ=255)
                                                                                                                                      CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc               <  2:112904‑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:5019474‑M2/150‑112 (MQ=255)
                                                                                                                                         tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgc            <  1:4950923‑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        <  2:3374847‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa        <  1:161857‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa        <  1:1998524‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa        <  1:1998512‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa        <  1:5103522‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa        <  1:5103459‑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)
                                                                                                                                                     |                                                                                                                                                   
TGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTT  >  NC_002947/1967075‑1967371

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: