Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,346 T→G 8.1% T13P (ACG→CCG)  trmE ← GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,3460TG 8.1% 237.5 / 11.6 87T13P (ACG→CCG) trmEGTPase
Reads supporting (aligned to +/- strand):  ref base T (33/47);  new base G (4/3);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCT  >  NC_002947/6198‑6495
                                                                                                                                                    |                                                                                                                                                      
tCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTg                                                                                                                                                       <  2:1775683/150‑1 (MQ=255)
  gTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGCGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGc                                                                                                                                                     <  1:116527/150‑1 (MQ=255)
     cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGAt                                                                                                                                                  >  1:1424381/1‑150 (MQ=255)
      aGCACAAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGAGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATg                                                                                                                                                 >  1:847764/1‑150 (MQ=255)
       gCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACATCGCCACGGCCCAGGGCCGGGGCGGTgt                                                                                                                                                <  2:511573/150‑2 (MQ=255)
           cAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCTCCGTCCCTGGGCCGTGGCGATggcggc                                                                                                                                            >  2:1259474/1‑150 (MQ=255)
           cAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATggcggc                                                                                                                                            >  2:1623429/1‑150 (MQ=255)
                 ttCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGt                                                                                                                                      <  2:1697460/150‑1 (MQ=255)
                  tCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGtt                                                                                                                                     <  1:1692025/150‑1 (MQ=255)
                     tcgtcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCa                                                                                                                                  >  1:314994/1‑150 (MQ=255)
                      cgtcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGTCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCAc                                                                                                                                 <  1:1325531/150‑1 (MQ=255)
                      cgtcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTTACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCAc                                                                                                                                 <  2:1183193/150‑1 (MQ=255)
                      cgtcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCAc                                                                                                                                 <  2:1866523/150‑1 (MQ=255)
                        tcgCGGAATGGACCGTAATGGGCATGCCGCGGGGCCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGc                                                                                                                               >  1:910508/1‑150 (MQ=255)
                           cgGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                                            <  2:1748644/150‑1 (MQ=255)
                               aTGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAgtgt                                                                                                                        >  1:388817/1‑150 (MQ=255)
                                tGGACCGTAATGGGCATGCCGCGGGGTGAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                                       >  1:1211313/1‑150 (MQ=255)
                                tGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                                       <  2:94669/150‑1 (MQ=255)
                                 ggACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                                      <  2:1642299/150‑1 (MQ=255)
                                 ggACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                                      <  2:54097/150‑1 (MQ=255)
                                   aCCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                                       >  2:1442502/1‑147 (MQ=255)
                                   aCCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                                       <  1:1442502/147‑1 (MQ=255)
                                     cGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATTGCGGCGATGGTTTCACGCACAGTGTTCATGc                                                                                                                  <  1:154085/150‑1 (MQ=255)
                                       tAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                                <  1:894473/150‑1 (MQ=255)
                                           gggCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaat                                                                                                            <  2:936050/150‑2 (MQ=255)
                                           gggCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTCGTCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                                            <  1:1550971/150‑1 (MQ=255)
                                             gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                                          <  1:23820/150‑1 (MQ=255)
                                              cATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGcc                                                                                                         <  1:1902827/150‑1 (MQ=255)
                                                 gCCGCGGGGTAAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCctct                                                                                                      <  2:1211313/150‑1 (MQ=255)
                                                  ccGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTa                                                                                                     <  1:1623429/150‑1 (MQ=255)
                                                         gTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACaaa                                                                                              >  2:609812/1‑150 (MQ=255)
                                                           aaGGGTGCGGCCGGTAGTCAACAGGCCGGCCTGGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTAGAAAGCCTTTACGACaaaga                                                                                            <  2:335984/149‑3 (MQ=255)
                                                            aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGCTGGAGGCGGTGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAAc                                                                                           <  1:800009/150‑1 (MQ=255)
                                                                      cGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGc                                                                                 >  2:1884773/1‑150 (MQ=255)
                                                                        gTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaa                                                                               >  1:1659241/1‑150 (MQ=255)
                                                                         tGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGGGGCggt                                                                                                                                                  <  2:892168/82‑3 (MQ=255)
                                                                         tGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGGGGCggt                                                                                                                                                  >  1:892168/1‑80 (MQ=255)
                                                                               aCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                                        >  1:848259/1‑150 (MQ=255)
                                                                                cAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGccc                                                                       >  1:1805725/1‑150 (MQ=255)
                                                                                     cggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTa                                                                  <  1:1426360/150‑1 (MQ=255)
                                                                                       gccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAgg                                                               >  1:1235318/1‑150 (MQ=255)
                                                                                         cTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAgggg                                                             <  1:246342/150‑1 (MQ=255)
                                                                                          ttGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAggggg                                                            >  1:1112131/1‑150 (MQ=255)
                                                                                          ttGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAggggg                                                            >  1:1395841/1‑150 (MQ=255)
                                                                                              cAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtt                                                        <  1:676784/150‑1 (MQ=255)
                                                                                              cAGCCAAAGGGCCGGACAAACTGACAATTCCCACACCGCCCCGGCCCTGGGCAGTGGCAATGGCGGCGCTGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGAAACAACCGACAGATAGCAACAACGCCCCAATAGGGGGCGtt                                                        >  2:411637/1‑150 (MQ=255)
                                                                                               aGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGttt                                                       <  2:1395841/150‑1 (MQ=255)
                                                                                               aGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCCTAGTGGGGCGttt                                                       <  2:1741131/150‑1 (MQ=255)
                                                                                                   agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTa                                                   <  2:910508/150‑1 (MQ=255)
                                                                                                   agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTa                                                   <  2:1491795/150‑1 (MQ=255)
                                                                                                      gggCCGGACAACCTGACAATTCCCACACCGCCCCGGCCGCGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTc                                                <  1:1259474/150‑1 (MQ=255)
                                                                                                       ggCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCATCTAAGAAAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTca                                               <  2:847764/150‑1 (MQ=255)
                                                                                                        gCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTcac                                              >  1:1579160/1‑150 (MQ=255)
                                                                                                             aCAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCt                                         <  2:1235318/150‑1 (MQ=255)
                                                                                                             aCAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCt                                         <  1:1884773/150‑1 (MQ=255)
                                                                                                                 ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                                >  1:201310/1‑76 (MQ=255)
                                                                                                                 ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                                <  2:201310/76‑1 (MQ=255)
                                                                                                                  cTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                                            >  1:275810/1‑79 (MQ=255)
                                                                                                                  cTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                                            <  2:275810/79‑1 (MQ=255)
                                                                                                                  cTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGt                                    >  2:670649/1‑150 (MQ=255)
                                                                                                                  cTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGt                                    >  1:268969/1‑150 (MQ=255)
                                                                                                                      gatctCCCACCCCGCCCCGGCCCTGGGACGTGGCTATGTCGGCGGTGTTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTa                                <  2:207732/146‑1 (MQ=255)
                                                                                                                      cAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTa                                <  2:1805725/150‑1 (MQ=255)
                                                                                                                       aaTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAg                               >  2:714425/1‑150 (MQ=255)
                                                                                                                       aaTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGAAAGCAACACCGGGTGCACGCAGATAATTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTaa                               >  2:65895/1‑149 (MQ=255)
                                                                                                                          tCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtttt                                                      >  2:1049844/1‑124 (MQ=255)
                                                                                                                          tCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtttt                                                      <  1:1049844/124‑1 (MQ=255)
                                                                                                                          tCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                            <  2:848259/150‑1 (MQ=255)
                                                                                                                            ccACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGt                                                         >  1:269428/1‑119 (MQ=255)
                                                                                                                            ccACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGt                                                         <  2:269428/119‑1 (MQ=255)
                                                                                                                            ccACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAgagatc                          >  1:207732/1‑145 (MQ=255)
                                                                                                                             cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAggg                                                              <  2:715581/113‑1 (MQ=255)
                                                                                                                             cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAggg                                                              >  1:715581/1‑113 (MQ=255)
                                                                                                                               catcGCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                            >  2:819051/4‑145 (MQ=255)
                                                                                                                               catcGCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                            <  1:819051/142‑1 (MQ=255)
                                                                                                                               cacCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCaa                       <  2:1659241/150‑1 (MQ=255)
                                                                                                                                acCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAgg                                                               <  1:1263162/109‑1 (MQ=255)
                                                                                                                                acCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAgg                                                               >  2:1263162/1‑109 (MQ=255)
                                                                                                                                acCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCAAGTCAag                      >  1:27264/1‑150 (MQ=255)
                                                                                                                                 tcGCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc                     >  1:527200/2‑150 (MQ=255)
                                                                                                                                  cGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                                <  1:817866/59‑1 (MQ=255)
                                                                                                                                  cGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                                >  2:817866/1‑59 (MQ=255)
                                                                                                                                    ccacGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagc                  >  2:1172101/4‑150 (MQ=255)
                                                                                                                                     cccGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGCCAGAAACCAAAAACGCCCCAATAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagca                 >  2:998828/1‑150 (MQ=255)
                                                                                                                                      ccGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcag                <  1:971459/150‑1 (MQ=255)
                                                                                                                                             tGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGTCCGTAGAACTGTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttttt         <  1:1377094/150‑1 (MQ=255)
                                                                                                                                                   gTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAgc   <  2:388817/150‑1 (MQ=255)
                                                                                                                                                    tGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAgct  <  1:884099/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCT  >  NC_002947/6198‑6495

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