Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,737,772	T→C	18.2%	V706A (GTG→GCG)	polB →	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,737,772	0	T	C	18.2%	45.4 / 5.3	22	V706A (GTG→GCG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (10/8);  new base C (2/2);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

TTATGTCCGCGAGTATGTACGGCGCACGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATCACCACCGCCGGGCCTGAGCCGCTGGAAAACCTGCAGGCACC  >  NC_002947/2737656‑2737901                                                                                                                     |                                                                                                                                  ttATGTCCGCGAGTATGTACGGCGCACGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGcc                                                                                                                >  2:396056/1‑136 (MQ=255) ttATGTCCGCGAGTATGTACGGCGCACGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGcc                                                                                                                >  2:393655/1‑136 (MQ=255)         gcgAGTATGTACGGCGCACGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTa                                                                                                        >  2:250873/1‑136 (MQ=255)             gTATGTACGGCGCACGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCcgac                                                                                                             >  1:229476/1‑127 (MQ=255)             gTATGTACGGCGCACGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCcgac                                                                                                             <  2:229476/127‑1 (MQ=255)             gTATGTACGGCGCACGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTacaac                                                                                                    >  2:94750/1‑136 (MQ=255)                     ggCGCACGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCt                                                                                            <  1:477164/136‑1 (MQ=255)                          aCGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCcgac                                                                                                             >  2:457615/1‑114 (MQ=255)                          aCGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCcgac                                                                                                             <  1:457615/114‑1 (MQ=255)                          aCGTTGGCGGGCAAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGCGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAgc                                                                                       >  2:112192/1‑136 (MQ=255)                          aCGTTGGCGGGCAAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGCGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAgc                                                                                       >  2:114726/1‑136 (MQ=255)                            gTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAgc                                                                                       <  1:229029/134‑1 (MQ=255)                            gTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAgc                                                                                       >  2:229029/1‑134 (MQ=255)                                       aGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTa                                                                          <  1:383116/136‑1 (MQ=255)                                                                           gTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATcaccac                                      <  1:94750/136‑1 (MQ=255)                                                                           gTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGCGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATcaccac                                      <  1:112192/136‑1 (MQ=255)                                                                           gTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGCGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATcaccac                                      <  1:114726/136‑1 (MQ=255)                                                                              gCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATCACCAccgc                                   >  2:1246/1‑136 (MQ=255)                                                                                  gctggctgATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATCACCACCGCCggg                               >  1:242930/1‑136 (MQ=255)                                                                                                 gcgcAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATCACCACCGCCGGGCCTGAGCCGCTGGaa                <  1:194219/136‑1 (MQ=255)                                                                                                               cgcATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATCACCACCGCCGGGCCTGAGCCGCTGGAAAACCTGCAGGCAcc  >  1:85817/1‑136 (MQ=255)                                                                                                               cgcATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATCACCACCGCCGGGCCTGAGCCGCTGGAAAACCTGCAGGCAcc  <  1:294090/136‑1 (MQ=255)                                                                                                                     |                                                                                                                                  TTATGTCCGCGAGTATGTACGGCGCACGTTGGCGGGCGAGCAGGACGAACTGCTGGTTTATCGCAAGCGCCTGCGTCGCCCGCTGGCTGATTACCAGCGCAACGTGCCGCCGCATGTGCGTGCCGCACGCCTTGCCGACGACTACAACAAGCGTCTGCAGCGCCCTCTGCAGTACCAGCGTGGGGGCTGGATCAGTTACGTGATCACCACCGCCGGGCCTGAGCCGCTGGAAAACCTGCAGGCACC  >  NC_002947/2737656‑2737901

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