Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	696,470	C→T	noncoding (35/75 nt)	glnX ←	tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	696,470	0	C	T	93.8%	43.9 / ‑3.1	17	noncoding (35/75 nt)	glnX	tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (9/6);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGAT  >  NC_000913/696229‑696681                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                     gCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATCAGAATCCg                                                                                                                                                                                                              <  1:101429/250‑1 (MQ=255)        tgtTGAAACGTCAGGCGCTGAAATCTGCGTTAGGGTTGTAGCCTTACAACATTGCCGGAAGCGTGCCAACCGTAGGACGGATAAGACGCGGCAGCGTCGCATACGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTAAGAGGATTCGACACCT‑GGGATACAGGAATTAGAATCCGGTGCCtt                                                                                                                                                                                                       >  2:320867/1‑250 (MQ=255)           tGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCttcacg                                                                                                                                                                                                   >  2:107940/1‑247 (MQ=255)               aCGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGCTGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTg                                                                                                                                                                                               >  1:152568/1‑251 (MQ=255)                 gTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGGCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGc                                                                                                                                                                                             <  2:178627/251‑1 (MQ=255)                                          gTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCCTCCGACATTGAAGGATAAGACGTGTCAACATCGCCTTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTAc                                                                                                                                                                     >  1:290494/1‑250 (MQ=255)                                                                          cAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGCAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGACTTGAATAAATTGACTGGGGt                                                                                                                                     >  1:12812/1‑250 (MQ=255)                                                                                                          ggaTCCGACATTGAAGGATAAGACGTGTCAACAGCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAAGTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACCGCGGAATGCCGGAATCAg                                                                                                   <  2:24881/249‑1 (MQ=255)                                                                                                                              aGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATTGGTTGAACTTACAAGGTAAGGGCCTTGCTTAAATTGCCTGGG‑TATATAGGATACGAACCTCGGTATGCGGGAAAAAGACTCCGGTGCCTTACCGCtt                                                                               >  1:162169/1‑251 (MQ=25)                                                                                                                                 cGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTTCCTTACCGCTTGGc                                                                            <  1:346003/251‑1 (MQ=255)                                                                                                                                                  cGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTg                                                           <  1:247723/251‑1 (MQ=255)                                                                                                                                                                       aaGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACata                                      <  1:164114/251‑1 (MQ=255)                                                                                                                                                                          cAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGGACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGCTGCACTTGCAAGGTACGCGTCTTGAATAAATTGACTGGGGTA‑CGAGCATTCGAAGCTCGGCATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGAGATACACCAAGAAAATGGTTTTGGATTTGcagcacgtattt                                   >  1:162213/1‑240 (MQ=25)                                                                                                                                                                               aGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACATCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATAc                              >  1:442153/1‑251 (MQ=255)                                                                                                                                                                                    gcggAAGCGAGATTTTTCAGTGTAGCTGGGGTACGAGGACTCGA‑AACATTGAAGGACGAAATTAGTATCCGGTGCCTTCCCGATTGAGGGTACGCCAACTGGGTGCACTTACAAGGTAAGCGGCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCa                         <  2:401892/248‑1 (MQ=21)                                                                                                                                                                                              gCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTtggtg               >  1:277470/1‑251 (MQ=255)                                                                                                                                                                                                           ggccggggTACGAGGATTCGA‑AACTCGGAATGCCGGAATGAGAAACCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGAt  <  2:6221/247‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                     GCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGA‑ACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTA‑CGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGAT  >  NC_000913/696229‑696681

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