Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,878,536	G→A	Y528Y (TAC→TAT)	gyrB ←	DNA gyrase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,878,536	0	G	A	100.0%	212.2 / NA	58	Y528Y (TAC→TAT)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/46);  total (12/46)

ACGGTCTGCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCT  >  NC_000913/3878243‑3878822                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                               aCGGTCTGCTCATCGGAAATGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCCCTTTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  2:257211/300‑1 (MQ=255) aCGGTCTGCTCATCATAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATGGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  2:430536/300‑1 (MQ=255) aCGGTCTGCTCATCAGAAGGGTCAGCTTCCGTCAACTTCGGCTGATAGATAAGCTCTTTCAGTATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGTCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  2:119324/300‑1 (MQ=255) aCGGTCTGCTCATCAGAACAGTCAGCTTCCGTCACCGTCGGCTGATAGAGAAGCTCTTTCAGCATTGCTTTCGGATAGCGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCATATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTACGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  2:119354/300‑1 (MQ=255) aCGGTCTGCTCATCAGAAAGGTCAGCTTCCTTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  1:56072/300‑1 (MQ=255) aCGGTCTGCTCATCAGAAAGGTCAGCTTCCGTTAACGTCTGCTGATAGATAGGCTCTGTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  2:350241/300‑1 (MQ=255) aCGGTCTGCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  1:23236/300‑1 (MQ=255) aCGGTCTGCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  2:277937/300‑1 (MQ=255)                                  aaCGTCGGCTGATAGATAAGCTCTTTCAGCGGTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCTTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCTCCAGCCAATGCCGGTGCACTGGCGGTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATATAGATCTGGTACTGATCCATCGCGTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  2:277768/267‑1 (MQ=255)                                               agataAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGt                                                                                                                                                                                                                                              <  2:379949/298‑1 (MQ=255)                                                    aGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTANTGATCCATCGCTTCGTCGTCTTTAATGTACTGNTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAAcagc                                                                                                                                                                                                                                       <  2:479301/300‑1 (MQ=255)                                                        ctTTCAGCATTGCTTTCGGATAATGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGTAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCAATTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGACAGAAGAAGGTCaaaagaagc                                                                                                                                                                                                                                    >  2:573144/1‑292 (MQ=255)                                                             aGCAGTGCTTTCGGATAACGACGCTCCCGCGCGATTGATCATTTTCTGTGTCGCGTTGTACTCAGATACCAGTTCCTCTAACGCGTCGCCAGCCAATGCCGGTGCACTTGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACg                                                                                                                                                                                                                               <  2:210876/300‑1 (MQ=255)                                                                aTTGCTTTCGGATAACGACACTCC‑ATACGATTGATCATTTTCTGCTTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAAtg                                                                                                                                                                                                                           <  2:546013/300‑1 (MQ=255)                                                                                               ttGATCATTTTCTGCGTCGCGTTGTACTCAGATACTAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAACGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCTTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                          <  2:119335/207‑1 (MQ=255)                                                                                                                 gcgTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGNTAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATaa                                                                                                                                                                           >  1:469117/1‑300 (MQ=255)                                                                                                                   gTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTAATGATGATGATGCTGTGATAAc                                                                                                                                                                          >  2:567406/1‑299 (MQ=255)                                                                                                                    ttGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGc                                                                                                                                                                        <  1:428389/300‑1 (MQ=255)                                                                                                                                   aGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGtt                                                                                                                                                          <  1:550802/299‑1 (MQ=255)                                                                                                                                      ttCTCTAACTCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTAc                                                                                                                                                      <  2:428824/300‑1 (MQ=255)                                                                                                                                      ttCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTAc                                                                                                                                                      <  1:391809/300‑1 (MQ=255)                                                                                                                                       tctcTAATGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTGAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACTTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGGCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACt                                                                                                                                                     <  1:445228/300‑1 (MQ=255)                                                                                                                                       tctcTAACGCTTCGCCATCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACt                                                                                                                                                     <  1:445083/300‑1 (MQ=255)                                                                                                                                        ctctAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTc                                                                                                                                                    >  2:430447/1‑300 (MQ=255)                                                                                                                                                 ttCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTTTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGAc                                                                                                                                            <  1:401834/299‑1 (MQ=255)                                                                                                                                                 ttCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGAc                                                                                                                                            <  1:440874/299‑1 (MQ=255)                                                                                                                                                       agccaATGCCTGTGCACTGGCTTTGGTGGGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTCCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCTGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATAcc                                                                                                                                     <  2:269548/300‑1 (MQ=255)                                                                                                                                                            aTGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTNGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGc                                                                                                                                <  1:131232/300‑1 (MQ=255)                                                                                                                                                            aTGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGc                                                                                                                                >  1:29729/1‑300 (MQ=255)                                                                                                                                                            aTGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTACGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGc                                                                                                                                <  1:522679/300‑1 (MQ=255)                                                                                                                                                            aTGCAGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGc                                                                                                                                >  1:29504/1‑300 (MQ=255)                                                                                                                                                                        tGGCGTTGGTGTGCAGCTTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGa                                                                                                                    <  1:163057/300‑1 (MQ=255)                                                                                                                                                                        tGGCGTTGGTGTGAAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCTTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATATAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATAATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGa                                                                                                                    <  1:163027/300‑1 (MQ=255)                                                                                                                                                                                      aatttggcgcCGTCCAGCGAGAGATGGATCTGGTACTGATCCACCGCTTCTTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCg                                                                                                              <  2:625075/286‑1 (MQ=255)                                                                                                                                                                                       gCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGNTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAgcgc                                                                                                                         >  1:459509/1‑280 (MQ=255)                                                                                                                                                                                       gCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACtt                                                                                                     >  1:94944/1‑300 (MQ=255)                                                                                                                                                                                       gCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCANCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAgcgc                                                                                                                         <  2:459509/280‑1 (MQ=255)                                                                                                                                                                                         gTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTcc                                                                                                   <  1:193156/300‑1 (MQ=255)                                                                                                                                                                                              gcCGTCCAGCTCGAGAGAGATCTGGTACTGATCCATCTCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTgag                                                                                               <  2:38185/299‑1 (MQ=255)                                                                                                                                                                                                 gTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCTTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGagaag                                                                                            <  2:534890/299‑1 (MQ=255)                                                                                                                                                                                                 gTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCTGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGagaag                                                                                            <  2:534873/299‑1 (MQ=255)                                                                                                                                                                                                 gTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCNTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGagaag                                                                                            <  2:469117/299‑1 (MQ=255)                                                                                                                                                                                                 gTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGagaag                                                                                            <  1:394407/299‑1 (MQ=255)                                                                                                                                                                                                 gTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGagaag                                                                                            <  2:534706/299‑1 (MQ=255)                                                                                                                                                                                                     aGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGc                                                                                       <  1:58264/300‑1 (MQ=255)                                                                                                                                                                                                     aGCGCGAGAGAAATCTGTTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACGGTGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGTTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGTTCATGATGATGATGCGGTGATAACACAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTCCCTTAGAAGAGAGc                                                                                       <  2:167116/300‑1 (MQ=255)                                                                                                                                                                                                             agagATCTGTTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATAGAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATc                                                                               <  2:320935/300‑1 (MQ=255)                                                                                                                                                                                                                 aTCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAg                                                                           <  1:643319/300‑1 (MQ=255)                                                                                                                                                                                                                 aTCTGGTACGGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTTAGAAGAGAGCATCTTATCGAAg                                                                           <  1:351088/300‑1 (MQ=255)                                                                                                                                                                                                                                      tcgtcgTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAgg                                                      <  1:567406/300‑1 (MQ=255)                                                                                                                                                                                                                                        gtcgtcTTTAATGTACTGTTCCTACTTGCCTTTCTTCGCTTTGTACAGCGGCGGCTGAGCGACATAGACATTACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTGCGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCAGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAACGCGCTTTCTCGACGTTGAGGAt                                                    <  2:1369/300‑1 (MQ=255)                                                                                                                                                                                                                                        gtcgccTTTAATGTAATGTTCCTACTTGCCGTTCTTCACTTTTTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCGACGATTGCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATCATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAATAGAGCATCTTATCGAAGCGCACTTTCTCGACGTTTAGGAt                                                    <  2:902/300‑1 (MQ=255)                                                                                                                                                                                                                                           gtcTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGAt                                                    <  2:1388/297‑1 (MQ=255)                                                                                                                                                                                                                                                         gTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGNATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCa                                   >  2:465980/1‑300 (MQ=255)                                                                                                                                                                                                                                                                      tttCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGtt                      <  1:30837/300‑1 (MQ=255)                                                                                                                                                                                                                                                                             aCTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAACGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGc                 >  2:63188/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                ttGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTa            >  1:440220/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                          gcggcTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTAACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCAACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCt  >  2:82443/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                               ACGGTCTGCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCC‑ATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCT  >  NC_000913/3878243‑3878822

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