Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,567,413	A→C	intergenic (‑122785/‑74351)	unknown ← / → unknown	tRNA‑Met‑CAT/tRNA‑Gly‑GCC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,567,413	0	A	C	88.9%	13.9 / ‑3.8	9	intergenic (‑122785/‑74351)	unknown/unknown	tRNA‑Met‑CAT/tRNA‑Gly‑GCC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (5/3);  minor base T (0/1);  total (5/4)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ACCTAACCACCACATGGCAAAATTGCCCTCGGGAACAACTTCGTTCTCAATTTCTTCATTTAACCAAGTGCCCCATTCAGGAAAAGTAGAAAGGATCCAACTTTCACGGGTGATTTCATTAACTTTAGCCATACAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAAT‑ATTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAATTCCAGACTTAGCTAACATTTTAACCTTAGATAACAACGCCAACTCTTGATTCGAAATTTCTAATCACACATAAGTATTTCAATAAAATAAATCGTGACATTGCTCACATTTTATCAAAAAACATCTTTGAAAAATCATTATAAATATTTATTATCCGCCGCCAACAAACGAAATTTACGCATGAATAT  >  1/1567203‑1567648                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            aCCTAACCACCACATGGCAAAATTGCCCTCGGGAACAACTTCGTTCTCAATTTCTTCATTTAACCAAGTGCCCCATTCAGGAAAAGTAGAAAGGATCCAACTTTCACGGGTGATTTCATTAACTTTAGCCATACAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAAT‑CTTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAAttt                                                                                                                                                                                                     >  1:713694/1‑251 (MQ=255)          ccacATGGCAAAATTGCCCTCGGGAACAACTTCGTTCTCAATTTCTTCATTTAACCAAGTGCCCCATTCAGGAAAAGTAGAAAGGATCCAACTTTCACGGGTGATTTCATTAACTTTAGCCATACAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAAT‑CTTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAAtt                                                                                                                                                                                             >  1:414676/1‑250 (MQ=255)                      aTTGCCCTCGGGAACAACTTCGTTCTCAATTTCTTCATTTAACCAAGTGCCCCATTCAGGAAAAGTAGAAAGGATCCAACTTTCACGGGTGATTTCATTAACTTTAGCCATACAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAAT‑CTTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAATTCCAGACTTAGCTa                                                                                                                                                                                <  1:189807/251‑1 (MQ=255)                                             tctcAATTTCTTCATTTAACCAAGTGCCCCATTCAGGAAAAGTAGAAAGGATCCAACTTTCACGGGTGATTTCATTAACTTTAGCCATACAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAAT‑CTTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAATTCCAGACTTAGCTAACATTTTAACCTTAGATaacaac                                                                                                                                                         <  1:394707/251‑1 (MQ=255)                                                          aTTTAACCAAGTGCCCCATTCAGGAAAAGTAGAAAGGATCCAACTTTCACGGGTGATTTCATTAACTTTAGCCATACAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAAT‑CTTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAATTCCAGACTTAGCTAACATTTTAACCTTAGATAACAACTCCAACTCTTGAt                                                                                                                                            >  1:39274/1‑251 (MQ=255)                                                                                      gTAGAAAGGATCCAACTTTCACGGGTGATTTCATTAACTTTAGCCATACAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAAT‑CTTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAATTCCAGACTTAGCTAACATTTTAACCTTAGATAACAACTCCAACTCTTGATTCGAAATTTCTAATCACACATAAGTAtt                                                                                                                >  1:557669/1‑251 (MQ=255)                                                                                                                                    tACAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAATCTTTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAATTCCAGACTTAGCTAACATTTTAACCTTAGATAACAACTCCAACTCTTGATTCGAAATTTCTAATCACACATAAGTATTTCaa                                                                                                            <  1:480584/210‑1 (MQ=255)                                                                                                                                      cAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAAT‑CTTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAATTCCAGACTTAGCTAACATTTTAACCTTAGATAACAACTCCAACTCTTGATTCGAAATTTCTAATCACACATAAGTATTTCAATAAAATAAATCGTGACATTGCTCACATTTTATCAAAAAACAt                                                                  >  1:413019/1‑249 (MQ=255)                                                                                                                                                                                                    tttCATTAATGGAAT‑CTTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAATTCCAGACTTAGCTAACATTTTAACCTTAGATAACAACTCCAACTCTTGATTCGAAATTTCTAATCACACATAAGTATTTCAATAAAATAAATCGTGACATTGCTCACATTTTATCAAAAAACATCTTTGAAAAATCATTATAAATATTTATTATCCGCCGCCAACAAACGAAATTTACGCATGAatat  <  1:795975/251‑1 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            ACCTAACCACCACATGGCAAAATTGCCCTCGGGAACAACTTCGTTCTCAATTTCTTCATTTAACCAAGTGCCCCATTCAGGAAAAGTAGAAAGGATCCAACTTTCACGGGTGATTTCATTAACTTTAGCCATACAATAACTCCTTTTTTTAGGTTTAGTAGAATAATTTAATCATTTATAATCATAAACAATCATTTTCATTAATGGAAT‑ATTTTTTTTTAAATTTTTCTATTAAAATATGTTAAAAATTTATTAAATTCCAGACTTAGCTAACATTTTAACCTTAGATAACAACGCCAACTCTTGATTCGAAATTTCTAATCACACATAAGTATTTCAATAAAATAAATCGTGACATTGCTCACATTTTATCAAAAAACATCTTTGAAAAATCATTATAAATATTTATTATCCGCCGCCAACAAACGAAATTTACGCATGAATAT  >  1/1567203‑1567648

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