Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,606,502	G→C	7.4%	noncoding (33/74 nt)	PP_t64 ←	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,606,502	0	G	C	7.4%	167.6 / 3.9	54	noncoding (33/74 nt)	PP_t64	tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (21/29);  new base C (2/2);  total (23/31)
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.74e-01

TAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTGGCTACAAGGCGCCGA  >  NC_002947/1606356‑1606640                                                                                                                                                   |                                                                                                                                             tAATTGGCCATGGGCCAGGTCAAGGCAGGATTAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGTTGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGaa                                                                                                                                            <  5:10385/149‑1 (MQ=255)       gCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCttt                                                                                                                                       <  1:315262/148‑1 (MQ=255)       gCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTCCCAAGCttta                                                                                                                                      >  3:257923/1‑148 (MQ=255)              gCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGc                                                                                                                               >  6:334280/1‑149 (MQ=255)              gCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGc                                                                                                                               <  5:69190/149‑1 (MQ=255)              gCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGc                                                                                                                               <  2:274599/149‑1 (MQ=255)              gCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGc                                                                                                                               <  2:220008/149‑1 (MQ=255)              gCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTCCCAAGCTTTAGttgagg                                                                                                                               <  4:257923/149‑5 (MQ=255)                   gTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTACTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACt                                                                                                                           >  2:353837/1‑148 (MQ=255)                    tCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTaa                                                                                                                         <  1:199251/149‑1 (MQ=255)                           aGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATAc                                                                                                                  >  3:115233/1‑149 (MQ=255)                             gATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATAcc‑c                                                                                                               <  8:389849/149‑1 (MQ=255)                                  aaaaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATc                                                                                                          >  6:13764/1‑149 (MQ=255)                                   aaaCGAAAAAAGGAGCCAAGGCTTCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCt                                                                                                         <  1:353837/149‑1 (MQ=255)                                   aaaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCt                                                                                                         <  5:334280/149‑1 (MQ=255)                                   aaaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCt                                                                                                         <  5:45370/149‑1 (MQ=255)                                     aCGTAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGa                                                                                                       <  8:295276/149‑1 (MQ=255)                                              ggAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATc                                                                                                          <  1:120500/137‑1 (MQ=255)                                              ggAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATc                                                                                                          >  2:120500/1‑137 (MQ=255)                                                aGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCAt                                                                                                           >  8:161944/1‑134 (MQ=255)                                                aGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCAt                                                                                                           <  7:161944/134‑1 (MQ=255)                                                    aaGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTAc                                                                                        >  2:192926/1‑148 (MQ=255)                                                     aGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTAc                                                                                        >  6:354566/1‑147 (MQ=255)                                                       gCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGa                                                                                    <  1:375189/149‑1 (MQ=255)                                                            ttttGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCtt                                                                               >  6:16238/1‑149 (MQ=255)                                                               tGACTGCGATGACGTTGCTGGTACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCa                                                                             <  8:85785/148‑1 (MQ=255)                                                               tGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCa                                                                             <  5:45185/148‑1 (MQ=255)                                                                     cGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGaa                                                                        <  5:13764/147‑1 (MQ=255)                                                                         gACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAg                                                                    >  7:357672/1‑147 (MQ=255)                                                                          aCGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGaa                                                                  <  6:9417/148‑1 (MQ=255)                                                                                cTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCt                                                                                                    <  4:306872/109‑1 (MQ=255)                                                                                cTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAg                                                           <  4:115233/149‑1 (MQ=255)                                                                                 tGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGt                                                          >  4:320673/1‑149 (MQ=255)                                                                                 tGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGt                                                          >  3:234170/1‑149 (MQ=255)                                                                                      cGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATaa                                                     <  1:192926/149‑1 (MQ=255)                                                                                           tGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCTCGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATaaaaaa                                                 <  1:328283/149‑1 (MQ=255)                                                                                             cGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCt                                              <  3:320673/149‑1 (MQ=255)                                                                                                       aTGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGTAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGcc                                    >  1:290910/1‑149 (MQ=255)                                                                                                            gAGCGGGTAGAGGTAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATaaaaaa                                                 <  7:210257/131‑1 (MQ=255)                                                                                                            gAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTTTAAAAAAGCTTGTCGGCGCCCCGTg                               >  5:258264/1‑149 (MQ=255)                                                                                                            gAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATaaaaaa                                                 >  8:210257/1‑131 (MQ=255)                                                                                                              gCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTgg                              <  4:234170/148‑1 (MQ=255)                                                                                                                 ggTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCAACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTg                          >  7:71994/1‑149 (MQ=255)                                                                                                                       gggAATCGCACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTa                    >  7:131607/1‑149 (MQ=255)                                                                                                                               aaCCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTGGCTa             >  3:101240/1‑148 (MQ=255)                                                                                                                                  ccTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTGGCTACAAg         >  3:227118/1‑149 (MQ=255)                                                                                                                                    tCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGGTCAGACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGACGCCGCCCCGTGGCCTGCGGGTATTGGCTACAAgg        >  1:112764/1‑149 (MQ=255)                                                                                                                                    aaaaCGAAAGCTTCCCAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTg                 >  7:19029/3‑139 (MQ=255)                                                                                                                                    aaaaCGAAAGCTTCCCAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTg                 <  8:19029/137‑1 (MQ=255)                                                                                                                                      aaCTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCTGGTATTGGCTACAAGgc       <  8:71994/147‑1 (MQ=255)                                                                                                                                      aaCTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTGGCTACAAGgcgc     <  2:290910/149‑1 (MQ=255)                                                                                                                                      aaCTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTGGCTACAAGgcgc     >  1:192306/1‑149 (MQ=255)                                                                                                                                       aCTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTGGCTACAAGGCGcc    <  8:357672/149‑1 (MQ=255)                                                                                                                                           aaGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTGGCTACAAGGCGCCGa  <  5:354566/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             TAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACC‑CGCATCTGAGCTC‑GACTTTTTACCAGAACCTTCAGCGAATTAGAAGCCCAAGTTATAAAAAAGCTTGTCGGCGCCCCGTGGCCTGCGGGTATTGGCTACAAGGCGCCGA  >  NC_002947/1606356‑1606640

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: ‑