Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,307,372	A→C	M143R (ATG→AGG)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,372	0	A	C	95.7%	110.1 / ‑0.5	46	M143R (ATG→AGG)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (27/17);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

GACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCA  >  NC_000913/4307143‑4307600                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     gACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTcccc                                                                                                                                                                                                                                 >  1:589344/1‑235 (MQ=255)     aaaGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCATGCCACACCGTGTGCACCCTTCCACGGCGCACCGTTCATCCCAACTGCGAACCCCCTCCCCGTAcc                                                                                                                                                                                                                            >  2:455739/1‑236 (MQ=255)         gggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCgag                                                                                                                                                                                                                         >  1:175745/1‑235 (MQ=255)            cGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAg                                                                                                                                                                                                                     >  1:561268/1‑236 (MQ=255)                   cGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCCCACTGCGAACCCCCTCCCCGTACCGAGATCAgccgcc                                                                                                                                                                                                               >  2:338280/1‑235 (MQ=255)                               aaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                                                                                                                                                                                  >  2:39992/1‑236 (MQ=255)                                    cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCGGAACCAGACGTTGctg                                                                                                                                                                                              >  2:603269/1‑233 (MQ=255)                                    cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCCGAACCAGTTGTTGCGTa                                                                                                                                                                                             >  1:430856/1‑236 (MQ=255)                                    cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCCGTTGTTGCGt                                                                                                                                                                                              >  1:206768/1‑235 (MQ=255)                                                                    tACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGTCGGTTTTCTACCACGTCCCAGGAGAGCTg                                                                                                                                                             >  2:12224/1‑236 (MQ=255)                                                                     aCCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGTCGGTTTTCTACTACGTCCCAGGagag                                                                                                                                                                >  1:325779/1‑232 (MQ=255)                                                                                    ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGCTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCgg                                                                                                                                             >  1:349397/1‑236 (MQ=255)                                                                                                  ttGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCCGTGTTTGCGTAAGGCGGCTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCgg                                                                                                                               >  2:689476/1‑236 (MQ=255)                                                                                                  ttGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTTGCGGGAAAACTCAACCg                                                                                                                                >  2:65455/1‑235 (MQ=255)                                                                                                                          ccACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCtt                                                                                                       >  1:226739/1‑236 (MQ=255)                                                                                                                           cacaCGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCCGTTGAGTTGCAGGTTAACGTCGCGGGCAAACTCAACCGGCCAAATCAGCGTAGTTTCGAGCATa                                                                                                      >  2:580199/1‑236 (MQ=255)                                                                                                                                           ggTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTAATCGGCggaatccatt                                                                                       >  2:353965/1‑226 (MQ=255)                                                                                                                                                tatCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTGGCGGGAACACTCAACCGGACAATTCATCGGATTTTCGAGCTTTTCGGCGCGAGCTAATAAATc                                                                                  >  2:184220/1‑235 (MQ=255)                                                                                                                                                tatCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGCCAATTCAACGTATTTTCCGGCTTTACGGCGAGATccacttaatc                                                                                  >  2:640757/1‑226 (MQ=255)                                                                                                                                                tatCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATcc                                                                                 >  1:17138/1‑236 (MQ=255)                                                                                                                                                tatCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATcc                                                                                 >  1:242974/1‑236 (MQ=255)                                                                                                                                                          ggggATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGCCAATTCAGCGGATTTTCGAGCTTATCGGGCGGATCATATAAATcagcccccggta                                                                       >  2:25405/1‑229 (MQ=255)                                                                                                                                                          ggggAAATGACCCAGGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGtt                                                                        <  1:304212/236‑1 (MQ=255)                                                                                                                                                                atGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCa                                                                 >  1:441994/1‑236 (MQ=255)                                                                                                                                                                      cA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAg                                                           <  1:455739/236‑1 (MQ=255)                                                                                                                                                                         gTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTACGCGTATTTTCGAGCTTATCGGCGAGATCAAATAAATCCGACGCTGTTAACGGCAGGGTTGGc                                                         >  2:113190/1‑236 (MQ=255)                                                                                                                                                                           tCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGt                                                       <  1:640757/236‑1 (MQ=255)                                                                                                                                                                           ccGCCTGCCATACCGGGGGCACCCGGTCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGt                                                       <  1:664580/235‑1 (MQ=255)                                                                                                                                                                                          tgtgCCCCCGTCCGCGGCGCAGCGTTAATCCACACTGCGAACCCCCTCCACGTACCGAGATAGGCCGGCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGGAGGGTGCAGGTTAACGGCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCg                                        <  2:126274/236‑1 (MQ=255)                                                                                                                                                                                          tgtgCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGGATTTTCGAGCTTATCGGCGAGAGCACATAAATCCGCCGCGGTTAACGGCAGGGTTGGGGGAGGAAAAatgatgtg                                        >  2:291739/1‑234 (MQ=255)                                                                                                                                                                                                cccGTCCGCGGCCGACCGTTCATCCACAAGTCGAACCCCCTCCCCGTCCCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTCCGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTAt                                   <  2:206768/235‑1 (MQ=255)                                                                                                                                                                                                  cGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGc                                                                                                 >  1:352959/1‑171 (MQ=255)                                                                                                                                                                                                  cGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGc                                                                                                 <  2:352959/171‑1 (MQ=255)                                                                                                                                                                                                  cGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCt                                 <  1:12224/235‑1 (MQ=255)                                                                                                                                                                                                    tCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAgttg                                                                                                                                                                                                    >  1:599599/1‑70 (MQ=255)                                                                                                                                                                                                    tCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAgttg                                                                                                                                                                                                    <  2:599599/70‑1 (MQ=255)                                                                                                                                                                                                         ggCGCACGGTTCATCACCACTGCGACCCCCCGCCCCGTACCGGGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGAc                         <  2:242974/236‑1 (MQ=255)                                                                                                                                                                                                                   tCATCCACCAGGCAAACCCCCTCCCGTTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCTGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGa               <  2:441994/236‑1 (MQ=255)                                                                                                                                                                                                                        cacaCTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATcc          <  1:25405/236‑1 (MQ=255)                                                                                                                                                                                                                        cacaCGGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAGCCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATcc          <  2:589344/236‑1 (MQ=255)                                                                                                                                                                                                                         acccTGCGAACCCCCTCCCCGTACCGAGATAGGCCACCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATccc         <  1:338280/233‑1 (MQ=255)                                                                                                                                                                                                                         acacTGCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATccc         <  1:587843/236‑1 (MQ=255)                                                                                                                                                                                                                              gCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTGCGTGCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCAc    <  2:430856/236‑1 (MQ=255)                                                                                                                                                                                                                              gCGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCAc    <  2:17138/236‑1 (MQ=255)                                                                                                                                                                                                                               cGAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGGGTAGAAATAATGGTGAGGTTATCTTTACTGACAAGCGCCGCGAatcacatcac    >  2:368174/1‑235 (MQ=255)                                                                                                                                                                                                                                gAACCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCa  <  1:39992/236‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     GACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCA‑GTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCA  >  NC_000913/4307143‑4307600

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