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.9%	62.2 / NA	33	noncoding (35/75 nt)	glnX	tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/11);  minor base G (1/2);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCC  >  NC_000913/696345‑696595                                                                                                                              |                                                                                                                              tGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGAAATCGCGACGAGGATTCGAACCTCGGAATGCCGGAATTAGa                                                                                                                            >  2:111894/1‑129 (MQ=38)       aTAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAATCTCGCCTAAGCTACCTTTTTTAATGTTTCTTTGTTACGATTATTCGAACATCGGAATGCCGGAATGAGAATCCgg                                                                                                                      <  1:89407/129‑1 (MQ=255)         aaGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCATAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTTGGTTTTGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTg                                                                                                                    <  1:419042/129‑1 (MQ=38)         aaGACGTGTCAACATCGCATTCGACAGTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGACTCCggag                                                                                                                    >  2:43831/1‑127 (MQ=255)                  cAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCCAAGCGCGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAACGAGAATCCGGTGCCTTACCGc                                                                                                           >  2:178190/1‑129 (MQ=255)                   aaCATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCt                                                                                                          >  1:15994/1‑129 (MQ=255)                       tCGCATTCGACACTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGTGTTACTTGTTTTCTATTCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGc                                                                                                      <  1:64495/129‑1 (MQ=25)                         gCATTCGACTCTGACTGAACTCAGAATAGCAAGCAGCTCGCCGAAGCGAGCTTTTTTAAGGTGGCTGGGGTACGAGTATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGa                                                                                                    <  2:346383/129‑1 (MQ=18)                          cATTTGACATTGAACGAACGCCGAAAAGCAAAGAGCTCGCCTAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTATAATCCGGTGCCTTACCGCTTGGCGAt                                                                                                   <  2:651707/129‑1 (MQ=21)                               gACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTTACTTACCGCTTGGCGATAcccc                                                                                              >  2:258480/1‑129 (MQ=255)                                 cATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACTTCGTAATGTCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCaa                                                                                            <  1:617226/129‑1 (MQ=255)                                 cATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAACTTTTTAAATGTGACAGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCGTACCGCTTGGCGATACCCCaa                                                                                            >  1:41561/1‑129 (MQ=21)                                    taaattaaCGCAGAAAAGCAAAAAGCTCGTCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTg                                                                                         <  2:83311/127‑1 (MQ=255)                                       atgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGt                                                                                      >  2:29422/1‑129 (MQ=255)                                                       aaaaaGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCGTCGGGTCAATGCCATACCCCAACAGGGTGCACTT‑CAAgggcaag                                                                     >  2:438206/1‑124 (MQ=18)                                                        aaaaGCTCGCCGAAGCGAGCTTGTTTAATGTTGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGTTTTTCCGTTTGGTGATATCCCAACTGGGTGCACTTACAAGGTAAGc                                                                     <  1:27978/129‑1 (MQ=21)                                                             cTCGCCGAAGCGAGCTTTTTTAATGTGGCTAAGAAACAAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCtt                                                                >  1:506257/1‑129 (MQ=25)                                                                   gAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATAAAAACATAGAAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGTTAAGAGTCTTgaagaa                                                          >  1:135280/1‑129 (MQ=21)                                                                     aGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGTCTTACTGTTTTGCGATACCCCAACTGTGTGCACTTACAAGGTAAGCGTCTGGAATAAAt                                                        <  1:402287/129‑1 (MQ=21)                                                                            ttttttAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTAAATATTAAGAGCATAGAATAAATTGGCTgg                                                 >  2:37285/1‑129 (MQ=17)                                                                            ttttttAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGCTGCACTTACAAGGTCAGCGTCATGAATATATTGGGTgg                                                 >  2:343385/1‑129 (MQ=18)                                                                                        gCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAAAATCCGATGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTATTGAATAAATTGGCTGGGGTACGAGGGtt                                     >  1:503288/1‑129 (MQ=21)                                                                                          tGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCg                                   <  2:15994/129‑1 (MQ=255)                                                                                            gggTACGAGGCTTCGAACCTCGGAATGCCGGAATTAGAATAAGAAGACATACAGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTCAGCGTCTTGAATAAAATGGCTGGGGTACGAGGATTCGaa                                 >  1:203801/1‑129 (MQ=25)                                                                                                 cGAGGATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGGACTTAAAAGGTAAGCGTCTTGAATAAATTGGCTGGAGTACGAGGATTCGAACCTCg                            >  1:562861/1‑129 (MQ=25)                                                                                                    ggATTCGAACCTCGGAATGACGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGaa                         <  1:111894/129‑1 (MQ=255)                                                                                                     gATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGAGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                        >  1:51916/1‑129 (MQ=255)                                                                                                     gATTCGAACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTAGGCGATACCACAACAAGCAACACATAAAAGCTAAGCTTCGAGAATAAATTGGCTGGGGTACCACGATTCGAACCCCGGAAt                        >  1:236846/1‑129 (MQ=12)                                                                                                        tcacaCCTCGGAATGCCGGAATGAGAATCCGGTGCCTTCCCGCTTGGCTATACGCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGcc                     <  2:570239/125‑1 (MQ=38)                                                                                                                cGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAGGCTGCACTTAAAAGGTCAGCGTCTAGACTAAACAGGCAGGGGTACGAGGCTACGACCGTAGGAGAACCGGAATCAg             >  2:177665/1‑129 (MQ=12)                                                                                                                     tGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGTTAAGCGTCTGGAATAAATTGGCTGGAGTACGAGGATGCGAACCTAGGAATAAAGGAAACAGAATcc        >  2:143022/1‑129 (MQ=25)                                                                                                                         ggAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGTATGCATTTACATGGTAAGCGTCGTCACTAAATCGGCCGAGGTACGAGCACTCGACCGTCGGCAACCCCAAATCCGAATcccgac    >  2:364253/1‑125 (MQ=12)                                                                                                                           aaTTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAAATGGGTGCACTTACAAGGTAAGCGTCTTTAATAAATTGGCTTGGGTACGATGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc  <  1:41102/129‑1 (MQ=37)                                                                                                                              |                                                                                                                              TGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCC  >  NC_000913/696345‑696595

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