Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,820,944	C→A	23.5%	P236Q (CCA→CAA)	PP_1624 →	group II intron‑encoding maturase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,820,944	0	C	A	23.5%	23.7 / 5.0	17	P236Q (CCA→CAA)	PP_1624	group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base C (5/8);  new base A (2/2);  total (7/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 = 5.60e-01

TGGAACTCGATCTGGAGAAGTTTTTCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCAGGAGGGGACGCCGCAAGGCGGCCCGCTCTCGCCGTTGCTGTCGAA  >  NC_002947/1820797‑1820995                                                                                                                                                    |                                                    tGGAACTCGATCTGGAGAAGTTTTTCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCA                                                    >  3:139559‑M1/1‑149 (MQ=255)   gAACTCGATCTGGAGAAGTTTTTCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACG                                                  <  1:97044‑M1/149‑1 (MQ=255)   gAACTCGATCTGGAGAAGTTTTTCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCAC                                                   <  1:272009‑M1/148‑1 (MQ=255)   gAACTCGATCTGGAGAAGTTTTTCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCAC                                                   <  4:311952‑M1/148‑1 (MQ=255)     aCTCGATCTGGAGAAGTTTTTCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACGG                                                 >  3:277307‑M1/1‑148 (MQ=255)      cTCGATCTGGAGAAGTTTTTCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACGG                                                 >  3:30894‑M1/1‑147 (MQ=255)           tCTGGAGAAGTTTTTCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGACAG‑CAACGGCgagagcg                                         >  7:61351‑M1/1‑142 (MQ=255)                     tttttCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAAAGACAAGTGCGTACTCAGGCTTATTCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCAACGGCgaaagcgggccgcctt                                >  5:303309‑M1/1‑133 (MQ=255)                                          gACATCCTGATGGCCTGTATCGAGCGCCGGGTCAAAGATAAGTGCGTACTCAGGCTTATTCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCgaaagcgggccgccttgcggcgtcccctcctgccgtg           <  5:364703‑M1/149‑38 (MQ=255)                                            cATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACG                                                  >  7:302663‑M1/1‑108 (MQ=255)                                            cATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACG                                                  <  8:302663‑M1/108‑1 (MQ=255)                                            cATCCTCATGGCCTGCATCGAACGTCGGAGCAAAGACAAGTGCGTACTCAGGCTTATTCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCgaaagcgggccgccttgcggcggcccctcctgccgtggg         <  2:145222‑M1/149‑40 (MQ=255)                                                cTCATGGCCTGCATCGAACGTCGGATCAAAGACAAGTGCGTACTCAGGCTTATTCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCAACGGCgaaagcgggccgccttgcggcgtcccctcctgccgtgggctga     <  6:303309‑M1/149‑44 (MQ=255)                                                   aTGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGACAG‑CAACGGCgagagcgggccgccttgcggcgtcccctcctgccgtgggctgacga  <  8:61351‑M1/148‑47 (MQ=255)                                                                              gaCAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCgagagcgggccgccttgcggcgtcccctcctgccgtg           >  5:340438‑M1/1‑76 (MQ=255)                                                                              gaCAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCgagagcgggccgccttgcggcgtcccctcctgccgtg           <  6:340438‑M1/113‑38 (MQ=255)                                                                                           tCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCaggaggggccgccttgcggcgtcccctc                    <  2:234930‑M1/91‑29 (MQ=255)                                                                                                                                                    |                                                    TGGAACTCGATCTGGAGAAGTTTTTCGATCGGGTCAATCACGACATCCTCATGGCCTGCATCGAACGTCGGATCAGAGACAAGTGCGTACTCAGGCTTATCCGCCGCTACCTTGAAGCCGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCAGGAGGGGACGCCGCAAGGCGGCCCGCTCTCGCCGTTGCTGTCGAA  >  NC_002947/1820797‑1820995

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