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%	37.1 / NA	16	noncoding (71/76 nt)	lysT	tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/7);  minor base G (0/1);  total (8/8)
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.28e-01

ATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAG  >  NC_000913/780389‑780857                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          aTGCAAAAAACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATACGCAGTTGAGTGCTCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGCGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAAtagcaacgg                                                                                                                                                                                                                            >  2:136266/1‑243 (MQ=38)     aTAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTacggagg                                                                                                                                                                                                                         >  2:99327/1‑244 (MQ=255)               aCCAGAAGTCGCATTATTTCTGGTCTTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGTTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACGAACCAATCGCTAAGGATTAAGCGgtagta                                                                                                                                                                                                              <  2:16813/251‑1 (MQ=38)                           attattTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGa                                                                                                                                                                                                  <  1:26236/251‑1 (MQ=255)                                  cTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGtt                                                                                                                                                                                           >  1:235539/1‑251 (MQ=255)                                    ggTCGTGGCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTg                                                                                                                                                                                          <  2:73659/250‑1 (MQ=255)                                                                                                                                 cACGGGATATCAAACAACCAGAAAGCAACGAAAAAGTGGGTCGATAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGACTCCTGCACGACTCACCAATCGCTAAGGAGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGACGAGTCATCCTGCACGACCCACCACTAACATAGTTCGTTGTAGTATCCAGCGTAGAATCggg                                                                                            >  1:199544/1‑251 (MQ=25)                                                                                                                                   cGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGATGACTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGACTCACCAATCGCTAAGGAGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTg                                                                                          <  2:199544/251‑1 (MQ=255)                                                                                                                                         atCAAACAAACCGAAAGCAACGAAAAAGGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACTAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGTCCCGAATGGCGATACGAAGTCGAGTCATCCTGCACGACCCGCCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGc                                                                                    <  1:108435/251‑1 (MQ=37)                                                                                                                                                                aaaaGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGAAAGCGGTAGTAAAACGTGACGGATAACGTTGCATGAGCAACGGCCCGCAGGGCGAGACGAAGTCGAGCCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCtt                                                             >  1:166699/1‑251 (MQ=255)                                                                                                                                                                              aGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCAGCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTc                                               <  1:175998/251‑1 (MQ=255)                                                                                                                                                                                 ccAGTTGGTAGAGCAGTTGACTTGTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcggc                                            <  2:235539/250‑1 (MQ=255)                                                                                                                                                                                   aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATAGCTAAGGAGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCACGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAGCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCAACTCCCTTACAAGGAGGGGGTcggcgg                                          >  1:90306/1‑251 (MQ=255)                                                                                                                                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGACGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGAAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGCGGGGGTCGGCGGTTCGATCCCGGCCTCACCCCCCACCgg               >  2:157639/1‑251 (MQ=38)                                                                                                                                                                                                                       cAGGTTCGAATACTGCACGACTCACCAATCGCTACGGTGGAAGCGGAAGTAAAACGTGAAGGATAACGTTGCATGAGCAAAGGCCCGAAGGGCGAGACGAGGTCGAGTCCTCCTGCACGACCCAACACTAACATAGGTGGTTGTAGTGTCCAGTGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGATCGGGCCCGTCGTcacccaccccccggtagtcagg      >  2:106134/1‑237 (MQ=25)                                                                                                                                                                                                                            tCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAg  <  1:99997/250‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          ATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAG  >  NC_000913/780389‑780857

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