Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,129,843	C→A	30.3%	P349P (CCC→CCA)	dmoA‑I →	dimethyl‑sulfide monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,129,843	0	C	A	30.3%	16.7 / 4.1	14	P349P (CCC→CCA)	dmoA‑I	dimethyl‑sulfide monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (3/6);  new base A (2/2);  total (5/9)
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.34e-01

TGGCGCTGATCTCGGGATGGACGGGTATCGACTTCGGTCAGTACCGGCCCGACCAGGTGCTCAAGCATATCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGACCGTCGCCGACGAGCTGCAGGCCTGGGTCGAGGAAACCGACGTGGACGGCTTCAACCTGGCTTAC  >  NC_002947/3129722‑3129990                                                                                                                          |                                                                                                                                                    tggcgctgATCTCGGGATGGACGGGTATCGACTTCGGTCAGTACCGGCCCGACCAGGTGCTCAAGCATATCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCTGAACGAGCAATGGACGGTGCAGGAACTgg                                                                                                                         <  1:19901/150‑1 (MQ=255)            tcGGGATGGACGGGTATCGACTTCGGTCAGTACCGGCCCGACCAGGTGCTCAAGCATATCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCAGACAAGCAATGGACGGTGCAgg                                                                                                                               >  1:132234/1‑133 (MQ=255)            tcGGGATGGACGGGTATCGACTTCGGTCAGTACCGGCCCGACCAGGTGCTCAAGCATATCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCAGACAAGCAATGGACGGTGCAgg                                                                                                                               <  2:132234/133‑1 (MQ=255)            tcGGGATGGACGGGTATCGACTTCGGTCAGTACCGGCCCGACCAGGTGCGCAAGCATATCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACGGCCGACCCCGACAAGCAATGGACGGTGCCGGAACTCGCCGGTTGGGtt                                                                                                              >  2:49060/1‑150 (MQ=255)                                  tCGGTCAGTACCGGCCCGACCAGGTGCTCAAGCATATCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGt                                                                                        <  2:452708/150‑1 (MQ=255)                                                        ggTGCTCAAGCATATCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGAc                                                                  <  2:692004/150‑1 (MQ=255)                                                                      tCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGTTGATCGTCGGCa                                                                            >  1:587712/1‑126 (MQ=255)                                                                      tCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGTTGATCGTCGGCa                                                                            <  2:587712/126‑1 (MQ=255)                                                                         aTACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCAGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGACCGTc                                                              >  1:378161/1‑137 (MQ=255)                                                                         aTACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCAGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGACCGTc                                                              <  2:378161/137‑1 (MQ=255)                                                                                           cGGCAGTGGAAACCTTTTCCACCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCTGCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGACCGTCGCCGACGAGCTGCAGGCCTGGGTCGAGGaaa                               <  2:546668/150‑1 (MQ=255)                                                                                                  tgAAACCCTTTCCACCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGCGGTCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGACCGTCGCCGACGAGCTGCAGGCCTGGGTCGAGGAAACCGACGt                        <  1:49060/149‑1 (MQ=255)                                                                                                        cTTTTCCACCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGGATTGGCGGCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGACCGTCGCCGACGAGCTGCAGGCCTGGGTCGAGGAAACCGACGTGGACgg                  >  1:280833/1‑150 (MQ=255)                                                                                                               aCCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGACCGTCGCCGACGAGCTGCAGGCCTGGGTCGAGGAAACCGACGTGGACGGCTTCAAc           <  2:280833/150‑1 (MQ=255)                                                                                                                        cccGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGACCGTCGCCGACGAGCTGCAGGCCTGGGTCGAGGAAACCGACGTGGACGGCTTCAACCTGGCTTAc  >  2:198250/1‑150 (MQ=255)                                                                                                                          |                                                                                                                                                    TGGCGCTGATCTCGGGATGGACGGGTATCGACTTCGGTCAGTACCGGCCCGACCAGGTGCTCAAGCATATCCATACCAATGCCATTCAGTCGGCAGTGGAAACCTTTTCCACCGCCGACCCCGACAAGCAATGGACGGTGCAGGAACTGGCCGATTGGGTTGGTATTGGCGGCTTCGGCCCGTTGATCGTCGGCAGTGCGCAGACCGTCGCCGACGAGCTGCAGGCCTGGGTCGAGGAAACCGACGTGGACGGCTTCAACCTGGCTTAC  >  NC_002947/3129722‑3129990

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