Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	24,977	C→T	L863L (CTG→TTG)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	24,977	0	C	T	95.8%	203.8 / ‑0.1	71	L863L (CTG→TTG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (33/35);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

GAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGAGTCAATCGATCTG  >  NC_000913/24731‑25223                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        gAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGtt                                                                                                                                                                                                                                                     <  1:346085/251‑1 (MQ=255) gAACGTGAAAAATACGTCTGCACCGGTTGGTGGTGCGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGtt                                                                                                                                                                                                                                                     <  2:442657/251‑1 (MQ=255)     gTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAAATGACCGCGCTGGGCGATGAATTACGATTTGTCttgtt‑gccc                                                                                                                                                                                                                                                >  2:395204/1‑248 (MQ=255)         aaaaTACGTCTTCACTGGTGAGTGGTTCGAAGGCCTGGTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGGGATGAATTACGATTTGTCTTGTT‑GACCTCCg                                                                                                                                                                                                                                            <  2:421603/251‑1 (MQ=255)         aaaaTACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCg                                                                                                                                                                                                                                            <  2:166254/251‑1 (MQ=255)          caagaCGTCTTCACCGGTGAGTGGTAAGAAGGCCTGTTAGGTATGGCCGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACACAGTCATAGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGATGGAAGCGGCAGTAAACTTGTATACAGAACAGGAACTGGCGGCGAAAATGAACGAGCTGGGCGAGGAATTACGATTTGTCTTGTT‑GACCTACgg                                                                                                                                                                                                                                           >  2:41218/5‑251 (MQ=255)                 tCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTAcc                                                                                                                                                                                                                                    >  1:78723/1‑251 (MQ=255)                 tCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTAcc                                                                                                                                                                                                                                    >  1:128780/1‑251 (MQ=255)                    gcACCGGTGAGTGGTACGAAGGCCTGGTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTAcc                                                                                                                                                                                                                                    <  2:392723/247‑1 (MQ=255)                       ccGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGAGCGGCGAAGTGAACCAAGTCATTGAGCCAGCGCGTGCCGCCAAGAAAGTGGGTGGCTCGATGGAAGCGGACGTAACCTGGTATGCCGAACCGGAACTGGCGGCGGAACATACCGCGCTGGGCGCAGAATTACGATTGGTCTTGTT‑GTCGTCCggggcgacggtcgc                                                                                                                                                                                                                               >  2:267092/1‑238 (MQ=255)                        cGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGTACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAg                                                                                                                                                                                                                             <  1:330562/251‑1 (MQ=255)                            gAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCGGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTTACTTGTTTGAC‑TCCGGCGCTACCggttgaggct                                                                                                                                                                                                                          >  2:40289/1‑241 (MQ=255)                            gAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACta                                                                                                                                                                                                                         <  1:23396/251‑1 (MQ=255)                             aGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGAGGGTGGCTCGCTGGAAGCGGCAGTAACAGTGTATGAAGAACCGGAACTAGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGCTTTGTCTTGTT‑GACCTCAGGCGCTCCCGTTGcggaatat                                                                                                                                                                                                                        >  2:408560/1‑244 (MQ=255)                                 gTACGAAGGCCTGTTTGGGCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGGGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAacg                                                                                                                                                                                                                    <  2:25592/251‑1 (MQ=255)                                      aaGGCCCGTTTGGTCTGGCAGACAGGGAAGCTATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCa                                                                                                                                                                                                               <  2:109675/251‑1 (MQ=255)                                       aGGCATGTTTGGTCTGGCAGACAGTGAAGCAATGAATGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAATAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCAc                                                                                                                                                                                                              <  2:54571/251‑1 (MQ=255)                                        ggCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCAcc                                                                                                                                                                                                             <  2:178712/251‑1 (MQ=255)                                         gCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCGGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACct                                                                                                                                                                                                            <  2:478911/251‑1 (MQ=255)                                                   tcTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGCCGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGAGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGCCTCTAACGACGCACCTGCTGATGCTc                                                                                                                                                                                                  >  1:42992/2‑251 (MQ=255)                                                      ggCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTTCCGTTGCAGTCTATAACGACGCACCTGCTGATGCTca                                                                                                                                                                                                 <  1:199089/249‑1 (MQ=255)                                                      ggCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagc                                                                                                                                                                                               <  2:383191/251‑1 (MQ=255)                                                               tGAAGCAATGTACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGACGTGAACAAAGTCATAGAGCCAGCGCGGGCCGACAAGAAAGTGGGTGGCTAGATGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGCTTTGTCTTGTT‑GACCTCAGGAGCTACCGTGGCCGACAATAACGACGCACCTGCTGATGCTAAGCAGAGCgc                                                                                                                                                                                        >  2:326395/1‑248 (MQ=255)                                                                 aaGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGTTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGt                                                                                                                                                                                     >  1:100406/1‑250 (MQ=255)                                                                    caaTGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTc                                                                                                                                                                                 <  2:42992/249‑1 (MQ=255)                                                                      aTGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGCTGCTCAGCAGAGCGCAGTCCTCa                                                                                                                                                                                >  2:246051/1‑250 (MQ=255)                                                                              tGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTg                                                                                                                                                                        >  1:203752/1‑250 (MQ=255)                                                                                     tGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTcgc                                                                                                                                                                >  1:102297/1‑251 (MQ=255)                                                                                      gggACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTcg                                                                                                                                                                 <  1:106631/249‑1 (MQ=255)                                                                                            tgCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGGGGCTCGCTGGAAGCGGCAGTAACCTGGTATGCAGAAACGGACCTGGCGGCGGTACTGGCCGCGCGTGGCGATGAATTACGATTTGTATTTTT‑GACCTCCGGCGCTACCGTTGCAGGCTATAACGACGCAACTGCTGATGCTCAGCCGAGCGAAGTACTCCACGGGCTGCAAGTAGCGTTGAGt                                                                                                                                                         >  2:134061/2‑251 (MQ=255)                                                                                              cggtggAAAGGGCGTGGCGAAGTGAACCAAGTCATGGAGCAAGCGCGTGACGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTTGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTTGCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCTGCAGAGCGAGGTACTCAAAGGGCTGAAAGTCGCGTTGAGTaa                                                                                                                                                       <  2:336199/246‑1 (MQ=255)                                                                                                gTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCCTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTaaa                                                                                                                                                      <  1:408560/250‑1 (MQ=255)                                                                                                        tGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGCTGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTTGGTAAAGCCGAGGGt                                                                                                                                             >  2:13077/1‑251 (MQ=255)                                                                                                                         aaaGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgctg                                                                                                                            <  1:502767/251‑1 (MQ=255)                                                                                                                         aaaGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGAGTTGAGTAAAGCCGGAGGTGAGAAGGGCCAACgctg                                                                                                                            >  2:113195/1‑251 (MQ=255)                                                                                                                           aGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgctgct                                                                                                                          >  1:444774/1‑251 (MQ=255)                                                                                                                                     cAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCCAAAGGCTGAAAGTCGCGTTTCGTAAAACCGAAGGTGAGAAGGGCCAACGCTGCTGGCCCTacac                                                                                                                >  2:189876/1‑251 (MQ=255)                                                                                                                                      aaGCGCGGGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAAGTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACAcc                                                                                                               <  2:83138/251‑1 (MQ=255)                                                                                                                                        gcgcGTGCAGACAAGAAAGGGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCa                                                                                                             <  2:409592/251‑1 (MQ=255)                                                                                                                                            gTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGAt                                                                                                         <  1:106345/251‑1 (MQ=255)                                                                                                                                             tGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCg                                                                                                                                                                                                                                            >  1:68793/1‑119 (MQ=255)                                                                                                                                             tGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCg                                                                                                                                                                                                                                            <  2:68793/119‑1 (MQ=255)                                                                                                                                               ccGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGGGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTc                                                                                                      >  2:418357/1‑251 (MQ=255)                                                                                                                                                                   tCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGTTTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCACCAGAGCGAAGTAATCAAAGGGCTGAAAGTCGGGTTGAGTAAAGCCGAAGGTGAGAGATGCCCACGCTGCTGGCACCACCCCCAGGATGTCGGCAAGGTGGCGGCACACGc                                                                                  >  2:378299/1‑251 (MQ=255)                                                                                                                                                                         gAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAAt                                                                            >  1:431894/1‑251 (MQ=255)                                                                                                                                                                                cAGTAACCTTGTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGc                                                                     <  2:100406/251‑1 (MQ=255)                                                                                                                                                                                     aCCTTGTATGCCGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGCTTTGTCCTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATCACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCCCTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGACGCtg                                                                >  2:456248/1‑251 (MQ=255)                                                                                                                                                                                     aCCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCtg                                                                <  1:456248/251‑1 (MQ=255)                                                                                                                                                                                          gtattaaGAACGGGAACTGTCGGCGAGACTGCCCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCGCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCGGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCa                                                           <  2:233292/245‑1 (MQ=255)                                                                                                                                                                                          gTATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGGGGCCGCttggtc                                                            >  2:495078/1‑245 (MQ=255)                                                                                                                                                                                           tATGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAGGTGCCCACGCTGCTGGAACTACACCCAGGATGTCGGCCAGGTGGCGGAAAACGCAGAAATCTGCGGCCGCTGAGTc                                                            >  2:182197/1‑249 (MQ=255)                                                                                                                                                                                            aTGCAGAACCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAg                                                          <  2:298451/250‑1 (MQ=255)                                                                                                                                                                                              gcagcaCCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTAAGATTTGTCCTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTACAGCCGCAGGTGAGCAGTGCCCACCATGCTGGCACTACACCCAGGATGTCGGCAAGGTGTCGGAACACGCATAAATATGCGGCCGcctggtacgca                                                        >  2:94457/1‑240 (MQ=255)                                                                                                                                                                                                   aCCGGCACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGGTGATAGGCGCGTTGTGTAAAGCCGAAGGTGCGAAGTGCCCACGCTGCTGGCACTCCACCCAGGAAGTCGGCGAGGGGGCGGAACACGCCGAAATCTGCTGCGGCTGGGTCAGAAACggcg                                                  >  2:115823/1‑248 (MQ=255)                                                                                                                                                                                                   aCCGGAACTGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCGCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGAGTTGAGTAAAGCCGAAGGTGAGAAGTGCTCACGCTGCTGGCACTACACCCAGGAGGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGTCCGCTGTGTCAGCAACGTCg                                                  <  1:197613/251‑1 (MQ=255)                                                                                                                                                                                                      ggtaCTGGCGGCGAAACTGACCGCGCTGGGCGATGCATTACGATTTGTCTTGTT‑GTCCGCCGGCGCTACCGTTGAAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGATCTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCAGATGGTGAGAAGTGCCCACGCTGCTGGCGATAACCCCAGGATGTCGGCAGGGGGGCGGAAACCCCCGAAATATGAGGCCGCTGTGTGGGCAACGTCGCcg                                               >  2:296191/4‑251 (MQ=255)                                                                                                                                                                                                                  gAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGaa                                    >  1:446701/1‑250 (MQ=255)                                                                                                                                                                                                                     aCTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAAc                                <  2:128780/251‑1 (MQ=255)                                                                                                                                                                                                                      cTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACg                               >  1:259007/1‑251 (MQ=255)                                                                                                                                                                                                                         aCCGCGCTGGGCGCTGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCt                                                                                                                                                                         >  1:157834/1‑110 (MQ=255)                                                                                                                                                                                                                         aCCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGCAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTCCACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGAGTCAGCAACGTCGCcggtg                                            >  2:433568/1‑235 (MQ=255)                                                                                                                                                                                                                         aCCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCt                                                                                                                                                                         <  2:157834/110‑1 (MQ=255)                                                                                                                                                                                                                         aCCGCGCTGGGCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCcggtg                                            <  1:433568/235‑1 (MQ=255)                                                                                                                                                                                                                          ccGCGCTGGGCGAGGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGACGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAg                           <  2:178174/251‑1 (MQ=255)                                                                                                                                                                                                                                 gggCGATGAATTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGc                      <  1:169625/249‑1 (MQ=255)                                                                                                                                                                                                                                       tGAATTAAGATTTGTCTTGTT‑GACCTCCAGCGCTACCGTTGCAGACTCTAACGACGCACATGCTGCTGCTCGGCAAAGCGCAGTCCTCAAACAGCAGCCAGACGCGATGCGTAAAGCCGAAGGCGAGCAGTGCCCACGCTGCTGGCACTACAACCTGGCTGTCGGCAAGGTGACGGATCACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGCCGCCGGTGACGGTGGCAACCGTAAGTTGGCCtcatgcg              >  2:255153/1‑249 (MQ=255)                                                                                                                                                                                                                                         aaTTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCTTTGCAGACTATAACGGCGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCCCTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGAGTc            <  2:113601/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaTTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGAGTc            >  2:411695/1‑251 (MQ=255)                                                                                                                                                                                                                                          aTTACGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGAGTcc           >  2:169452/1‑250 (MQ=255)                                                                                                                                                                                                                                             aCGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGAGTCAAt         <  1:379165/250‑1 (MQ=255)                                                                                                                                                                                                                                              cGATTTGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGAGTCAATCg       <  1:84183/251‑1 (MQ=255)                                                                                                                                                                                                                                                   tGTCTTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCCGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGAGTCAATCGATCTg  >  2:374258/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        GAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTT‑GACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCGGTGACGGTGAAAAACGTAAGTTTGCCTGATGAGTCAATCGATCTG  >  NC_000913/24731‑25223

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