Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,967,224 A→G 5.9% 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 5.9% 367.0 / 8.3 121intergenic (+51/+68)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (53/57);  new base G (2/5);  total (58/63)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 122 total aligned reads displayed.

TGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATC  >  NC_002947/1967075‑1967367
                                                                                                                                                     |                                                                                                                                               
tGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCa                                                                                                                                                 <  1:2669942/150‑1 (MQ=255)
 gCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCAc                                                                                                                                                >  1:677424/1‑150 (MQ=255)
    cgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGATCACCGc                                                                                                                                             <  1:2947464‑M1/41‑1 (MQ=255)
    cgcagaacaacccgctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGATCACCGc                                                                                                                                             <  1:2947438‑M1/41‑1 (MQ=255)
            gACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCcaca                                                                                                                                     >  2:1749011/1‑150 (MQ=255)
             aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                    <  1:5373755/150‑1 (MQ=255)
             aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                    <  2:4894649/150‑1 (MQ=255)
             aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                    <  2:2802188/150‑1 (MQ=255)
               tGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAATCCttgagg                                                                                                                                  <  2:1183292/150‑6 (MQ=255)
                gCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGa                                                                                                                                 <  1:4827935/150‑1 (MQ=255)
                   tgttcctgttccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTg                                                                                                                              <  1:4515583/150‑1 (MQ=255)
                   tgttcctgttccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTg                                                                                                                              <  1:4704831/150‑1 (MQ=255)
                   tgttcctgttccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTg                                                                                                                              <  1:4515495/150‑1 (MQ=255)
                    gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                                                                                   >  2:3219105/1‑144 (MQ=255)
                    gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                                                                                   <  1:491378/144‑1 (MQ=255)
                    gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                                                                                   <  1:3219105/144‑1 (MQ=255)
                    gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                                                                                   >  2:491378/1‑144 (MQ=255)
                        gaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCTCACCGCAACCCTCAAGGATTGCGccc                                                                                                                         >  1:2679525‑M1/90‑150 (MQ=255)
                         tgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                                                                                        <  1:2159497/150‑1 (MQ=255)
                         aagagatcatgcctgtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTACCACAGGCAGTCCGCAACCCTCAAGGATTGCGCcct                                                                                                                        >  2:5388137‑M1/89‑150 (MQ=255)
                         aagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                                                                                        >  2:3729967‑M1/89‑150 (MQ=255)
                           gagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaaccttcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                      >  1:9469‑M1/87‑150 (MQ=255)
                           gagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                      >  1:5240922‑M1/87‑150 (MQ=255)
                             ccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcct                                                                                                                    >  1:4141020‑M2/1‑147 (MQ=255)
                               tGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                                                                                  >  1:3160343‑M2/1‑145 (MQ=255)
                               tGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                                                                                  <  2:129472‑M2/150‑6 (MQ=255)
                                 gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgg                                                                                                                >  1:1672148‑M2/1‑143 (MQ=255)
                                 atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                >  1:663805‑M1/81‑150 (MQ=255)
                                 atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                >  1:479094‑M1/81‑150 (MQ=255)
                                 atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                >  1:4610377‑M1/81‑150 (MQ=255)
                                 atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                >  2:3228606‑M1/81‑150 (MQ=255)
                                 atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                >  2:4018635‑M1/81‑150 (MQ=255)
                                   gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                                                                              >  1:582655‑M2/1‑141 (MQ=255)
                                   gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                                                                              >  2:677424‑M2/1‑141 (MQ=255)
                                   gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                                                                              >  2:3140530‑M2/1‑141 (MQ=255)
                                   gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                                                                              >  1:1006619‑M2/1‑141 (MQ=255)
                                   gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                                                                              >  1:1423016‑M2/1‑141 (MQ=255)
                                   gccggtgtacctgtacctgatggggccggacagcccaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGCCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                              >  2:3037968‑M1/79‑150 (MQ=255)
                                    cgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCAGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggag                                                                                                             >  2:2528114‑M2/1‑140 (MQ=255)
                                    ccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCTGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAg                                                                                                             >  1:3139835‑M1/78‑150 (MQ=255)
                                      ggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCg                                                                                                           <  1:852422‑M1/75‑1 (MQ=255)
                                       tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                          >  1:129472‑M2/1‑137 (MQ=255)
                                       tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                          >  1:510943‑M2/1‑137 (MQ=255)
                                       gtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttctacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                                          >  1:4503072‑M1/75‑150 (MQ=255)
                                       gtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaatacacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                                          >  1:1718426‑M1/75‑150 (MQ=255)
                                        gcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcggg                                                                                                         >  1:4305268‑M2/1‑136 (MQ=255)
                                        gcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcggg                                                                                                         >  1:4305437‑M2/1‑136 (MQ=255)
                                           tgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgggcaa                                                                                                      >  1:1292428‑M2/1‑133 (MQ=255)
                                           acctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaatacacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCaa                                                                                                      >  2:2906420‑M1/71‑150 (MQ=255)
                                             ctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGc                                                                                                    >  2:2947464‑M1/69‑150 (MQ=255)
                                             ctgtaccggatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGGATCACCGCAACCCTCACGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGc                                                                                                    >  2:2947438‑M1/69‑150 (MQ=255)
                                                 tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                                                                              >  1:2768206‑M2/1‑127 (MQ=255)
                                                 tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                                                                              <  2:2768206‑M2/136‑10 (MQ=255)
                                                   gAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGGGTGGCCGCTACCCCCCGCGCCCCGCAACCCAAAAGGCATTCGCCCCGCCGGGGGGGGAGGGCCAGCCcgcgc                                                                                              >  2:3808223/1‑149 (MQ=255)
                                                   gAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCGGccggtgggtgcgggcaagcccgcga                                                                                              <  2:3808260‑M2/150‑26 (MQ=255)
                                                      atggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                                                                                        <  1:2103438‑M1/62‑1 (MQ=255)
                                                      atggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                                                                                        >  2:2103438‑M1/60‑121 (MQ=255)
                                                               gacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                                          >  2:1318704‑M1/51‑126 (MQ=255)
                                                                               cggcggGGGGTGCGGGCTGAGCGCTGGCCGGGCCCGGCCTCTGCGCGGGCTGGCCCGCCCCCACCGGAGCCCCGCACCCCCCACGGATTGCGCCCGGCCTGGGGGAGCGGGCACGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGtt                                                                  <  2:2978498/150‑1 (MQ=255)
                                                                                  cGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCATTCCCCCACCCCCAGGATGAACCCCACCCAAAAGGCTTTGGCCCGGCCCGGGGGGGCGGGCAAGCCCGCGGACAGGCCCGTGGCGGGGTATTGGGCGGTGac                                                               >  2:4577478/1‑149 (MQ=255)
                                                                                             aacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                                          >  1:2428254‑M1/21‑96 (MQ=255)
                                                                                             aacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                                          <  2:2428254‑M1/76‑1 (MQ=255)
                                                                                                tGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                                                                                  >  1:1853041‑M2/1‑80 (MQ=255)
                                                                                                tGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                                                                                  <  2:1853041‑M2/85‑6 (MQ=255)
                                                                                                                    CCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCACGGCTGGCGCCCTGCCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGcccc                             >  2:4243916‑M1/1‑150 (MQ=255)
                                                                                                                    CCTCTTCGCGGGCTGGCCCGCTCCCACAGGAACCCCGCAACCCACAAGGATTGCGCCCGGCCTGTGGGTGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGcccc                             >  2:4243972‑M1/1‑150 (MQ=255)
                                                                                                                        tgCGCGGGCGGGCCCGCCCCCACCGGATCACCGCACCCCACACGGATGGCGCCCTGCCTGGGGGCGCGGGCACGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATca                         >  1:1749011‑M1/3‑150 (MQ=255)
                                                                                                                             GGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTGCAGGTTTATTGGGCGTTGACCGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtac                    <  2:5348254‑M1/150‑1 (MQ=255)
                                                                                                                                    CCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAggtacaggtaca             <  1:3271320‑M1/150‑1 (MQ=255)
                                                                                                                                    CCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAggtacaggtaca             <  2:1750531‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac            <  2:5422080‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac            <  1:4503062‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac            <  1:914787‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac            <  1:233662‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac            <  1:690283‑M1/150‑1 (MQ=255)
                                                                                                                                      CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc           <  1:1680917‑M1/150‑1 (MQ=255)
                                                                                                                                         tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgc        <  1:5172459‑M2/150‑112 (MQ=255)
                                                                                                                                         tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgc        <  1:4102123‑M2/150‑112 (MQ=255)
                                                                                                                                         tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgccctt                         <  1:4992872‑M2/133‑95 (MQ=255)
                                                                                                                                         tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgccctt                         >  2:4992872‑M2/1‑39 (MQ=255)
                                                                                                                                         TCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCCGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGc        <  2:2994100‑M1/150‑1 (MQ=255)
                                                                                                                                         TCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGc        <  2:1490525‑M1/150‑1 (MQ=255)
                                                                                                                                          cccACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgct       <  2:4974738‑M2/150‑113 (MQ=255)
                                                                                                                                           ccACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccgatgcaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctc      <  1:1280934‑M2/150‑114 (MQ=255)
                                                                                                                                           ccACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccgatgcaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctc      <  1:1280915‑M2/150‑114 (MQ=255)
                                                                                                                                           CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt      <  2:2225367‑M1/150‑1 (MQ=255)
                                                                                                                                           CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt      <  2:4256749‑M1/150‑1 (MQ=255)
                                                                                                                                           CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt      <  1:3066284‑M1/150‑1 (MQ=255)
                                                                                                                                            cacaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagct                   <  2:4009602‑M2/136‑101 (MQ=255)
                                                                                                                                            cacaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagct                   >  1:4009602‑M2/1‑36 (MQ=255)
                                                                                                                                            CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATg     <  1:1418527‑M1/150‑1 (MQ=255)
                                                                                                                                            CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATg     <  1:2021543‑M1/150‑1 (MQ=255)
                                                                                                                                             acaGGATCACCGCAACCCTCAAGGATTGCGGTGATcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcca    <  1:3123692‑M2/150‑116 (MQ=255)
                                                                                                                                             acaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcca    <  1:1794367‑M2/150‑116 (MQ=255)
                                                                                                                                             acaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcca    <  1:2020968‑M2/150‑116 (MQ=255)
                                                                                                                                             ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAgg                                                                               <  1:3002362‑M1/75‑1 (MQ=255)
                                                                                                                                             ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa    <  1:954308‑M1/150‑1 (MQ=255)
                                                                                                                                             ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa    <  1:954306‑M1/150‑1 (MQ=255)
                                                                                                                                              ccGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAt   >  2:4225428‑M1/3‑150 (MQ=255)
                                                                                                                                               AGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAtc  <  1:1531636‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATC  >  NC_002947/1967075‑1967367

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