Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	883,396	C→T	intergenic (‑8/‑277)	ybjG ← / → mdfA	undecaprenyl pyrophosphate phosphatase/multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	883,396	0	C	T	94.8%	188.5 / ‑0.4	58	intergenic (‑8/‑277)	ybjG/mdfA	undecaprenyl pyrophosphate phosphatase/multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (24/31);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

CGGAAAAAGATGTCCCATCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAAC  >  NC_000913/883152‑883641                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       cGGAAAAAGATGTCCCATCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGc                                                                                                                                                                                                                                                   >  1:444016/1‑250 (MQ=255)   gAAAAAGATGTCCCATCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATACAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGCGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCaa                                                                                                                                                                                                                                                 >  2:157567/1‑250 (MQ=255)      aaaGATGTCCCATCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTg                                                                                                                                                                                                                                             >  2:214031/1‑251 (MQ=255)          atgccccAGCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTa                                                                                                                                                                                                                                         <  2:117660/247‑1 (MQ=255)            gTCCCATCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTGCCAGCATAACTTCCTGACGCAAAGTGATTaaa                                                                                                                                                                                                                                       <  1:681261/251‑1 (MQ=255)                 ctCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGCTTTTTATCACCAGTTGCCGTAGTGCTGTAAGCCCCCACAACCAAAGTACAACGGCCAGCAACGACACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTACTGAGTGCGAATGGAGAGAGATTCAAATTTTCCAGCATAACTCCCTGACGCAAAGTGTTTAAAAgggg                                                                                                                                                                                                                                  >  2:22772/2‑250 (MQ=255)                   cgcccAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTGAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGAATAACTTCCTGACGCAAAGTGATTAAAAGGGAGc                                                                                                                                                                                                                                <  1:368930/248‑1 (MQ=255)                     tCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGcc                                                                                                                                                                                                                               >  2:501190/1‑250 (MQ=255)                      ccAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCa                                                                                                                                                                                                                              >  1:381145/1‑250 (MQ=255)                      ccAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCa                                                                                                                                                                                                                              <  1:605471/250‑1 (MQ=255)                        aGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑Ta                                                                                                                                                                                                                          <  1:748577/251‑1 (MQ=255)                            cacaAACAGGCTGACCGCCAGCGCGATAGGGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGGACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTATCAATAAAAATCGCCAACGAGATCATCCACGGGGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGAGTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAgg                                                                                                                                                                                                                      <  1:360194/251‑1 (MQ=255)                              aaacaaggctacctgcccgcccgATAGCGATTGTTATCCCCAATGGCCGTTTGGCTGGAAGCCCCCCAAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCAGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTGGCTGACGCAAAGTGATTAAAAGGGTGCCAA‑GACAGGCa                                                                                                                                                                                                                    <  2:335640/230‑1 (MQ=255)                                       tGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAGGTCGTTg                                                                                                                                                                                                            >  1:570108/1‑250 (MQ=255)                                        gACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTg                                                                                                                                                                                                            <  1:214031/249‑1 (MQ=255)                                        gACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGAAGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTg                                                                                                                                                                                                            <  1:543913/249‑1 (MQ=255)                                         aCCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTgaga                                                                                                                                                                                                         <  1:719579/251‑1 (MQ=255)                                          ccGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGaa                                                                                                                                                                                                        <  2:337457/251‑1 (MQ=255)                                                   cgctAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAAAGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGATAATAAGAGAGAATAGAGAGAGATTCAAATTTTACAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCAAAGTACAG‑CACGTCATTGAGAATAAAGTGc                                                                                                                                                                                                >  2:555283/4‑250 (MQ=255)                                                                          tGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGCGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCAt                                                                                                                                                                        >  1:283960/1‑251 (MQ=255)                                                                             cGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTACAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGt                                                                                                                                                                     >  2:499437/1‑251 (MQ=255)                                                                                    cTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTAttt                                                                                                                                                              >  2:634456/1‑251 (MQ=255)                                                                                       tAAGCCCCCACACCCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTACTAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATCACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGCGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTccc                                                                                                                                                           >  2:52545/1‑251 (MQ=255)                                                                                       tAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGGGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTccc                                                                                                                                                           >  1:389345/1‑251 (MQ=255)                                                                                       tAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTccc                                                                                                                                                           >  1:80030/1‑251 (MQ=255)                                                                                                         gTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCTAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACg                                                                                                                                         <  1:16614/251‑1 (MQ=255)                                                                                                                          ggCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACgg                                                                                                                        >  2:374321/1‑251 (MQ=255)                                                                                                                                     aaTCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGTAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTTAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTGCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTtaa                                                                                                             <  1:342289/251‑1 (MQ=255)                                                                                                                                              tttAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATCGTAACGGGCCTAGGTTAAATAAAACtt                                                                                                    >  2:564535/1‑251 (MQ=255)                                                                                                                                               ttAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAAc                                                                                                      <  1:669980/248‑1 (MQ=255)                                                                                                                                               gtAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACt                                                                                                     <  1:374321/248‑1 (MQ=255)                                                                                                                                                 aGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTaa                                                                                                  <  1:287525/250‑1 (MQ=255)                                                                                                                                                  gCAATAAAAATCGCCAAAGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTaaag                                                                                                <  2:570108/251‑1 (MQ=255)                                                                                                                                                    aaTAAAAAGCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGGGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTaaagaa                                                                                              <  2:389345/251‑1 (MQ=255)                                                                                                                                                     aTAAAAATCGCCAACGAGATCAGCCACGGAGCCGAGTCGGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTGTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTaaagaaa                                                                                             <  1:706060/251‑1 (MQ=255)                                                                                                                                                          aaTCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGt                                                                                         >  2:51372/1‑250 (MQ=255)                                                                                                                                                             cGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGc                                                                                      <  1:378464/250‑1 (MQ=255)                                                                                                                                                                    gagaTCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAC‑TACAGGCACGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGttt                                                                                                                                            >  2:390463/1‑189 (MQ=255)                                                                                                                                                                    gagaTCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAATTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCGGttt                                                                                                                                            <  1:390463/189‑1 (MQ=255)                                                                                                                                                                            ccACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATg                                                                      >  1:173368/1‑251 (MQ=255)                                                                                                                                                                                      gAGTCTGGCGGCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAAtgtg                                                            <  2:80030/251‑1 (MQ=255)                                                                                                                                                                                             gCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAAcc                                                     <  1:738505/251‑1 (MQ=255)                                                                                                                                                                                                   cgTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGtt                                               >  1:473169/1‑251 (MQ=255)                                                                                                                                                                                                    gTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGttt                                              <  1:157567/251‑1 (MQ=255)                                                                                                                                                                                                                 aaTAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGTCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAAtt                                 >  1:454793/1‑251 (MQ=255)                                                                                                                                                                                                                  aTAGAGAGAGATGCAAATTTTCCAGCATAACTTCCTGACGCAAAGAGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTc                                <  2:588335/251‑1 (MQ=255)                                                                                                                                                                                                                    agagagagATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCg                               <  1:110760/250‑1 (MQ=255)                                                                                                                                                                                                                        agagaTTCAAATTTTCCAGCAGAACTACCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTtc                          <  2:242616/251‑1 (MQ=255)                                                                                                                                                                                                                         gagaTTCAAATTTTCCAGCAGAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAA‑TTCAGACAAGTCGTTGAGAATAAAGAGCAGGTTAAACGGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATATAACTTAAAGGAAGCGTGGCTATACTCGTAATAGTGTAAGAATGTGCTTAACAGGGGTTTAAGCTACAATATTAGCTTtcg                         >  2:226122/1‑250 (MQ=255)                                                                                                                                                                                                                                   ttttCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGAGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTgc               >  1:306085/1‑251 (MQ=255)                                                                                                                                                                                                                                    tttCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTgc               <  1:634456/250‑1 (MQ=255)                                                                                                                                                                                                                                     ttCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTTGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTgcg              <  2:381145/250‑1 (MQ=255)                                                                                                                                                                                                                                     ttCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTgcg              <  1:612937/250‑1 (MQ=255)                                                                                                                                                                                                                                     ttCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGGTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTgcg              <  1:52545/250‑1 (MQ=255)                                                                                                                                                                                                                                      tCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAATAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTgcgc             <  1:499437/250‑1 (MQ=255)                                                                                                                                                                                                                                       ccAGCATAACTTACTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAACAATGTGCTTAACCGTGGTTTCAGCTACAAAATTAGCTTTCTCGTTAGCTGCGCtt           >  2:195992/1‑251 (MQ=255)                                                                                                                                                                                                                                                cTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAAc  >  1:533124/1‑251 (MQ=255)                                                                                                                                                                                                                                                cTTCCTGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAAc  >  1:218584/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CGGAAAAAGATGTCCCATCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAA‑TACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAAAC  >  NC_000913/883152‑883641

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑