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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	696,592	0	T	C	97.9%	157.5 / ‑5.0	47	noncoding (25/75 nt)	glnV	tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (21/25);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCA  >  NC_000913/696479‑696724                                                                                                                  |                                                                                                                                     gTGCCTTAACGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCcaac                                                                                                            >  1:493020/1‑140 (MQ=255)    ccTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACaaa                                                                                                         <  1:514879/140‑1 (MQ=255)        aCCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCcaa                                                                                                             <  1:316780/132‑1 (MQ=255)        aCCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCcaa                                                                                                             >  2:316780/1‑132 (MQ=255)           gCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGttt                                                                                                  >  2:45730/1‑140 (MQ=255)            cTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGtttt                                                                                                 >  2:154952/1‑140 (MQ=255)            cTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGtttt                                                                                                 >  1:363205/1‑140 (MQ=255)                    tACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGc                                                                                         <  2:712729/140‑1 (MQ=255)                       cccAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGa                                                                                      <  2:760537/140‑1 (MQ=255)                         cAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAAc                                                                                    <  1:45730/140‑1 (MQ=255)                           aCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACa                                                                                   >  1:402801/1‑139 (MQ=255)                           aCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACa                                                                                   <  2:402801/139‑1 (MQ=255)                             tGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACatat                                                                                <  2:271749/140‑1 (MQ=255)                               ggTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTc                                                                              >  1:124944/1‑140 (MQ=255)                                gTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCg                                                                             <  2:179241/140‑1 (MQ=255)                                 tGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGa                                                                            >  2:573311/1‑140 (MQ=255)                                 tGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGa                                                                            <  1:98080/140‑1 (MQ=255)                                 tGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGa                                                                            >  2:98080/1‑140 (MQ=255)                                      ttACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGAt                                                                           >  1:322011/1‑136 (MQ=255)                                      ttACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGAt                                                                           <  2:322011/136‑1 (MQ=255)                                        aCAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTc                                                                     >  2:97934/1‑140 (MQ=255)                                              tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAAtt                                                               >  2:638324/1‑140 (MQ=255)                                              tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAAtt                                                               >  1:39945/1‑140 (MQ=255)                                                aGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTTGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTtg                                                             <  2:493020/140‑1 (MQ=255)                                                aGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTtg                                                             >  2:99818/1‑140 (MQ=255)                                                  cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTGGCAGAACATATTCGATACATTCAGAATTtggt                                                           >  2:560896/1‑140 (MQ=255)                                                    tCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTtggtgg                                                         >  2:525608/1‑140 (MQ=255)                                                    tCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTtggtgg                                                         >  2:86221/1‑140 (MQ=255)                                                          aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTacga                                                   >  2:494081/1‑140 (MQ=255)                                                              aTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACggg                                               >  1:519967/1‑140 (MQ=255)                                                                  gCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTc                                           <  2:519967/140‑1 (MQ=255)                                                                   cTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCg                                          <  1:99818/140‑1 (MQ=255)                                                                   cTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCg                                          <  1:494081/140‑1 (MQ=255)                                                                   cTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCg                                          <  1:573311/140‑1 (MQ=255)                                                                   cTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCg                                          <  1:341642/140‑1 (MQ=255)                                                                   cTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCg                                          <  2:124944/140‑1 (MQ=255)                                                                              ggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGAcc                               <  1:638324/140‑1 (MQ=255)                                                                                aTTCGAACCTCGGAATGCCGGAATCCGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATATGGTGGCTACGACGGGAGTCGAACCTGTGAcccc                             >  2:603777/1‑140 (MQ=255)                                                                                aTTCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGAcccc                             <  1:86221/140‑1 (MQ=255)                                                                                  tCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGAcc                               <  2:497846/136‑1 (MQ=255)                                                                                  tCGAACCTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGAcc                               >  1:497846/1‑136 (MQ=255)                                                                                        cTCGGAATGCCGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTAt                     <  1:525608/140‑1 (MQ=255)                                                                                                  cGGAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGAtgtg           <  1:97934/140‑1 (MQ=255)                                                                                                    gAATCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGct         <  1:326169/140‑1 (MQ=255)                                                                                                     aaTCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGAt              >  2:582704/1‑134 (MQ=255)                                                                                                     aaTCAGAATCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGAt              <  1:582704/134‑1 (MQ=255)                                                                                                           aaTCCGGCGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCa  <  2:691276/140‑1 (MQ=255)                                                                                                                  |                                                                                                                                     GTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCA  >  NC_000913/696479‑696724

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