Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,012,835	T→A	7.5%	L416Q (CTG→CAG)	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,835	0	T	A	7.5%	153.7 / 3.0	53	L416Q (CTG→CAG)	prfC	peptide chain release factor 3
Reads supporting (aligned to +/- strand):  ref base T (21/28);  new base A (2/2);  total (23/30)
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.95e-01

ATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGAGCCTATCACC  >  NC_002947/1012696‑1012980                                                                                                                                              |                                                                                                                                                  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                                                                                                                   <  6:360997/149‑1 (MQ=255)                        cGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTcc                                                                                                                   >  3:201746/1‑149 (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                                                                  >  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)                                                                                                                                              |                                                                                                                                                  ATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCG‑CGCCGGAGCTGTTCCGCCGCGTACGCCTGAA‑AGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGAGCCTATCACC  >  NC_002947/1012696‑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: ‑