Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	855,535	A→G	intergenic (‑5878/‑66683)	unknown ← / → unknown	tRNA‑Met‑CAT/tRNA‑Leu‑GAG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	855,535	0	A	G	94.1%	37.9 / ‑4.3	17	intergenic (‑5878/‑66683)	unknown/unknown	tRNA‑Met‑CAT/tRNA‑Leu‑GAG
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (4/12);  total (4/13)
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.71e-01

GTTAGTAGTCCGTGATTAATAAACTCACCAAGAATTTTACCCTCTTTATCACGATATTTACTTGCCTTATAGGAAAATATATCCTGTAAAACGATTTGACCGTTTTCCATTCCCATAACCTCAGTGATATTCATAATTTTTCTTGAACCATCATTTAAACGGGACGCCTGCACTATTATGTTTACTGCTGAAGAGATATTACGGCGAATTGCCTCTAACGGTAGGGATGCATTCGACATCATAACCATACTTTCCAAGCGAGATGTTGCATCGCGAGGGCTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCGGTATTCATTGCTTGTAACATTTGGAATGCTTCTCCACCACGGCATTCGCCCACAATA  >  1/855401‑855775                                                                                                                                       |                                                                                                                                                                                                                                                 gTTAGTAGTCCGTGATTAATAAACTCGCCAAGAATTTTACCCTCTTTATCACGATATTTACTTGCCTTATAGGAAAATATATCTTGCAAAACGATTTGACCGTTTTCCATTCCCATAACTTCAGTGATATTCATAATTTTTCTTGAACCATCGTTTAAGCGGGACGCCTGCACTATTATGTTTACTGCTGAAGAAATATTACGGCGAATTGCCTCCAACGGTAGGGATGCATTCGACATCATAAcc                                                                                                                                   <  1:507049/246‑1 (MQ=255)                                                          ttACTTGCCTTATAAGAAAAAATATCTTGTAAAACAATTTGCCCATTTTCCATACCCATTATTTCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAAcc                                                                     <  1:323313/251‑1 (MQ=255)                                                          ttACTTGCCTTATAAGAAAAAATATCTTGTAAAACAATTTGCCCATTTTCCATACCCATTATTTCAGTAATATTCATGATTCTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAAcc                                                                     <  1:438528/251‑1 (MQ=255)                                                              ttGCCTTATAAGAAAAAATATCTTGTAAAACAATTTGCCCATTTTCCATACCCATTATTTCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCg                                                                 <  1:833597/251‑1 (MQ=255)                                                              ttGCCTTATAAGAAAAAATATCTTGTAAAACAATTTGCCCATTTTCCATACCCATTATTTCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCg                                                                 <  1:105586/251‑1 (MQ=255)                                                                         gAAAAAATATCTTGTAAAACAATTTGCCCATTTTCCATACCCATTATTTCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGAccagtgtt                                                      >  1:401941/1‑245 (MQ=255)                                                                                              aTTTGCCCATTTTCCATACCCATTATTTCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGaa                                 <  1:145787/251‑1 (MQ=255)                                                                                              aTTTGCCCATTTTCCATACCCATTATTTCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGaa                                 <  1:775030/251‑1 (MQ=255)                                                                                                    ccattttCCATACCCATTATTTCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGAATGCTtc                           <  1:390491/248‑1 (MQ=255)                                                                                                        tttCCATACCCATTATTTCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGAATGCTTCTcca                       >  1:462289/1‑251 (MQ=255)                                                                                                         ttccatacccattatttCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGAATGCTTCTcca                       >  1:676093/17‑250 (MQ=255)                                                                                                          tccatacccattatttCAGTAATATTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCAtc                                                  <  1:737220/207‑2 (MQ=255)                                                                                                                            gtaatatTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGaa                                 <  1:846452/218‑1 (MQ=255)                                                                                                                            gtaatatTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGAATGCTTCTCCACCACGGCATTcacctacgat   >  1:563338/4‑242 (MQ=255)                                                                                                                             taatatTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGAATGCTTCTCCACCACGGCATTcacctacgata  <  1:405504/249‑11 (MQ=255)                                                                                                                             taatatTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGAATGCTTCTCCACCACGGCATTcacctacgat   <  1:629222/248‑10 (MQ=255)                                                                                                                             taatatTCATGATTTTTCTTGTTCCATCATTTAAACGGGAGGCTTGTACGATGATATTTACTGCAGATGCAATATTACGGCGAATTGCTTCCAACGGTAGGGACGCATTCGACATCATAACCATACTCTCCAATCGAGATGTTGCATCCCGAGGGGTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCAGTGTTCATCGCTTGTAACATCTGGAATGCTTCTCCACCACGGCATTcacctacgat   <  1:418153/248‑10 (MQ=255)                                                                                                                                       |                                                                                                                                                                                                                                                 GTTAGTAGTCCGTGATTAATAAACTCACCAAGAATTTTACCCTCTTTATCACGATATTTACTTGCCTTATAGGAAAATATATCCTGTAAAACGATTTGACCGTTTTCCATTCCCATAACCTCAGTGATATTCATAATTTTTCTTGAACCATCATTTAAACGGGACGCCTGCACTATTATGTTTACTGCTGAAGAGATATTACGGCGAATTGCCTCTAACGGTAGGGATGCATTCGACATCATAACCATACTTTCCAAGCGAGATGTTGCATCGCGAGGGCTATTTGCGTGCAACGTTGACATAGAACCATCGTGACCGGTATTCATTGCTTGTAACATTTGGAATGCTTCTCCACCACGGCATTCGCCCACAATA  >  1/855401‑855775

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