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	93.8%	39.1 / ‑6.0	16	noncoding (71/76 nt)	lysT	tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/5);  minor base G (1/0);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCT  >  NC_000913/780379‑780854                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       aCGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTTGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACTACtcacca                                                                                                                                                                                                                                   <  2:205501/251‑1 (MQ=255)                tAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGa                                                                                                                                                                                                                     >  1:18809/1‑250 (MQ=255)                             gAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAAc                                                                                                                                                                                                       <  1:113381/251‑1 (MQ=255)                                  cGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTACGGTGGAAGCGGTAGTAAAACGTGa                                                                                                                                                                                                   >  1:239893/1‑250 (MQ=255)                                                                         tGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGAAGTAAAACGTGAAGGATAACGTTGGATGAGCAACGGCCCGAAGGGCGAGACg                                                                                                                                                           >  1:252573/1‑251 (MQ=255)                                                                             tGATCTACCTCGGAGTTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGGAACGGAAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGt                                                                                                                                                       >  1:41413/1‑251 (MQ=37)                                                                                              tttGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTACGGTGGAATCAGTAGTTAAACGTGAAGGATAACGATGCATGAGCAACGGCCCGTATGGCGAGACAAAGTCGAGAGATCCTGCACg                                                                                                                                       >  1:192971/1‑250 (MQ=37)                                                                                                         cAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGCAGGGCGAGACGAAGTCGAGTAATCCTGCACGCCccacccata                                                                                                                            >  1:124081/1‑246 (MQ=255)                                                                                                                              atGCCTCCGTTGCCACGGGATATCATACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAACCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGTTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGATGAAGTCGAGTCATCCTGCACGACCCACCACTAACAGAGTTAGTTGTAGTATCCa                                                                                                      <  2:213105/251‑1 (MQ=25)                                                                                                                                           cACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGAAGGTACGAATCTCGCACGACGCAACAGTCGCAAAGGCGGAGGCGGACGTAAAACGAGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGCCGAAGAAGAGGCATTCTGCACGACCCACCACAAAAATATGTGGTTGCAGGATCCAGCgtagaataggg                                                                                         >  2:241059/1‑244 (MQ=19)                                                                                                                                                 atatCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTa                                                                                   <  2:18809/251‑1 (MQ=255)                                                                                                                                                                           aaaGTGGGTCGTTCGCTCAGTTGGTAGAGCAGTTGACTTTTAATCCATTGGTCGCAGGTTAGAATCCTGCACGACTCACCACTATCAACGGTGGACGCGGTAGTCAACCGTGAAGGATAACGTTGCATAAGCAACGGCCCGCAGGGCGAGCCGAAGGAGAGTCATCCTGCACCACCCACCACAAACATCGTTAGTTGCAGTATCCAGCGTAGTAGCGGGTGATTAACTCAGCTTGGAGCGCACCTACCtt                                                          >  2:324440/1‑250 (MQ=16)                                                                                                                                                                                    cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAggg                                                >  1:106350/1‑251 (MQ=255)                                                                                                                                                                                                               cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCCCTAACATAGGTAGTTGGAGGATCCAGCGTAGAATCGGCAGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAGGGAGGGGCTaggcggggtgaaccccgggttacc                     >  2:37073/1‑233 (MQ=25)                                                                                                                                                                                                                  ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcaccca                   <  2:106350/250‑1 (MQ=255)                                                                                                                                                                                                                            tgTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGTCGAAGTCGAGTCATCCTGCACGACCCACGACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGGTTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTc         <  1:194217/249‑1 (MQ=255)                                                                                                                                                                                                                                  aGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCt  >  1:294656/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       ACGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCT  >  NC_000913/780379‑780854

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