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	98.2%	165.7 / ‑4.0	55	Y528Y (TAC→TAT)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (28/25);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

CATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACC  >  NC_000913/3878305‑3878778                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    ccttgttttCGGATAACGACGCTCCATCCGATTGATCATTTTCTGCGTCGCGTTGTACTC‑AGAGCCCAGTTTCTTTACCGCGTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCaa                                                                                                                                                                                                                                  <  2:387774/245‑1 (MQ=255)     gCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATTTAGACGTGACCGCGTTCAACGAt                                                                                                                                                                                                                              >  2:104800/1‑251 (MQ=255)         tCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTcc                                                                                                                                                                                                                          >  1:255751/1‑251 (MQ=255)           ggATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCgg                                                                                                                                                                                                                        >  1:201632/1‑251 (MQ=255)           ggATAACGACGCTCCATACGATTGATCATTGTCTGCGTCGCGTTGTACTC‑AGATACCAGTTGCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCgg                                                                                                                                                                                                                        <  2:403514/251‑1 (MQ=255)                acgacgCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCt                                                                                                                                                                                                                   <  1:67205/251‑1 (MQ=255)                      cTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCTCGTTCAACGATTTCCGGCATCTGACGAt                                                                                                                                                                                                             >  1:315860/1‑251 (MQ=255)                             cGATTGATCATTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGTGCGGTATAGCAGTGACCCCATTCAACGATTTCCGGCATCTGTCGTTagaagcg                                                                                                                                                                                                      >  2:112553/1‑249 (MQ=255)                                 ggATCATTTTCGGCGTCGCGTTGT‑CTCCAGACACCAGTTTCTCTACCGCTTCCCCAGCCAATGCCTGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCGTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGt                                                                                                                                                                                                   <  2:359501/249‑1 (MQ=255)                                      aTTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAAc                                                                                                                                                                                               >  1:362888/1‑249 (MQ=255)                                      aTTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAAc                                                                                                                                                                                               >  1:46175/1‑249 (MQ=255)                                         ttCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGCGCGCTATAGACGTGACCGCGGTCACCGATTTCCGGCATCTGACGATAGAAGAAGGTCAAcagc                                                                                                                                                                                            >  2:108755/1‑249 (MQ=255)                                            tGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGt                                                                                                                                                                                       >  1:258550/1‑251 (MQ=255)                                                    gTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGACGAAGGTCAACAGCAGCGTACG‑AAtgt                                                                                                                                                                              >  1:299125/1‑251 (MQ=255)                                                      tGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGc                                                                                                                                                                            >  2:375640/1‑251 (MQ=255)                                                            c‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGACAACAGCAGCGTACG‑AATGTGCGAGCCg                                                                                                                                                                      >  2:158858/1‑251 (MQ=255)                                                                aTACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCg                                                                                                                                                                   >  1:81502/1‑251 (MQ=255)                                                                        ttCTCTAACGCTCCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGAGCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTGCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCg                                                                                                                                                            <  2:57182/250‑1 (MQ=255)                                                                            ctAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCAt                                                                                                                                                         <  1:103412/249‑1 (MQ=255)                                                                                    tCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCatg                                                                                                                                                 <  1:83091/249‑1 (MQ=255)                                                                                       ccagccaATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTTCTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCTCTTTGTACAGCGGCGGCTGAGCGATATATACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCatgatgat                                                                                                                                            <  2:392548/251‑1 (MQ=255)                                                                                       ccagccaAGGCCGGTGCACTGGCGTGGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTGCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCatgatgat                                                                                                                                            <  2:203998/251‑1 (MQ=255)                                                                                                gCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCTGCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGGGACCGCGTTCAACGATTGCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCtgtg                                                                                                                                   <  2:254851/251‑1 (MQ=255)                                                                                                  cGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGAACGCGTTCACCGATTTCCGGCCTCTGACGATAGAAGAAGGACAACAGCAGCGTACG‑AAAGTGCGAGCCGTCGACGTCCGCATCGGCCATGATGATGCAGCTGTGAt                                                                                                                                 >  2:386997/1‑251 (MQ=255)                                                                                                     tGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGCGCGATATAGACGTTACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCCGCGTACG‑CATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGCTGATGCtgtgtaaac                                                                                                                              >  2:15796/1‑246 (MQ=255)                                                                                                             ggTTGGTGTGCCGCGTTGCGCCGGCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCTCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTg                                                                                                                      <  2:315860/250‑1 (MQ=255)                                                                                                                 ggTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCg                                                                                                                  <  1:356100/251‑1 (MQ=255)                                                                                                                  gtgtGCAGCGTTGCGACGTCCATCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCgg                                                                                                                 <  1:289703/251‑1 (MQ=255)                                                                                                                        aGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGt                                                                                                            >  1:301536/1‑250 (MQ=255)                                                                                                                        aGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGt                                                                                                            >  1:128624/1‑250 (MQ=255)                                                                                                                                               agagATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGAGATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGc                                                                                    <  1:209623/251‑1 (MQ=255)                                                                                                                                               agagATCTGGGACTGATCCTTCGCTTCGGCGTCGTTAATGTCCTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAg                                                                                     <  2:102044/250‑1 (MQ=255)                                                                                                                                                       ggTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTTTTCACTTTGTACAGCGGCGGATGAGCGATATCGACGTGAACGCGTTCAACGATTTTCGGCATATTCAGCTAGAAGAAGGCCAACAGAAGCGTACG‑ACTGTGCGAGCCGTCGACGTCCGCATCGGTCCTGATGATGATGCTGCGATAACGCAGTTTGTCCGGGTTGTCAGCGTCACGACCGATACCACAGCCAcgcgcg                                                                            >  1:405930/1‑251 (MQ=255)                                                                                                                                                          aCTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGt                                                                          <  1:160499/250‑1 (MQ=255)                                                                                                                                                            tGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGCCCGATACCCCAGCCAAGCGGGGTGAt                                                                       >  2:194264/1‑251 (MQ=255)                                                                                                                                                                 cATCGTTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCg                                                                  <  2:255751/251‑1 (MQ=255)                                                                                                                                                                  aTCGCTTAGTCGGCTTTAATGTACTGTCCCTGCTTGCCTTTCTTCACTGTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGt                                                                 <  2:411956/251‑1 (MQ=255)                                                                                                                                                                     gCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGCGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCCGTTTGTCCGGGTTGTTCTCGTCACGACCGATACCACAGCCAAGCGCggcgaacagcaccgc                                                              >  2:72607/1‑238 (MQ=255)                                                                                                                                                                       ttCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCCAGCGCGGTGATAAGCGTCGccc                                                            >  1:186949/1‑250 (MQ=255)                                                                                                                                                                                 ttAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTg                                                     <  1:105918/248‑1 (MQ=255)                                                                                                                                                                                  tAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCCGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTgaga                                                  >  1:145173/1‑250 (MQ=255)                                                                                                                                                                                      gTACTGTTCCTGCTTGCCTTTCTTCTCTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCTGTTATAAGCGTCGCCACTTCCTGAGAagag                                              <  1:360250/250‑1 (MQ=255)                                                                                                                                                                                      gTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAagag                                              <  1:104800/250‑1 (MQ=255)                                                                                                                                                                                           gttccgttttGCCTGTCTTCACTTTGTCCAGCGGTGGCTGAGCGATATAGACGTTACCGCGTTCAACGATTTCCGGCATCTGACGAGAGAAGAAGGTCAACAGCAGCGTACG‑AATTTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTTTGATAACGCAGTTTGTCTTGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTTATAAGCGTCGTCACTTCCTGAGAAGAGAGCAt                                         <  2:230155/242‑1 (MQ=255)                                                                                                                                                                                              ccTGCTTGCCTTTCTTCTCTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTa                                     <  1:158858/251‑1 (MQ=255)                                                                                                                                                                                                    tGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCTTTCAACGATTGCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCTTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGTGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGaa                                <  1:209744/250‑1 (MQ=255)                                                                                                                                                                                                       cTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGCGAGCATCTTATCGAAgcgg                            >  1:41091/1‑250 (MQ=255)                                                                                                                                                                                                        tttCATCACTTTGTACAGCGTCGGCTGATCGATATATTCGTGACCTCGTTCCACGATTTCCGGCATCTGACGATAGAACAAGGGCAACAGCCACGTACGCAAT‑TGCGAGCCGTCGACGTCCGCCTCGTTCATGCTGATGAAGCTGTGATAACGCCGTTTGTCCGGGTTGTAATCGTCACTAACGGTACCACCGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTTGGAAGAGAACCTCATCTCGAAgcgcc                           >  2:321035/1‑250 (MQ=255)                                                                                                                                                                                                            ttcACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCttt                       >  1:34562/1‑251 (MQ=255)                                                                                                                                                                                                                  ttGTACAGCGGCGGCTGAGCGAGATAGGCGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCg                   <  2:34562/249‑1 (MQ=255)                                                                                                                                                                                                                       cAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCACGCGCGGTGATAAGCGACGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCCCTTACTCGCCGTTg             >  2:338106/1‑250 (MQ=255)                                                                                                                                                                                                                             cggcTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGAttt      >  1:305329/1‑251 (MQ=255)                                                                                                                                                                                                                             cggcTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCCGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGAttt      >  1:142794/1‑251 (MQ=255)                                                                                                                                                                                                                              ggcTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGAtttt     >  1:155194/1‑251 (MQ=255)                                                                                                                                                                                                                                 tGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTAcc  >  2:188916/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    CATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTC‑AGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACG‑AATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACC  >  NC_000913/3878305‑3878778

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