Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,624	C→T	noncoding (71/76 nt)	lysT →	tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,624	0	C	T	91.7%	29.6 / ‑5.4	12	noncoding (71/76 nt)	lysT	tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (6/5);  minor base A (0/1);  total (6/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGG  >  NC_000913/780433‑780817                                                                                                                                                                                                |                                                                                                                                                                                                  gtgCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAg                                                                                                                                            <  1:16372/247‑1 (MQ=255)                               aTCGAAATTTTTGTTGCGCTCAAGTCTGAAATCATTAATATATGCCTCCGTTGCCACGAGATAGCAAACAAACCGAAAGCAACGAAAAAGGGGGTCGTTAGCTCAGTTGGTAGAGTAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATTCTGCACGACTCTCCAAACGCTAAGGTGGAAGCGGTAGTAAAGCGGGAAGTATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCAt                                                                                                          <  1:135359/251‑1 (MQ=37)                                 cGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCTCAGGTTCGAATCCTGCACGACACACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGTATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATcc                                                                                                        <  1:143592/251‑1 (MQ=255)                                  naaaTTTTTGTTGCGCTCAAGTCTGAACTCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAACACGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGCATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGCGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCt                                                                                                       >  1:31086/2‑251 (MQ=255)                                              gcgcTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGCAAGCAACGAAAACGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGACGGATAACGTTGCATGAGCAACGGCCCGCAGGGCGAGACGAAGTCGAGTCATCCTGCACGACccacc                                                                                           >  1:82707/1‑251 (MQ=255)                                               cgcTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACccacc                                                                                           >  1:10520/1‑250 (MQ=255)                                                        ntGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACAt                                                                                   >  1:212383/2‑249 (MQ=255)                                                                    atatatGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCNc                                                                                                   <  1:277719/221‑1 (MQ=255)                                                                                                   caaacCGAAAGCAACGCAAAAGTGGGTCGTTAGCTACGTCGGTAGAGCAGTTGACTTTTAATCTCTTGGCCGACGGATCGAATCATGCACGGCTCCCCAATCGCTAAGGTGGGAGCGGTAGTACAACGTGAAGGATAACGTTGCATGAGCAACGGCCAGAAGGGCGAGACGAAGTCGAGTCATCCTGCCCGACCCACCACTAAAATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCt                                      >  1:226103/1‑251 (MQ=25)                                                                                                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcggcgn  <  1:284197/232‑2 (MQ=255)                                                                                                                                                               aTCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTATCTCAGCTGGGAGAGCAc                            <  1:96068/201‑1 (MQ=255)                                                                                                                                                                             ggTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGgaga                                >  1:117815/1‑183 (MQ=255)                                                                                                                                                                                                |                                                                                                                                                                                                  GTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGG  >  NC_000913/780433‑780817

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