Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,102,999	C→G	21.9%	noncoding (35/76 nt)	PP_t23 →	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,102,999	0	C	G	21.9%	53.4 / 12.7	32	noncoding (35/76 nt)	PP_t23	tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base C (12/13);  new base G (4/3);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

CTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTCAAATTGCGGGAATAGCTCA  >  NC_002947/2102862‑2103139                                                                                                                                          |                                                                                                                                             cTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCCCGACCTTGGCAAGGTCgtg                                                                                                                                    <  4:219375/148‑3 (MQ=255)          tGAAAAGTGTTGTGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                                                                                                                          >  2:196930/1‑149 (MQ=255)          tGAAAAGTGTTGTGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                                                                                                                          >  5:272711/1‑149 (MQ=255)          tGAAAAGTGTTGTGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                                                                                                                          >  7:48274/1‑149 (MQ=255)          tGAAAAGTGTTGTGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTGTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGAGCTTGGCAAGGTCGGGGTCGCGAGtt                                                                                                                          >  5:96215/1‑149 (MQ=255)                 tgttgtGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                               >  1:36599/1‑137 (MQ=40)                 tgttgtGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                               <  2:36599/137‑1 (MQ=40)                         gtCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCCCGACCTTGGCAAGGTCgtg                                                                                                                                    >  1:7238/1‑122 (MQ=255)                         gtCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCCCGACCTTGGCAAGGTCgtg                                                                                                                                    <  2:7238/124‑3 (MQ=255)                                       gTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTaaa                                                                                              >  6:17250/1‑148 (MQ=37)                                          gCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGATTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTCCGAGTCTCGTTTCCCGCTCCAGTTTAAATACt                                                                                          <  1:92808/149‑1 (MQ=255)                                                       cccTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGc                                                                             <  5:312616/149‑1 (MQ=255)                                                               cggcggTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGc                                                                                     >  3:243745/1‑133 (MQ=255)                                                               cggcggTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGc                                                                                     <  4:243745/133‑1 (MQ=255)                                                                        cAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGaa                                                            <  8:48274/149‑1 (MQ=255)                                                                         aGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAg                                                           >  3:226468/1‑149 (MQ=255)                                                                                  aGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAg                                                  >  4:56884/1‑149 (MQ=255)                                                                                         tAGGGGACTCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGAGTCTATGCGGTGCAGTGGTGGAGAAGGTAGTAAAGGAGCtt                                            >  4:304797/1‑148 (MQ=21)                                                                                           ggggACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTc                                           <  4:226468/147‑1 (MQ=255)                                                                                                  ccAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                       >  6:280423/1‑64 (MQ=35)                                                                                                     aaTGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATAATGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTg                               <  1:196930/149‑1 (MQ=255)                                                                                                      aTGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTgg                              <  3:56884/149‑1 (MQ=255)                                                                                                                  cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGc                        <  5:334770/143‑1 (MQ=255)                                                                                                                  cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGc                        >  6:334770/1‑143 (MQ=255)                                                                                                                         ggTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTc                     <  1:329750/139‑1 (MQ=255)                                                                                                                         ggTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTc                     >  2:329750/1‑139 (MQ=255)                                                                                                                             gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTCAAATTGCGGGAATa       <  5:17250/149‑1 (MQ=255)                                                                                                                             gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTCAAATTGCGGGAATa       <  3:304797/149‑1 (MQ=255)                                                                                                                             gaccccGACCTTGGCAAGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTATATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTCAAATTGCGGGaa         >  6:91092/6‑147 (MQ=255)                                                                                                                                  cGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTCAAATTGCGGGAATa       <  7:88436/144‑1 (MQ=255)                                                                                                                                  cGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTCAAATTGCGGGAATa       >  8:88436/1‑144 (MQ=255)                                                                                                                                  cGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTCAAATTGCGGGAATAGCTCa  <  6:272711/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                             CTGGTGCACTGAAAAGTGTTGTGTGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAGTCCCCTAGGGGACGCCAAATGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAGTTTAAATACTGTGGCGTCTATGCGGTGCAGTGGTGGTGAAGGTCGTAAAGGCGCTTCACGCCGAACGTGGTAAAGCTTCAAATTGCGGGAATAGCTCA  >  NC_002947/2102862‑2103139

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