Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,339,161	T→C	D87G (GAC→GGC)	gyrA ←	DNA gyrase subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,339,161	0	T	C	98.1%	186.2 / ‑3.9	52	D87G (GAC→GGC)	gyrA	DNA gyrase subunit A
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (24/27);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

CACCAGCAGGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGTTGACCGGTGTAATTTCTCTC  >  NC_000913/2338917‑2339406                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cACCAGCAGGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCAt                                                                                                                                                                                                                                                    >  1:77294/1‑248 (MQ=255)   ccAGCAGGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAcc                                                                                                                                                                                                                                               >  2:191158/1‑251 (MQ=255)      gcagGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAccgcc                                                                                                                                                                                                                                            <  2:113508/251‑1 (MQ=255)                     tggttgtctggACGTCCTGAATTTTTCCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑CCATGGGCAATTTTCGCCAGACGGATTTCC‑TTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATgggg                                                                                                                                                                                                                              <  2:239592/240‑1 (MQ=255)                      ggttggCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTGCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGAt                                                                                                                                                                                                                            <  1:327487/251‑1 (MQ=255)                                  gTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCg                                                                                                                                                                                                                 >  1:68351/1‑250 (MQ=255)                                         aTTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTc                                                                                                                                                                                                         >  2:111892/1‑251 (MQ=255)                                            ttttCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGGCCGTCTACCCGCATATAACGCAGCGCGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCAcc                                                                                                                                                                                                      >  1:64541/1‑251 (MQ=255)                                              ttCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGCTTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCa                                                                                                                                                                                                     >  1:56734/1‑250 (MQ=255)                                                                aTCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGAtttt                                                                                                                                                                                   <  1:251069/250‑1 (MQ=255)                                                                   aaCGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGAATTGCCGCCGCAGAGTAGCCGTCGATAGAACCGAAGTTACCATGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTtat                                                                                                                                                                               >  1:43469/1‑251 (MQ=255)                                                                     cGAAATCGACCTTCTCTTTTTCGAGTTCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCATAGTCGCCGTCGATAGAACCGAAGTTACCCTGCCCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAg                                                                                                                                                                             <  1:191158/251‑1 (MQ=255)                                                                      gAAATCGACCGTTTCTTTTTCGTGAGCGGTCA‑CCGTGTCATGGGCAATTTGCGCCAGACG‑ATTTCCCGTATACCGCAGTGCCGCCGCAGAGTCGCCGTCGAGAGAACCGACGTTACCCTGACCGTCTACCAGAATATAACGCAGCGAAAATGGGTGCGCGATGCGGACGATCTGGCCATAGACCGCCGAGTCACCATGTGGATGGTATTTACCGATTACGTCACCAACGACACCGTCAGATTTGTTATAgg                                                                                                                                                                            <  2:177798/251‑1 (MQ=255)                                                                      gAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAgg                                                                                                                                                                            <  1:162584/251‑1 (MQ=255)                                                                            gACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTCCCGATTACGTCACCAACGACACGGGCAGATTTTTTAGTGGCTTTGt                                                                                                                                                                      >  2:288009/1‑251 (MQ=255)                                                                              ccGTCTCGTTGTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTc                                                                                                                                                                    <  2:266304/251‑1 (MQ=255)                                                                                gTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTcc                                                                                                                                                                   >  1:298960/1‑250 (MQ=255)                                                                                gTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTcc                                                                                                                                                                   >  2:154549/1‑250 (MQ=255)                                                                                          cGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCTGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCt                                                                                                                                                        <  1:292227/251‑1 (MQ=255)                                                                                          cGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCTGTCGATAGAACCGAAGTTACCCTGACCGTCCACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCt                                                                                                                                                        <  1:292221/251‑1 (MQ=255)                                                                                                       cAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACg                                                                                                                                                                                             >  2:10237/1‑201 (MQ=255)                                                                                                       cAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACg                                                                                                                                                                                             <  1:10237/201‑1 (MQ=255)                                                                                                          t‑tCATGGGCAATTTTCGCCAGCCGGATTTCC‑TTATAACGCATTGCCGCCGCATAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCTAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGTGCAGATTTTTTATAAGCTTTGTTCCAGTCATTGCCTCGTACGTTCATGGCGt                                                                                                                                        >  2:271745/1‑251 (MQ=255)                                                                                                                gggCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGt                                                                                                                                   <  2:193074/251‑1 (MQ=255)                                                                                                                     attgtCGCCAGACGGATTTCC‑GGATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTAcgcg                                                                                                                              <  2:305100/247‑1 (MQ=255)                                                                                                                          cGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGgt                                                                                                                         <  2:27348/251‑1 (MQ=255)                                                                                                                           gCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGCGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATgg                                                                                                                                           >  2:365776/1‑232 (MQ=255)                                                                                                                           gCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATgg                                                                                                                                           <  1:365776/232‑1 (MQ=255)                                                                                                                                        cc‑GTATACCGCCTTGCCGCCGATGAGGCGCCGTCGATAGAACCGAAGTTACCCTGACCTTCTACCCGCATATAACGGAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTGTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCTTCTGTGTACCCGCTTCAgg                                                                                                           <  2:231715/251‑1 (MQ=255)                                                                                                                                            tataACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTAc                                                                                                                                                   <  1:21691/208‑1 (MQ=255)                                                                                                                                            tataACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTAc                                                                                                                                                   >  2:21691/1‑208 (MQ=255)                                                                                                                                             ataACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGCCACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCAt                                                                                                       >  1:247625/1‑251 (MQ=255)                                                                                                                                               aaCGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGTATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCAtct                                                                                                     >  2:55943/1‑251 (MQ=255)                                                                                                                                                  gCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCgg                                                                                                  >  1:184208/1‑251 (MQ=255)                                                                                                                                                        ccgccgCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCt                                                                                            >  2:225801/1‑251 (MQ=255)                                                                                                                                                           ccgcAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGc                                                                                         <  2:277054/251‑1 (MQ=255)                                                                                                                                                                 agTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCa                                                                                   <  1:39558/251‑1 (MQ=255)                                                                                                                                                                           gATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCcaaca                                                                          >  1:111815/1‑250 (MQ=255)                                                                                                                                                                           gATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCcaaca                                                                          >  1:53863/1‑250 (MQ=255)                                                                                                                                                                            aTAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAAt                                                                        >  2:115106/1‑251 (MQ=255)                                                                                                                                                                              ataaCCGAAGTTACCCTTACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCAATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCCGGCAGCGCACGGCCAACAATGa                                                                      <  1:134241/249‑1 (MQ=255)                                                                                                                                                                              aGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATg                                                                       <  1:68148/250‑1 (MQ=255)                                                                                                                                                                                  ccGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCg                                                                   >  2:192413/1‑250 (MQ=255)                                                                                                                                                                                            ccTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATaa                                                        <  1:333505/251‑1 (MQ=255)                                                                                                                                                                                               ggccGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGGACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAAt                                                       <  2:64541/247‑1 (MQ=255)                                                                                                                                                                                                             atatAACGCGGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGagctctt                                       <  2:29287/251‑1 (MQ=255)                                                                                                                                                                                                                         gagaATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCCACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTc                           >  2:17822/1‑251 (MQ=255)                                                                                                                                                                                                                             aTGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGCGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATg                       >  2:156095/1‑251 (MQ=255)                                                                                                                                                                                                                               ggCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGt                      <  1:193264/250‑1 (MQ=255)                                                                                                                                                                                                                                gCTGCGCCATGCGGGGGTTCGTGCCATAGGCCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGTTg                    <  2:112332/251‑1 (MQ=255)                                                                                                                                                                                                                                            ggACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTCCGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGAGTAACGTCCGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCCATGTTGACCGGTGTAAtt        >  1:24480/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGTTGACCGGTGTAATTtctctc  <  1:273996/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CACCAGCAGGTTAGGAATTTTGGTTGGCATGACGTCCGGAATTTTTTCCGTGCCGTCATAGTTATCAACGAAATCGACCGTCTCTTTTTCGAGATCGGCCATCAGT‑TCATGGGCAATTTTCGCCAGACGGATTTCC‑GTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACCGGCTTCAGGCCATCTCGGACATCTGGCAGCGCACGGCCAACAATGACCGACATCGCATAATCCAGATAGGAGCTCTTCAGCTCTTCCTCAATGTTGACCGGTGTAATTTCTCTC  >  NC_000913/2338917‑2339406

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