Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	696,470	C→T	noncoding (35/75 nt)	glnX ←	tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	696,470	0	C	T	90.3%	80.6 / NA	31	noncoding (35/75 nt)	glnX	tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/13);  minor base G (0/3);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGT  >  NC_000913/696226‑696709                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                cGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGAAGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGGACGAGGATTCGAACCTCGGAATGACGGAATTAGAATc                                                                                                                                                                                                                                           >  2:579308/1‑251 (MQ=255)       ggCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGc                                                                                                                                                                                                                                     >  2:274727/1‑251 (MQ=255)                     tCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACACCATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGACTCCGGTGCCTTACCGCTTGGCg                                                                                                                                                                                                                       >  2:552514/1‑251 (MQ=255)                       aGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGa                                                                                                                                                                                                                      <  2:17787/250‑1 (MQ=255)                                                       tACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCg                                                                                                                                                                                      <  2:192096/250‑1 (MQ=255)                                                                gCCGGATGCGTGCCAACCGTAGGGCGGATAAGACTCGGCAGCGTCGGTTCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATa                                                                                                                                                                             <  2:148572/250‑1 (MQ=255)                                                                  cGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAAt                                                                                                                                                                          <  2:328049/251‑1 (MQ=255)                                                                   ggATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGACGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGAGTGCATTTACAATGACAGCGTCTAGACTAAAtt                                                                                                                                                                         >  2:421786/1‑251 (MQ=38)                                                                                    aGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGAt                                                                                                                                                        >  1:321352/1‑251 (MQ=255)                                                                                              aGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCg                                                                                                                                              >  1:569649/1‑251 (MQ=255)                                                                                              aGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGGGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACGGGGTGCACTTACAAGTGAAGCGTCTTGAATAAATTGTCTGGGGGACGAGGATTCGAACCTCg                                                                                                                                              <  2:558081/251‑1 (MQ=38)                                                                                                     gCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAACGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATACGGTGCCTTACCGCTTGGCGATACCCCAACTAGGTGCAATTACAAGGTAAGAGTCTTGAATACATTGGCTGGGGTACGAAGATTTGGACCTCGGAAAGcc                                                                                                                                       >  2:452277/1‑251 (MQ=38)                                                                                                          gTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGCCATTGCCTGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTATTAATGTGGCTGGGGTAAGAGGATTCGAACATCGGAATGCCGGCATTAGAATCCGGTGCATTACCGCTTGGCGATACCCCAACTGGGTGCACTTCCAAGGTAAGTGTCTTGAATAAATTGTCTGGGGTACGAGGATACGAACCTAGGAATGCCGGAAt                                                                                                                                  >  1:351591/1‑251 (MQ=255)                                                                                                               aTCCGACATTGAAGGAAAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTATGCGTCTTGAATAAATTGGCTGGGTTATGAGGCTTCGCACCTCGCAATGCCGGTATCAGa                                                                                                                              >  1:7074/1‑250 (MQ=38)                                                                                                                  cGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCATCTGGGTGAAATTACAAGGTAAGCGTATTGAATAAATTGGCTGGGGTACGCGGATTTGAACCGCGGAATGCCGGAATCAGAATcc                                                                                                                          >  1:116685/1‑251 (MQ=38)                                                                                                                                  gACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGAATTAGAAACCGGGGCCTTACCGCTTGGCGATACCCCACCTGTCAGAACTTACAATGGAATCGTCTTGACTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGACTCCGGTGCCTTACCGCTTg                                                                                                          >  1:380100/1‑251 (MQ=37)                                                                                                                                              gcGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaa                                                                                              <  1:360934/250‑1 (MQ=255)                                                                                                                                               cGCATTCGGTATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCTGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaa                                                                                              <  1:579308/250‑1 (MQ=255)                                                                                                                                                  aTTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACaaa                                                                                          <  1:470287/251‑1 (MQ=255)                                                                                                                                                    tcgttaTTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTGTTTAATGTGGCTGGGGTACGAGGAATCGAACCTCGGAATGGCGGAATGAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTTTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAGGCCTGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                                                                                         <  1:290437/245‑1 (MQ=38)                                                                                                                                                     gggcttggaatgaaCGCAGAAAAGCAAAAAGCACGCCGAAGCGAGCTTGTTTAATGGGGCTGGGGCAGGAGGACTCGAAACTCGGAATGACAGGACGAGAAACCGGTGCCTTACCGATTGGAGGTACGCCGACTGGGTGCACGTACGAGGTAGGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTg                                                                                       <  2:156560/244‑1 (MQ=21)                                                                                                                                                              atgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTTGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATCCCCCAACAACTTGGTTTTGAAt                                                                              >  1:20271/1‑251 (MQ=255)                                                                                                                                                                gaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGAGGCTGGGGTACGAAGATTCGAACATCGGACTGCCGGAATTAGAATCCGGCGCCTTCCCGCTTGGCGATAACCAAACTGGGAGCCCTTACCAGGTAAGCGCATTGAATAAATTACCTGGGGTACCAGGAATCGAACCTCGGAATGACGGAGTCAGAATCCGGTGCCTTACCGCTTGCCGATACCCCAACTAAtcagatttaaattt                                                                            >  1:94091/1‑238 (MQ=255)                                                                                                                                                                            gaaaacaGCTCGCCGAAGCAATGTTTGTTAAAGTGGCTGGGGTACGAGTGCTCCAGACACGGAATGCTGGGATTAGAACCCAGCGCCTTTCCGCTTGACTATACCCCGAATGGGTGCACTCACAAGGTATGCGTCTTGAATAAATTGGCTGGGGTACTAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATAtt                                                                <  2:552984/245‑1 (MQ=16)                                                                                                                                                                             cAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGACTCGAACCTCGGAATGCCGGAATGAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTc                                                               <  1:421786/251‑1 (MQ=255)                                                                                                                                                                                      cgcggaagagagatttgtgaattgggCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCAGTGCATTACCGCTTGGCTATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAAGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTc                                                      <  2:323964/227‑1 (MQ=37)                                                                                                                                                                                                                    ggTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAAc                        <  1:552514/251‑1 (MQ=255)                                                                                                                                                                                                                      tACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCCAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCt                      >  1:262805/1‑251 (MQ=255)                                                                                                                                                                                                                       aCGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCtg                     <  1:151142/251‑1 (MQ=255)                                                                                                                                                                                                                                    ccTCGGAATGCCGGAATTAGAATCCAGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCTCTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGTCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCGGTGACCCCATCAtt        <  1:95946/251‑1 (MQ=255)                                                                                                                                                                                                                                          aaTGCCGGAATTCGAATCCGGAGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGCACCTGTTACCCCTTCATTCTGAGt  >  1:425081/1‑251 (MQ=38)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                CGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGT  >  NC_000913/696226‑696709

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