Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,012,829	A→T	8.3%	Q414L (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,829	0	A	T	8.3%	138.3 / 3.6	48	Q414L (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 = 9.66e-01

ACAACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGAGCCT  >  NC_002947/1012682‑1012974                                                                                                                                                      |                                                                                                                                                  acaACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcag                                                                                                                                                  <  2:351140/149‑1 (MQ=255)   aaCCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcagca                                                                                                                                                <  6:23287/149‑1 (MQ=255)   aaCCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcag                                                                                                                                                  >  1:40678/1‑147 (MQ=255)   aaCCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcag                                                                                                                                                  >  2:134849/1‑147 (MQ=255)               aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaag                                                                                                                                     >  2:89979/1‑148 (MQ=255)               aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGccttgac                                                                                                                                    >  8:148538/1‑144 (MQ=255)               aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGccttgac                                                                                                                                    <  7:148538/149‑6 (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                                                                           <  3:349455/110‑1 (MQ=255)                                                                                                                 gACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcggt                                                                           >  4:349455/1‑110 (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                                                            <  6:368260/103‑1 (MQ=255)                                                                                                                                       tGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGtt                                                            >  5:368260/1‑103 (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)                                                                                                                                                      |                                                                                                                                                  ACAACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGAGCCT  >  NC_002947/1012682‑1012974

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