Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 842,254 G→T intergenic (‑198/‑78) mcbA ← / → rlmF DUF1471 domain‑containing protein McbA/23S rRNA m(6)A1618 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913842,2540GT98.4% 226.0 / NA 64intergenic (‑198/‑78)mcbA/rlmFDUF1471 domain‑containing protein McbA/23S rRNA m(6)A1618 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (41/22);  minor base C (1/0);  total (42/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGCGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCC  >  NC_000913/841959‑842546
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    
agagacaggTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATACTGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGaa                                                                                                                                                                                                                                                                                                  <  2:558913/292‑1 (MQ=255)
 gACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAAt                                                                                                                                                                                                                                                                                                 <  1:42397/300‑1 (MQ=255)
        gTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTCAGGTTTTGTGGAATGATAAAc                                                                                                                                                                                                                                                                                          >  1:530087/1‑300 (MQ=255)
        gTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGaactgtctat                                                                                                                                                                                                                                                                                          >  1:568819/1‑292 (MQ=255)
        gTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATacac                                                                                                                                                                                                                                                                                          >  2:648324/1‑297 (MQ=255)
        gTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAAc                                                                                                                                                                                                                                                                                          >  2:283550/1‑300 (MQ=255)
          ccAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGc                                                                                                                                                                                                                                                                                        >  2:397117/1‑300 (MQ=255)
            aGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGCTTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCATGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAg                                                                                                                                                                                                                                                                                      >  1:322551/1‑300 (MQ=255)
                  cTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGGGGTACTTTTTCATAATGACATCCTCATAAATCCAACGGGTAGTGCCTGATACGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATNTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCGTAAACCATAGCCAGGTTAAGGTTTTGTGGAATGGTAAACGCAGCCAggt                                                                                                                                                                                                                                                                                >  2:179032/1‑298 (MQ=255)
                  cTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGat                                                                                                                                                                                                                                                                                >  1:717794/1‑300 (MQ=255)
                  cTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGat                                                                                                                                                                                                                                                                                >  2:82119/1‑300 (MQ=255)
                    cATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCAATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGatat                                                                                                                                                                                                                                                                              >  1:107778/1‑300 (MQ=255)
                    cATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAAGGGCCATAAACCATAGCCAGGTGAAGGTTTTGTGGAAGGATagaggaggggagaggt                                                                                                                                                                                                                                                                              <  1:405834/300‑16 (MQ=255)
                                 gCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTcgcg                                                                                                                                                                                                                                                                 >  1:138902/1‑300 (MQ=255)
                                      aGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCAATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATg                                                                                                                                                                                                                                                            >  1:606270/1‑300 (MQ=255)
                                                 agaTTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTa                                                                                                                                                                                                                                                 >  1:114197/1‑300 (MQ=255)
                                                    ttCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTAcgc                                                                                                                                                                                                                                              <  1:315978/300‑1 (MQ=255)
                                                        gCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCgct                                                                                                                                                                                                                                          >  1:628467/1‑300 (MQ=255)
                                                        gCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCAGGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCgct                                                                                                                                                                                                                                          >  2:437780/1‑300 (MQ=255)
                                                                 ggTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGTTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAgg                                                                                                                                                                                                                                 <  1:325761/300‑1 (MQ=255)
                                                                 ggTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAgg                                                                                                                                                                                                                                 <  1:325566/300‑1 (MQ=255)
                                                                 ggTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAgg                                                                                                                                                                                                                                 <  2:397222/300‑1 (MQ=255)
                                                                  gTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGga                                                                                                                                                                                                                                >  1:14014/1‑300 (MQ=255)
                                                                       aaTCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATNCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGACATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATg                                                                                                                                                                                                                           >  2:158918/1‑300 (MQ=255)
                                                                       aaTCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGNATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATCATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATg                                                                                                                                                                                                                           >  1:523231/1‑300 (MQ=255)
                                                                       aaTCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATg                                                                                                                                                                                                                           >  1:99408/1‑300 (MQ=255)
                                                                       aaTCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATg                                                                                                                                                                                                                           >  2:395388/1‑300 (MQ=255)
                                                                       aaTCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGACAAGTGTAGGTCGCCTGCGCGGTGCAAATGACGCAAAGAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATNCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACACTGTCGCGTTATGAGGGTATATCACTCCGCTGCCTAAAGGAGCATg                                                                                                                                                                                                                           >  2:158690/1‑300 (MQ=255)
                                                                         tCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAGTGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGCTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCt                                                                                                                                                                                                                         >  2:222501/1‑300 (MQ=255)
                                                                          cAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTGTGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTa                                                                                                                                                                                                                        >  1:134831/1‑300 (MQ=255)
                                                                           agtagtGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACANAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTAt                                                                                                                                                                                                                       <  2:630702/300‑1 (MQ=255)
                                                                                                     gACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTGAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATc                                                                                                                                                                                             <  2:2158/300‑1 (MQ=255)
                                                                                                           cTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGAACGCGCGATGTAAATGACACAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCACTGACATGCATAATGCTGATTTTTATGCTGCATTGCCATAACCCATAGCCATGTTCCGGTTTTGTGGCATGATAAACGCAGCCAGATATCTACCCCGTCGCGTTATGCGGGTATTTTATGCCGCTGCTTAAAGGAGAATGATATGTCCGCCCAGAACACGGGGTTGCATCCGCGCa                                                                                                                                                                                       >  2:287478/1‑300 (MQ=255)
                                                                                                                  aTACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAANAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTNAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGtca                                                                                                                                                                                >  2:673408/1‑300 (MQ=255)
                                                                                                                      aaCGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATcac                                                                                                                                                                            <  1:386897/300‑1 (MQ=255)
                                                                                                                          ggTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGcc                                                                                                                                                                        >  1:714404/1‑300 (MQ=255)
                                                                                                                          ggTAGTGCCTGATAAGTGTAGGTCGCATGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATACACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGcc                                                                                                                                                                        >  1:400007/1‑300 (MQ=255)
                                                                                                                                 ccTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTaaga                                                                                                                                                                 >  2:125114/1‑297 (MQ=255)
                                                                                                                                        aGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCNGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGcc                                                                                                                                                          >  1:482332/1‑300 (MQ=255)
                                                                                                                                        aGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCAATGTTAAGGTTTTGTGGAATGATAAACGCAGCAAGATATTTACCATGCCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTAGcc                                                                                                                                                          >  1:74346/1‑300 (MQ=255)
                                                                                                                                         gtgtAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATNAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCAAAGCCGCTACGAtcccgaac                                                                                                                                                         >  2:158203/1‑294 (MQ=255)
                                                                                                                                         gtgtAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCNAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCa                                                                                                                                                         >  2:536889/1‑300 (MQ=255)
                                                                                                                                         gtgtAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCa                                                                                                                                                         >  1:419152/1‑300 (MQ=255)
                                                                                                                                         gtgtAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTAAAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCAATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCCCAGCCGCTACGATCTCGCCa                                                                                                                                                         >  1:293787/1‑300 (MQ=255)
                                                                                                                                                                 aTGACGCAAAAAATTGACATAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATNAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAAc                                                                                                                                 <  1:519075/300‑1 (MQ=255)
                                                                                                                                                                 aTGACGCAAAAAATTGACATAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATACTATGTCCGCCCAGAAACCGGAGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAAc                                                                                                                                 <  1:149679/300‑1 (MQ=255)
                                                                                                                                                                                    gagaTGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCt                                                                                                              <  2:701994/300‑1 (MQ=255)
                                                                                                                                                                                    gagaTGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCt                                                                                                              >  1:237774/1‑300 (MQ=255)
                                                                                                                                                                                         gtgCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTacac                                                                                                         <  1:21393/300‑1 (MQ=255)
                                                                                                                                                                                          tgCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACAcc                                                                                                        <  2:702918/300‑1 (MQ=255)
                                                                                                                                                                                            cTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACccg                                                                                                      >  2:214552/1‑300 (MQ=255)
                                                                                                                                                                                                           atatatTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCg                                                                                       >  1:724006/1‑300 (MQ=255)
                                                                                                                                                                                                                 ttttACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATTTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACt                                                                                 <  1:93352/300‑1 (MQ=255)
                                                                                                                                                                                                                 ttttACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACt                                                                                 <  1:705483/300‑1 (MQ=255)
                                                                                                                                                                                                                                         gACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAgg                                                         <  1:82338/300‑1 (MQ=255)
                                                                                                                                                                                                                                              gCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTc                                                    >  1:124897/1‑300 (MQ=255)
                                                                                                                                                                                                                                                      cTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGc                                            <  2:349922/300‑1 (MQ=255)
                                                                                                                                                                                                                                                      cTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGc                                            <  1:285348/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gCCATAAACCATAGCCATGTTAAGGTTTTGTNGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACNGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGCGTTGCTGGCCCATTTTTAc                         >  2:506918/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                         gCCATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGCGTTGCTGGCCCATTTTTAc                         >  1:468017/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                           cATAAACCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGCGTTGCTGGCCCATTTTTACGc                       <  1:660304/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aaCCATAGCCATGTTAAGGTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGCGTTGCTGGCCCATTTTTACGCCGTa                   >  1:580189/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                         aTGTTCAGGTTTTGCGGAATGATAAACGCAGCCAGGTATTTACCCTGTCGCGTTATGTGGGTATCTTACGCAGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGAGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATTCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAACAAGGCGTTGCTTGCCCATTCTTACGCAGTatcgcaccggg         >  2:301392/1‑290 (MQ=255)
                                                                                                                                                                                                                                                                                                ggTTTTGTGGAATGATAAACGCAGCCAGATATTTACCCTGGCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAATCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGCGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATccc  <  2:617251/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    
TGACCGCTGTCCAGGTTACTCATCGGTTGCGCGGCATAAGCCGTGAGAGAGATTCCGCTCAGGACGGTGGCAATCAGTAGTGTGAGGCACTTTTTCATAATGACATCCTCATAAATACAACGGGTAGTGCCTGATAAGTGTAGGTCGCCTGCGCGATGCAAATGACGCAAAAAATTGACAGAGATGTGCTTATTTTTTGAACGATATATTTTTACAAAATAAGACAAATCAATGACATGCATAATGCTGATTTTTATGATGAATTGCCATAAACCATAGCCATGTTAAGGTTTTGGGGAATGATAAACGCAGCCAGATATTTACCCTGTCGCGTTATGCGGGTATCTTACGCCGCTGCTTAAAGGAGAATGCTATGTCCGCCCAGAAACCGGGGTTGCATCCGCGCAACCGTCATCACAGCCGCTACGATCTCGCCACGCTTTGTCAGGTCAATCCTGAACTCAGGCAATTCCTCACGCTTACACCCGCCGGGGAGCAAAGCGTAGACTTTGCCAATCCGCTGGCGGTGAAGGCGCTCAATAAGGCGTTGCTGGCCCATTTTTACGCCGTAGCGAACTGGGATATCCC  >  NC_000913/841959‑842546

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