Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 838,718 A→C 8.7% intergenic (‑18/+41) 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,7180AC 8.7% 463.0 / 28.2 161intergenic (‑18/+41)rplY/prs50S ribosomal protein L25/general stress protein Ctc/ribose‑phosphate pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (70/77);  new base C (9/5);  total (79/82)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 165 total aligned reads displayed.

GACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCG  >  NC_002947/838571‑838867
                                                                                                                                                   |                                                                                                                                                     
 aCGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAg                                                                                                                                                    >  1:3428425/1‑150 (MQ=255)
 aCGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAg                                                                                                                                                    <  1:4960625/150‑1 (MQ=255)
 aCGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAg                                                                                                                                                    >  1:4612018/1‑150 (MQ=255)
       gcggcggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTggg                                                                                                                                              >  1:5129840/1‑150 (MQ=255)
            ggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACCTGTGTTggggcggg                                                                                                                                         <  2:1612767/150‑1 (MQ=255)
            ggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACCTGTGTTggggcggg                                                                                                                                         <  2:4790592/150‑1 (MQ=255)
            ggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGACGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcggg                                                                                                                                         <  2:4633915/150‑1 (MQ=255)
             gcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcgggg                                                                                                                                        <  2:2629260/150‑1 (MQ=255)
              cTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcggggc                                                                                                                                       <  1:2999051/150‑1 (MQ=255)
              cTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcggggc                                                                                                                                       <  1:4746080/150‑1 (MQ=255)
                 gCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTGTGTTTCTCCAAAGTAGCCCCCCGAGGGCGGGTGCGCCCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCACCACCGGTTGGGGCGGGGCGCt                                                                                                                                    >  1:4709741/1‑150 (MQ=255)
                  cACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCtt                                                                                                                                   >  2:140882/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                                                                                                                                  <  1:3381331/147‑6 (MQ=39)
                      tttCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACCTGTGTTggggcggggcttttt                                                                                                                                  >  2:3381331/1‑142 (MQ=39)
                      tttCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTcaca                                                                                                                               >  1:600408/1‑150 (MQ=255)
                         ccccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCa                                                                                                                            >  1:688663/1‑150 (MQ=255)
                          cccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCACCACAGGTTGGGGCGGGGCGCTTCACATCaa                                                                                                                           >  1:4670661/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)
                            cAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAaca                                                                                                                         <  1:64626/150‑1 (MQ=255)
                               gTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACgt                                                                                                                      <  2:699954/150‑1 (MQ=255)
                                   gTACGCGCTTGGGCGTTCAGAGTGATTTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTcc                                                                                                                  <  1:1085306/150‑1 (MQ=255)
                                       gcgcTTGGGCGTTCAGAGTTAATTCAGTCTTTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCTTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCt                                                                                                              <  1:747/150‑1 (MQ=255)
                                       gcgcTTGGGCGTTCAGAGTGAATTCATTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCt                                                                                                              <  1:730/150‑1 (MQ=255)
                                       gcgcTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACAt                                                                                                                              >  1:2971595/1‑134 (MQ=255)
                                                  ttCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACAt                                                                                                   >  2:2282099/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                                                                                               >  1:2405799/1‑150 (MQ=255)
                                                      gagTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATcgcg                                                                                               >  2:1700898/1‑150 (MQ=255)
                                                       agTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCgcgc                                                                                              >  1:440011/1‑150 (MQ=255)
                                                        gTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCt                                                                                             >  2:2324973/1‑150 (MQ=255)
                                                          gAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGa                                                                                           <  2:748000/150‑1 (MQ=255)
                                                          gAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGa                                                                                           <  2:3904410/150‑1 (MQ=255)
                                                            aTTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATc                                                                                         >  2:4982387/1‑150 (MQ=255)
                                                            aTTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATc                                                                                         >  1:97466/1‑150 (MQ=255)
                                                                aGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACgt                                                                                                                      >  2:1869689/1‑117 (MQ=255)
                                                                aGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGAtt                                                                                     <  1:3901265/150‑1 (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)
                                                                    aTTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGAttcttc                                                                                 >  1:3153092/1‑150 (MQ=255)
                                                                     ttttGTTTCTCCAAAGTAGCCCACCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGCCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACAACGCGCTGATCGATTCTTCg                                                                                >  2:2630494/1‑150 (MQ=255)
                                                                              tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCgcgc                                                                                              >  1:931647/1‑127 (MQ=255)
                                                                              tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATg                                                                       >  1:1640716/1‑150 (MQ=255)
                                                                              tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATg                                                                       >  1:302968/1‑150 (MQ=255)
                                                                              tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATg                                                                       >  2:2740438/1‑150 (MQ=255)
                                                                               ccAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACCGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgc                                                                      <  1:39581/150‑1 (MQ=255)
                                                                               ccAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgc                                                                      >  2:286748/1‑150 (MQ=255)
                                                                               ccAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgc                                                                      <  1:2766913/150‑1 (MQ=255)
                                                                               ccAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgc                                                                      <  1:3488066/150‑1 (MQ=255)
                                                                                cAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATgcg                                                                     >  1:4153109/1‑150 (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:1408493/1‑150 (MQ=255)
                                                                                    gTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGt                                                                 <  2:1193618/150‑1 (MQ=255)
                                                                                    gTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGt                                                                 >  1:3267647/1‑150 (MQ=255)
                                                                                      aGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTAc                                                               <  1:402483/150‑1 (MQ=255)
                                                                                        ccccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCg                                                             <  2:4612018/150‑1 (MQ=255)
                                                                                               aGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGc                                                      <  2:1808813/150‑1 (MQ=255)
                                                                                               aGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGc                                                      <  2:4938032/150‑1 (MQ=255)
                                                                                               aGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGGTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGc                                                      <  2:3625436/150‑1 (MQ=255)
                                                                                                     tttGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAAc                                                >  2:4762253/1‑150 (MQ=255)
                                                                                                     tttGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAAc                                                <  1:1217580/150‑1 (MQ=255)
                                                                                                       tGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCg                                              >  2:4831365/1‑150 (MQ=255)
                                                                                                          gACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTg                                           <  1:3347760/150‑1 (MQ=255)
                                                                                                              aGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCg                                                                                <  2:310526/109‑1 (MQ=255)
                                                                                                                    aaCGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAg                                 >  1:5311047/1‑150 (MQ=255)
                                                                                                                     aCGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGc                                >  2:1961258/1‑150 (MQ=255)
                                                                                                                       ggggATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTg                                           >  1:1113050/1‑137 (MQ=255)
                                                                                                                         ggATGGCATAATAGCCCCGCCCCAACACATGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGc                            <  2:2082400/150‑1 (MQ=255)
                                                                                                                               cAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATAc                      <  1:621201/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             >  2:2550009/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc            >  2:3155448/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc            >  2:4198619/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCGGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc            >  2:3155446/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGTGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc            >  2:5101734/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc            >  1:2843827/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc            >  2:3525503/1‑150 (MQ=255)
                                                                                                                                         ccgccccAACACAGGTTGGGGCGGGGCGCATCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGc            <  2:3912777/150‑1 (MQ=255)
                                                                                                                                          cgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCt           <  1:140882/150‑1 (MQ=255)
                                                                                                                                           gccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCtt          <  2:1895098/150‑1 (MQ=255)
                                                                                                                                           gccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCtt          <  2:839020/150‑1 (MQ=255)
                                                                                                                                              ccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGag       <  2:1817393/150‑1 (MQ=255)
                                                                                                                                              ccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGag       <  1:2740414/150‑1 (MQ=255)
                                                                                                                                                   acaGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCg  >  1:1138485/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCG  >  NC_002947/838571‑838867

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: