Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,012,832	A→T	8.3%	Q415L (CAG→CTG)	prfC →	peptide chain release factor 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,012,832	0	A	T	8.3%	132.3 / 3.3	48	Q415L (CAG→CTG)	prfC	peptide chain release factor 3
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base T (2/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

AACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGAGCCTATCACC  >  NC_002947/1012684‑1012980                                                                                                                                                       |                                                                                                                                                     aaCCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcagca                                                                                                                                                      <  6:23287/149‑1 (MQ=255)             aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaag                                                                                                                                           >  2:89979/1‑148 (MQ=255)             aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGccttgac                                                                                                                                          <  7:148538/149‑6 (MQ=255)             aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGccttgac                                                                                                                                          >  8:148538/1‑144 (MQ=255)             aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGc                                                                                                                                                <  8:76439/143‑1 (MQ=255)             aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGc                                                                                                                                                >  7:76439/1‑143 (MQ=255)              tCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGGGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaag                                                                                                                                           >  5:9605/1‑147 (MQ=255)               ccAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGaagaag                                                                                                                                        >  2:197542/1‑149 (MQ=255)                cAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGaagaag                                                                                                                                        <  5:294071/148‑1 (MQ=255)                 aGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGt                                                                                                                                      <  1:134849/149‑1 (MQ=255)                  gATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTg                                                                                                                                     <  3:357731/149‑1 (MQ=255)                       gCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACTCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACt                                                                                                                                <  8:161213/149‑1 (MQ=255)                                  aCCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGttctt                                                                                                                     <  1:197542/149‑1 (MQ=255)                                    cGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTcc                                                                                                                   >  3:201746/1‑149 (MQ=255)                                    cGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTcc                                                                                                                   <  6:360997/149‑1 (MQ=255)                                         gCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGTCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAg                                                                                                              <  8:132627/149‑1 (MQ=255)                                          cGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAgc                                                                                                             >  7:78334/1‑149 (MQ=255)                                                  tGGGCTTCACCGGTATCCCGCACTTCGTCGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGAAACTGCGTCAAGGCCTGCAGCAGCTGGCCGaagaag                                                                                                                                        >  1:16620/1‑115 (MQ=255)                                                  tGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGaagaag                                                                                                                                        <  2:16620/114‑1 (MQ=255)                                                    ggCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAAc                                                                                                     <  7:20116/147‑1 (MQ=255)                                                        tCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAAGAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCACCGAcat                                                                                                >  1:96456/1‑149 (MQ=255)                                                        tCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGAcat                                                                                                >  6:156187/1‑148 (MQ=255)                                                          aCCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGAcatcat                                                                                             >  8:176778/1‑149 (MQ=255)                                                                   ccGCACTTCG‑CGCCGGGGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTg                                                                                     <  6:9605/148‑1 (MQ=255)                                                                   ccGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGc                                                                                    >  5:191276/1‑149 (MQ=255)                                                                   ccGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGc                                                                                    >  8:171549/1‑149 (MQ=255)                                                                       aCTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCACCGACATCATCCTgg                                                                                        <  7:319907/141‑1 (MQ=255)                                                                       aCTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTgg                                                                                        >  8:319907/1‑141 (MQ=255)                                                                           cg‑cgCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcggtcggt                                                                             <  2:40678/148‑1 (MQ=255)                                                                              cgcCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGtgtg                                                                          <  7:171549/149‑1 (MQ=255)                                                                                ccGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGtgct                                                                        <  4:293833/149‑1 (MQ=255)                                                                                  ggAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGtgctgc                                                                      <  5:156187/149‑1 (MQ=255)                                                                                           ccgccgCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGAc                                                              <  2:96456/148‑1 (MQ=255)                                                                                                        cTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTgg                                                                                        >  6:378664/1‑109 (MQ=255)                                                                                                        cTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCg                                                >  8:218717/1‑149 (MQ=255)                                                                                                           aa‑aGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCtt                                             <  8:78334/149‑1 (MQ=255)                                                                                                           aa‑aGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCt                                              <  4:218615/148‑1 (MQ=255)                                                                                                               gACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcggt                                                                                 >  4:349455/1‑110 (MQ=255)                                                                                                               gACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcggt                                                                                 <  3:349455/110‑1 (MQ=255)                                                                                                                   cGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTgaaggaa                                      <  7:218717/149‑1 (MQ=255)                                                                                                                     cTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAaga                                    <  4:201746/149‑1 (MQ=255)                                                                                                                          aTCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTAc                                >  5:192279/1‑148 (MQ=255)                                                                                                                           tCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTAc                                >  3:357985/1‑147 (MQ=255)                                                                                                                             cAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGa                                                               <  7:45752/114‑1 (MQ=255)                                                                                                                             cAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGa                                                               >  8:45752/1‑114 (MQ=255)                                                                                                                              aaGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGt                           >  2:153827/1‑149 (MQ=255)                                                                                                                                     tGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGtt                                                                  >  5:368260/1‑103 (MQ=255)                                                                                                                                     tGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGtt                                                                  <  6:368260/103‑1 (MQ=255)                                                                                                                                             ggCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGa            <  4:357985/149‑1 (MQ=255)                                                                                                                                                   cagcagCTGGCCGAAGAAGGTGCTACGCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGGTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGAGCCt        <  6:185007/147‑1 (MQ=255)                                                                                                                                                    agcagcTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTa                                 >  7:74775/1‑121 (MQ=255)                                                                                                                                                    agcagcTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTa                                 <  8:74775/121‑1 (MQ=255)                                                                                                                                                       agcTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGAGCCTATCAcc  <  7:176778/149‑1 (MQ=255)                                                                                                                                                       agcTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGAGCCTATCAcc  <  1:153827/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     AACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGAGCCTATCACC  >  NC_002947/1012684‑1012980

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