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%	141.8 / NA	40	Y528Y (TAC→TAT)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/21);  total (19/21)

TCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGG  >  NC_000913/3878301‑3878770                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           tCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTTGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCTCTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACcgcg                                                                                                                                                                                                                              <  1:265364/250‑1 (MQ=255)   aGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTc                                                                                                                                                                                                                           <  1:137718/251‑1 (MQ=255)    gCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCa                                                                                                                                                                                                                          <  1:228653/251‑1 (MQ=255)                gATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGc                                                                                                                                                                                                              >  1:40527/1‑251 (MQ=255)                      gacgCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTg                                                                                                                                                                                                         >  1:85632/1‑250 (MQ=255)                       acgCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGAc                                                                                                                                                                                                       >  1:84512/1‑251 (MQ=255)                                aCGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATagaaga                                                                                                                                                                                              <  1:7253/251‑1 (MQ=255)                                aCGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATagaag                                                                                                                                                                                               <  2:25902/250‑1 (MQ=255)                                                tGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGt                                                                                                                                                                              >  1:132315/1‑251 (MQ=255)                                                 gCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTGCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCGGCTTGTCTTTCTTCTCTTTGTACAGCGGCGGCTGAGCGATATAGACGTGAGCGCGTTCTACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGt                                                                                                                                                                              <  1:129445/250‑1 (MQ=255)                                                   gTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACg                                                                                                                                                                           <  1:1293/251‑1 (MQ=255)                                                                 aGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGt                                                                                                                                                             >  2:28313/1‑251 (MQ=255)                                                                         gTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTcc                                                                                                                                                     >  1:203929/1‑251 (MQ=255)                                                                            tctcTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGc                                                                                                                                                   >  1:268477/1‑250 (MQ=255)                                                                             ctctAACGCTTCGCCAGACAATGCCGGTGCACTGGCGTTGGAGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCAATCGCTTCGTAGTCTTTAATGTACTGTTCCTGCTTGCCTTTATTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAGTGTGCGAGCCGGCGACGTCCGCAt                                                                                                                                                 >  1:106735/1‑251 (MQ=255)                                                                                 aaCGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATc                                                                                                                                                <  2:227084/248‑1 (MQ=255)                                                                                 aaCGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGt                                                                                                                                             <  2:30252/251‑1 (MQ=255)                                                                                          ccagccaATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCatgatgat                                                                                                                                    >  2:217161/1‑251 (MQ=255)                                                                                                        tGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAAc                                                                                                                      >  1:258137/1‑251 (MQ=255)                                                                                                              ggCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAg                                                                                                                  <  1:41290/249‑1 (MQ=255)                                                                                                              ggCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGtt                                                                                                                >  1:247973/1‑251 (MQ=255)                                                                                                                    ggTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTcc                                                                                                           <  2:219012/250‑1 (MQ=255)                                                                                                                           aGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTa                                                                                                   <  2:268477/251‑1 (MQ=255)                                                                                                                           aGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTa                                                                                                   <  1:98013/251‑1 (MQ=255)                                                                                                                                                   gagaGCTGGTACTGAGCCATCGCTTCGGCGGCTTTAATGTACTGTGCCTGCTTTCCTTTCTTCTCTTTGTACAGCTGCGGCTGAGCGATATAGACGTGACCGCGTTCACCGATTGCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAgac                                                                           <  2:206111/251‑3 (MQ=255)                                                                                                                                                                gctccttCGCTTCGTCGTCGTTAATGTACTGTGCCTGCTGGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATTCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACGGCCAAGCGCGGTGa                                                                <  2:215045/243‑1 (MQ=255)                                                                                                                                                                 aTCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGCGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGGGATa                                                              >  2:122505/1‑250 (MQ=255)                                                                                                                                                                   ccATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGc                                                           <  1:142230/251‑1 (MQ=255)                                                                                                                                                                               cgtcTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTcc                                               >  1:217955/1‑251 (MQ=255)                                                                                                                                                                                  cTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTg                                             >  1:145694/1‑250 (MQ=255)                                                                                                                                                                                        tGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGagaag                                        >  1:189071/1‑249 (MQ=255)                                                                                                                                                                                         gTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGagaag                                        >  1:147018/1‑248 (MQ=255)                                                                                                                                                                                              gTTCCTGATTGCCTTTTTTCACTTTGTACAGCAGCGGCTGAGCGATATAGCCGTGACCGCGCTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGAAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCGGTGATAACGCAGTTTGTCCGGGTTGACCTCGTGACGACCACCACCACAGCACCGCGCGGCGATAATCGTCGCCCCTTCCTgcgagacgcgcatc                                >  2:33882/1‑238 (MQ=255)                                                                                                                                                                                                 ccTGCTTGCCTTTCTTCACTTGGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTa                             <  2:204836/251‑1 (MQ=255)                                                                                                                                                                                                         ccTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAgcg                     <  1:232133/251‑1 (MQ=255)                                                                                                                                                                                                            ttcttcACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAgcgc                    <  1:122207/249‑1 (MQ=255)                                                                                                                                                                                                                tcACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTtc              >  2:224846/1‑251 (MQ=255)                                                                                                                                                                                                                     ttGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCg           <  2:85632/249‑1 (MQ=255)                                                                                                                                                                                                                            gcggcgGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCAGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAgg  <  1:2912/251‑1 (MQ=255)                                                                                                                                                                                                                             cggcggCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGCAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGGTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGCCGTTGAgg  >  2:63876/1‑250 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           TCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGG  >  NC_000913/3878301‑3878770

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