Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,011,852	C→T	15.6%	T87T (ACC→ACT)	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,011,852	0	C	T	15.6%	52.2 / 3.0	26	T87T (ACC→ACT)	prfC	peptide chain release factor 3
Reads supporting (aligned to +/- strand):  ref base C (13/8);  new base T (2/2);  total (15/11)
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.85e-01

TACCGTGAAGTCGCGCAAGTCCGACCGCCACGCCACCTCCGACTGGATGGAAATGGAGAAGCAGCGCGGCATCTCCATCACCACCTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATGGACGTGTGCCGCCTGCGCGACACGCCCATCGTCAG  >  NC_002947/1011717‑1012001                                                                                                                                        |                                                                                                                                                      tACCGTGAAGTCGCGCAAGTCCGACCGCCACGCCACCTCCGACTGGATGGAAATGGAGAAGCAGCGCGGCATCTCCATCACCACCTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGa                                                                                                                                         <  1:3739/150‑1 (MQ=255)                      cGACCGCCACGCCACCTCCGACTGGATGGAAATGGAGAAGCAGCGCGGCATCTCCATCACCACCTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCg                                                                                                                    >  1:621747/1‑150 (MQ=255)                                                                 cgcgGCATCTCCATCACCACCTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGAcggc                                                                         <  1:65837/150‑1 (MQ=255)                                                                            cATCACCACCTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGgtgt                                                              >  2:547435/1‑150 (MQ=255)                                                                            cATCACCACCTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGGGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGgtgt                                                              >  1:607446/1‑150 (MQ=255)                                                                                  caccTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGCCGGCGGTAACGGTGTAGAGcc                                                        >  1:27210/1‑150 (MQ=255)                                                                                     cTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTa                                                                                                                       >  1:14293/1‑84 (MQ=255)                                                                                     cTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTa                                                                                                                       <  2:14293/84‑1 (MQ=255)                                                                                     cTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGgcg                                                                                           <  1:163538/112‑1 (MQ=255)                                                                                     cTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGgcg                                                                                           >  2:163538/1‑112 (MQ=255)                                                                                              gCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGc                                            >  2:635432/1‑150 (MQ=255)                                                                                              gCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGc                                            >  1:586894/1‑150 (MQ=255)                                                                                              gCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACACCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGc                                            <  2:806217/150‑1 (MQ=255)                                                                                                 gTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACTCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGa                                                                                       <  2:216977/104‑1 (MQ=255)                                                                                                 gTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACTCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGa                                                                                       <  2:216981/104‑1 (MQ=255)                                                                                                 gTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACTCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGa                                                                                       >  1:216981/1‑104 (MQ=255)                                                                                                 gTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACTCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGa                                                                                       >  1:216977/1‑104 (MQ=255)                                                                                                 gTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCt                                         >  2:657001/1‑150 (MQ=255)                                                                                                   tctcggaCCGCGAGCACATGATCAACCTGCTCGACACCACCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACTCTGACTGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGa                                       <  2:365625/144‑1 (MQ=255)                                                                                                      cGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATgg                                    <  2:621747/150‑1 (MQ=255)                                                                                                       atACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGCTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATAGa                                   >  1:365625/2‑150 (MQ=255)                                                                                                        tACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATGGAc                                  >  2:1647/1‑150 (MQ=255)                                                                                                                    aTGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATGGACGTGTGCCGCCTg                      >  2:18836/1‑150 (MQ=255)                                                                                                                       aTCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATGGACGTGTGCCGCCTgcgc                   <  1:736909/150‑1 (MQ=255)                                                                                                                       aTCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATGGACGTGTGCCGCCTgcgc                   >  1:135467/1‑150 (MQ=255)                                                                                                                                cTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATGGACGTGTGCCGCCTGCGCGACACGccc          <  2:135467/150‑1 (MQ=255)                                                                                                                                       cccccGGCCACGCAGACCTCTCGGAAGACACCTACCGCACCCTGCCCGCGGTGGACTCGGCGCTGATCGTGCTCGACGGCGGTAAAGGGGTAGAGCCGCGCACCATCGCCCTGATGGACGTGTGCCGCCTGCGCGACACGCCCATCGTCa   <  2:27210/150‑1 (MQ=255)                                                                                                                                        ccccGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATGGACGTGTGCCGCCTGCGCGACACGCCCATCGTCAg  >  2:132606/1‑150 (MQ=255)                                                                                                                                        |                                                                                                                                                      TACCGTGAAGTCGCGCAAGTCCGACCGCCACGCCACCTCCGACTGGATGGAAATGGAGAAGCAGCGCGGCATCTCCATCACCACCTCGGTGATGCAGTTCCCGTACCGCGAGCACATGATCAACCTGCTCGACACCCCCGGCCACGAAGACTTCTCGGAAGACACCTACCGCACCCTGACCGCGGTGGACTCGGCGCTGATGGTGCTCGACGGCGGTAAAGGTGTAGAGCCGCGTACCATCGCCCTGATGGACGTGTGCCGCCTGCGCGACACGCCCATCGTCAG  >  NC_002947/1011717‑1012001

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