Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	851,960	2 bp→AC	intergenic (‑2303/‑70257)	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	851,960	0	G	A	97.1%	81.2 / ‑4.5	34	intergenic (‑2303/‑70258)	unknown/unknown	tRNA‑Met‑CAT/tRNA‑Leu‑GAG
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (12/21);  total (12/22)
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.98e-01
*	1	851,961	0	T	C	100.0%	96.5 / NA	34	intergenic (‑2304/‑70257)	unknown/unknown	tRNA‑Met‑CAT/tRNA‑Leu‑GAG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/22);  total (12/22)

GGTTTGCTGCTTCTGAGCGTGGTGACAAATAAGAGAGATTATTATCCAACTTACGATAAGGCTGGCAGGTACTTAGAGTCGTTAATGATTCAGGACTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGGCTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAATCAGTACTTTTTAGTGGTTTGCTCCCATTATCATTGTATAGTTGATTACGTTCACGAAGAACATTGACCAGTGAGTATGAGGCATTGTCCAGTTTCCCTTGGGTCGTACGAACAATAACTAAGTCTGTTAAAAAAGCGAAAATAAATACTAAAAGTATAAGCATAAAAAGGAACTCAAGACTAACTGCGCCCCTTTTAT  >  1/851723‑852118                                                                                                                                                                                                                                              ||                                                                                                                                                              ggTTAGCTGATTCTGAGCGTGGCGACAAATAAGATAGATTATCATCCAACTTACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACgct                                                                                                                                                     <  1:385261/249‑3 (MQ=255)           ttCTGAGCGTGGCGACAAATAAGATAGATTATCATCCAACTTACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCgaagaatat                                                                                                                                         <  1:727878/251‑4 (MQ=255)                           aaaTAAGATAGATTATCATCCAACTTACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCgaagaatat                                                                                                                                         <  1:225663/235‑4 (MQ=255)                              tAAGATAGATTATCATCCAACTTACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGc                                                                                                                      <  1:794139/251‑1 (MQ=255)                                agatagaTTATCATCCAACTTACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGAc                                                                                                                                     >  1:670121/5‑234 (MQ=255)                                     gaTTATCATCCAACTTACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAGCTGATTACg                                                                                                                                                       <  1:539316/211‑1 (MQ=255)                                          tcatCCAACTTACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGttt                                                                                                          <  1:18778/249‑1 (MQ=255)                                          tcatCCAACTTACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGttt                                                                                                          <  1:675512/249‑1 (MQ=255)                                                 aCTTACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGttttccctgg                                                                                                   >  1:276791/1‑244 (MQ=255)                                                   ttACGATAAGGCCGACATGTGCCTAGATTCGTTAACATTTCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGt                                                                                                 <  1:254629/251‑1 (MQ=255)                                                                    tgtgcctagattcgttaacatttCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAATGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTa                                                                                <  1:720347/230‑1 (MQ=255)                                                                            gattcgttaacatttCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGtt                                                                        >  1:651093/14‑251 (MQ=255)                                                                              ttcgttaacatttCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACt                                                                                 <  1:524724/229‑1 (MQ=255)                                                                                        tttCAGGATTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGaatacg                                                                                                                          <  1:126619/188‑6 (MQ=255)                                                                                                ttttttttCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGCCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAAc                                                                                  <  1:431347/220‑1 (MQ=255)                                                                                                        ctctTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGTTAAAAAGGCGAAAATAAATACTAAAAgaa                                            <  1:779739/251‑3 (MQ=255)                                                                                                        ctctTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGTTAAAAAGGCGAAAATAAATACTAAAAgaa                                            <  1:574359/251‑3 (MQ=255)                                                                                                                   gCTCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTATCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGTTAAAAAGGCGAAAATAAATACTAAAAGAATAAGCATaaaa                                 >  1:442647/1‑251 (MQ=255)                                                                                                                     tCAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCAGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGTTAAAAAGGCGAAAATAAATACTAAAAGAATAAGCATAAAAAg                               <  1:67448/251‑1 (MQ=255)                                                                                                                      cAATAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTc                                                                            <  1:664346/205‑1 (MQ=255)                                                                                                                        aTAGTGACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGTTAAAAAGGCGAAAATAAATAc                                                    >  1:147517/1‑227 (MQ=255)                                                                                                                             gACACCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGt                                                                                                 <  1:682741/177‑1 (MQ=255)                                                                                                                                acCGACTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGTTAAAAAGGCGAAAATAAATACTAAAAGAATAAGCATAAAAAGAAGCTCAAgcg                    <  1:779644/251‑3 (MQ=255)                                                                                                                                    aCTTTATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGt                                                                                                            <  1:568382/159‑1 (MQ=255)                                                                                                                                       ttATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTc                                                                            >  1:892138/1‑188 (MQ=255)                                                                                                                                        tATTTGGACTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGc                                                                                                                      <  1:33869/145‑1 (MQ=255)                                                                                                                                                  gTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGt                                                                                                 <  1:577717/156‑1 (MQ=255)                                                                                                                                                  gTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTAc                                                                                            >  1:360430/1‑161 (MQ=255)                                                                                                                                                  gTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTa                                                                                >  1:589346/1‑173 (MQ=255)                                                                                                                                                  gTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGTTaaaaagg                                                                 <  1:281579/188‑3 (MQ=255)                                                                                                                                                  gTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGTTAAAAAGGCGAAAATAAATACTAAAAGAATAAGCATAAAAAga                              >  1:459000/1‑222 (MQ=255)                                                                                                                                                  gTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAACTAGGTCTGTTAAAAAGGCGAAAATAAATACTAAAAGAATAAGCATAAAAAGAAGCTCAAGCGTAACTGCGCCCTTTttat  >  1:884263/1‑251 (MQ=255)                                                                                                                                                         tctcCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGAATACGAGGCATTGTCCAGTTTTCCCTGGGTAGTACGAACAATAAc                                                                                  >  1:158211/1‑164 (MQ=255)                                                                                                                                                                 caaaATAAGTTTTGCCATTTTTTCATATTCAGTTAAGTCAGTACTTTTTAGTTTTTCGGTCCCATTATCATTGTATAACTGATTACGCTCGCGAAGAATATTGACCAGTGaatacg                                                                                                                          >  1:790982/1‑111 (MQ=255)                                                                                                                                                                                                                                              ||                                                                                                                                                              GGTTTGCTGCTTCTGAGCGTGGTGACAAATAAGAGAGATTATTATCCAACTTACGATAAGGCTGGCAGGTACTTAGAGTCGTTAATGATTCAGGACTTTTTTTCTCTTCCCAGTGCTCAATAGTGACACCGACTTTATTTGGGCTGTCTTTATCTCCAAACAAAATAAGTTTTGCCATTTTTTCATATTCAGTTAAATCAGTACTTTTTAGTGGTTTGCTCCCATTATCATTGTATAGTTGATTACGTTCACGAAGAACATTGACCAGTGAGTATGAGGCATTGTCCAGTTTCCCTTGGGTCGTACGAACAATAACTAAGTCTGTTAAAAAAGCGAAAATAAATACTAAAAGTATAAGCATAAAAAGGAACTCAAGACTAACTGCGCCCCTTTTAT  >  1/851723‑852118

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