Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 838,720 A→G 8.7% intergenic (‑20/+39) rplY ← / ← prs 50S ribosomal protein L25/general stress protein Ctc/ribose‑phosphate pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947838,7200AG 8.7% 457.5 / 24.5 161intergenic (‑20/+39)rplY/prs50S ribosomal protein L25/general stress protein Ctc/ribose‑phosphate pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (69/75);  new base G (9/5);  total (79/82)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

Only 100 of 166 total aligned reads displayed.

GACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGC  >  NC_002947/838571‑838869
                                                                                                                                                     |                                                                                                                                                     
 aCGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAg                                                                                                                                                      >  1:3428425/1‑150 (MQ=255)
 aCGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAg                                                                                                                                                      <  1:4960625/150‑1 (MQ=255)
   gCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGt                                                                                                                                                    <  1:3179701/150‑1 (MQ=255)
   gCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGt                                                                                                                                                    >  2:2180694/1‑150 (MQ=255)
     aGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTg                                                                                                                                                  >  1:3494375/1‑150 (MQ=255)
     aGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTg                                                                                                                                                  <  1:738490/150‑1 (MQ=255)
            ggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACCTGTGTTggggcggg                                                                                                                                           <  2:4790592/150‑1 (MQ=255)
            ggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcggg                                                                                                                                           <  1:3913729/150‑1 (MQ=255)
            ggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGACGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcggg                                                                                                                                           <  2:4633915/150‑1 (MQ=255)
              cTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcggggc                                                                                                                                         <  1:2999051/150‑1 (MQ=255)
                 gCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTGTGTTTCTCCAAAGTAGCCCCCCGAGGGCGGGTGCGCCCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCACCACCGGTTGGGGCGGGGCGCt                                                                                                                                      >  1:4709741/1‑150 (MQ=255)
                    ccTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACCTGTGTTggggcggggctttttt                                                                                                                                   <  1:1721773/150‑7 (MQ=39)
                     cTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACCTGTGTTggggcggggcttttttg                                                                                                                                  <  1:3471449/150‑8 (MQ=39)
                     cTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTcac                                                                                                                                  <  1:2180694/150‑1 (MQ=255)
                      tttCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACCTGTGTTggggcggggcttttt                                                                                                                                    >  2:3381331/1‑142 (MQ=39)
                      tttCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTcaca                                                                                                                                 <  1:4809069/150‑1 (MQ=255)
                      tttCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTcaca                                                                                                                                 >  1:600408/1‑150 (MQ=255)
                       ttCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACAt                                                                                                                                >  1:679927/1‑150 (MQ=255)
                         ccccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCa                                                                                                                              >  1:688663/1‑150 (MQ=255)
                          cccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCACCACAGGTTGGGGCGGGGCGCTTCACATCaa                                                                                                                             >  1:4670661/1‑150 (MQ=255)
                          cccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCaa                                                                                                                             >  2:1452311/1‑150 (MQ=255)
                          cccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCaa                                                                                                                             >  1:3116486/1‑150 (MQ=255)
                          cccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCaa                                                                                                                             >  1:2473169/1‑150 (MQ=255)
                           ccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAac                                                                                                                            >  2:2006236/1‑150 (MQ=255)
                           ccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATAAac                                                                                                                            >  1:3537289/1‑150 (MQ=255)
                               gTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACgt                                                                                                                        <  2:699954/150‑1 (MQ=255)
                                   gTACGCGCTTGGGCGTTCAGAGTGATTTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTcc                                                                                                                    <  1:1085306/150‑1 (MQ=255)
                                      cgcgCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGc                                                                                                                 <  1:4648635/150‑1 (MQ=255)
                                       gcgcTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACAt                                                                                                                                >  1:2971595/1‑134 (MQ=255)
                                               gCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGaa                                                                                                        >  1:3396335/1‑150 (MQ=255)
                                                  ttCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGAGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACac                                                                                                     <  2:4364031/150‑2 (MQ=255)
                                                   tCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATc                                                                                                    >  1:1193618/1‑150 (MQ=255)
                                                   tCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATc                                                                                                    <  1:2042718/150‑1 (MQ=255)
                                                      gagTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATcgcg                                                                                                 >  2:1700898/1‑150 (MQ=255)
                                                      gagTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATcgcg                                                                                                 >  1:2405799/1‑150 (MQ=255)
                                                       agTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCgcgc                                                                                                >  1:440011/1‑150 (MQ=255)
                                                        gTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCt                                                                                               >  2:2324973/1‑150 (MQ=255)
                                                          gAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGa                                                                                             <  2:3904410/150‑1 (MQ=255)
                                                          gAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGa                                                                                             <  2:748000/150‑1 (MQ=255)
                                                           aaTTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGAt                                                                                            >  2:2855781/1‑150 (MQ=255)
                                                            aTTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATc                                                                                           >  2:4982387/1‑150 (MQ=255)
                                                                aGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACgt                                                                                                                        >  2:1869689/1‑117 (MQ=255)
                                                                aGTCATGTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGAtt                                                                                       <  1:3866133/150‑1 (MQ=255)
                                                                 gTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGAttc                                                                                      >  1:1817393/1‑150 (MQ=255)
                                                                  tCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGaa                                                                                                        <  2:2926162/131‑1 (MQ=255)
                                                                  tCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGAttct                                                                                     >  2:400644/1‑150 (MQ=255)
                                                                    aTTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGAttcttc                                                                                   >  1:1420136/1‑150 (MQ=255)
                                                                     ttttGTTTCTCCAAAGTAGCCCACCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGCCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACAACGCGCTGATCGATTCTTCg                                                                                  >  2:2630494/1‑150 (MQ=255)
                                                                      tttGTTTCTCCAAATTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGTGTATTTCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCAGCTTAGCGGAACATCGCGCTGATCTATTCTTCGt                                                                                 <  2:1954615/150‑1 (MQ=255)
                                                                              tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCgcgc                                                                                                >  1:931647/1‑127 (MQ=255)
                                                                              tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATg                                                                         >  2:2740438/1‑150 (MQ=255)
                                                                              tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATg                                                                         >  1:302968/1‑150 (MQ=255)
                                                                              tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATg                                                                         >  1:1640716/1‑150 (MQ=255)
                                                                               ccAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACCGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgc                                                                        <  1:39581/150‑1 (MQ=255)
                                                                               ccAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgc                                                                        <  1:2766913/150‑1 (MQ=255)
                                                                                 aaaGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgcgg                                                                      >  2:414044/1‑150 (MQ=255)
                                                                                  aaGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgcggc                                                                     <  2:1590116/150‑1 (MQ=255)
                                                                                  aaGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgcggc                                                                     <  2:600408/150‑1 (MQ=255)
                                                                                   aGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAATGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgcggcg                                                                    >  1:4422761/1‑150 (MQ=255)
                                                                                    gTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGt                                                                   <  2:1193618/150‑1 (MQ=255)
                                                                                    gTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGt                                                                   >  2:1408493/1‑150 (MQ=255)
                                                                                    gTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGt                                                                   >  1:3267647/1‑150 (MQ=255)
                                                                                      aGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTAc                                                                 <  1:402483/150‑1 (MQ=255)
                                                                                        ccccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCg                                                               <  2:4612018/150‑1 (MQ=255)
                                                                                          ccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCt                                                             >  1:5313980/1‑150 (MQ=255)
                                                                                          ccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCt                                                             >  2:1825477/1‑150 (MQ=255)
                                                                                               aGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGc                                                        <  2:4938032/150‑1 (MQ=255)
                                                                                               aGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGc                                                        <  2:1808813/150‑1 (MQ=255)
                                                                                               aGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGGTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGc                                                        <  2:3625436/150‑1 (MQ=255)
                                                                                                     tttGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAAc                                                  <  1:1217580/150‑1 (MQ=255)
                                                                                                     tttGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAAc                                                  >  2:4762253/1‑150 (MQ=255)
                                                                                                       tGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCg                                                >  2:4831365/1‑150 (MQ=255)
                                                                                                          gACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTg                                             <  1:3347760/150‑1 (MQ=255)
                                                                                                           aCCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                            <  2:1640716/150‑1 (MQ=255)
                                                                                                           aCCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                            <  1:5294239/150‑1 (MQ=255)
                                                                                                            ccAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCg                                           >  2:2563005/1‑150 (MQ=255)
                                                                                                              aGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCg                                                                                  <  2:310526/109‑1 (MQ=255)
                                                                                                                       ggggATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTg                                             >  1:1113050/1‑137 (MQ=255)
                                                                                                                       ggggATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTg                                             <  2:1113050/137‑1 (MQ=255)
                                                                                                                         ggATGGCATAATAGCCCCGCCCCAACACATGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGc                              <  2:2082400/150‑1 (MQ=255)
                                                                                                                         ggATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGc                              <  2:3428425/150‑1 (MQ=255)
                                                                                                                         ggATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGc                              <  1:622793/150‑1 (MQ=255)
                                                                                                                            tGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGa                           <  2:3147627/150‑1 (MQ=255)
                                                                                                                              gCAAAAAAGCCCCGCCCCAACAACGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATa                         <  1:4228397/150‑1 (MQ=255)
                                                                                                                                aaaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACg                       >  2:1530552/1‑150 (MQ=255)
                                                                                                                                   aaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTc                    >  1:1939170/1‑150 (MQ=255)
                                                                                                                                   aaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTc                    <  2:1939170/150‑1 (MQ=255)
                                                                                                                                       ccccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAg                >  1:1757251/1‑150 (MQ=255)
                                                                                                                                        cccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAgg               >  1:858159/1‑150 (MQ=255)
                                                                                                                                        cccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAgg               >  2:2550009/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc              >  2:4198619/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc              >  2:3155448/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGTGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc              >  2:5101734/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc              >  1:2843827/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACACAGGTTGGGGCGGGGCGCATCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc              <  2:3912777/150‑1 (MQ=255)
                                                                                                                                          cgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCt             <  1:140882/150‑1 (MQ=255)
                                                                                                                                           gccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCtt            <  2:839020/150‑1 (MQ=255)
                                                                                                                                              ccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGag         <  1:2740414/150‑1 (MQ=255)
                                                                                                                                               cAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGagc        >  2:1907626/1‑150 (MQ=255)
                                                                                                                                                   acaGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCg    >  1:1138485/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGC  >  NC_002947/838571‑838869

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