Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,171,552	C→T	22.2%	A172V (GCG→GTG)	PP_2783 →	3‑oxoacyl‑ACP reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,171,552	0	C	T	22.2%	28.8 / 5.1	18	A172V (GCG→GTG)	PP_2783	3‑oxoacyl‑ACP reductase
Reads supporting (aligned to +/- strand):  ref base C (8/6);  new base T (2/2);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

CTGCATGATGCGCCAACGCAGCGGGTGCATCGTCAACCTCAGTTCGGTGGCCGCGCAAAAGCCCGGCAAGGGCCAGAGCAACTACGCCGCCGCCAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGAGCAGGAAATCCAGTCGCGGCTGTTGATCAAACGGTTCGCCCGGCCTGAAGAAATTGC  >  NC_002947/3171412‑3171690                                                                                                                                             |                                                                                                                                           cttcATGATGCGCCAACGCAGCGGGTGCATCGTCAACCTCAGTTCGGTGGCCGCGCAAAAGCCCGGCAAGGGCCAGAGCAACTACGCCGCCGCCAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCaa                                                                                                                                   >  2:388115/4‑150 (MQ=255)                                gtaaaCCTCAGTTCGGTGGCCGCGCAAAAGCCCGGCAAGGGCCAGAGCAACTACGCCGCCGCCAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATc                                                                                                    >  1:795946/4‑150 (MQ=255)                                                                 ggccaaagCCAGCCCAACCAAGCCGCCTCCAAAGGCGGTTTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGgc                                                                   >  2:49058/8‑150 (MQ=255)                                                                    aaGGGCCAGAGCAACTACGCCGCCGCCAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGcc                                                                >  2:301890/1‑150 (MQ=255)                                                                         ccAGAGCAACTACGCCGCCGCCAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGa                                                           <  2:795946/150‑1 (MQ=255)                                                                                   tACGCCGCCGCCAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGAGCAGGAAATc                                                 <  2:579039/150‑1 (MQ=255)                                                                                             ccAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGTGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCgg                                                                                                        <  1:380222/85‑1 (MQ=255)                                                                                             ccAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGTGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCgg                                                                                                        <  1:380230/85‑1 (MQ=255)                                                                                             ccAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGTGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCgg                                                                                                        >  2:380222/1‑85 (MQ=255)                                                                                             ccAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGTGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCgg                                                                                                        >  2:380230/1‑85 (MQ=255)                                                                                                    cGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGAGCAGGAAATCCAGTCGCGGCTGTTGAt                                <  1:388115/150‑1 (MQ=255)                                                                                                              gCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGAGCAGGAAATCCAGTCGCGGCTGTTGATCAAACGGTTc                      <  1:301890/150‑1 (MQ=255)                                                                                                                 ttGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGCTGGAGCACGGCATCGTACGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGAGCAGGAAATCCAGTCGCGGCTGTTGATCAAACGGTTCGcc                   >  1:297393/1‑150 (MQ=255)                                                                                                                     cacGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGAGCAGGAAATCCAGTCGCGGCTGTTGATCAAACGGTTCGCCCGGc               <  1:373106/150‑1 (MQ=255)                                                                                                                        gcgcACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGAGCAGGAAATCCAGTCGCGGCTGTTGATCAAACGGTTCGCCCGGCCTg            >  1:373052/1‑150 (MQ=255)                                                                                                                        gcgcACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCAGCATCGTCAGCACCGACATGAGCCAAGCCCTg                                                                         >  1:671831/1‑89 (MQ=255)                                                                                                                        gcgcACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCAGCATCGTCAGCACCGACATGAGCCAAGCCCTg                                                                         <  2:671831/89‑1 (MQ=255)                                                                                                                                  ggTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGAGCAGGAAATCCAGTCGCGGCTGTTGATCAAACGGTTCGCCCGGCCTGAAGAAATTGc  >  1:159174/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                           CTGCATGATGCGCCAACGCAGCGGGTGCATCGTCAACCTCAGTTCGGTGGCCGCGCAAAAGCCCGGCAAGGGCCAGAGCAACTACGCCGCCGCCAAAGGCGGTGTAGAAGCATTGACACGCGCACTGGCGGTGGAGTTGGCGCCGCGCAACATCCGGGTCAACGCGGTGGCGCCCGGCATCGTCAGCACCGACATGAGCCAAGCCCTGGTCGGCGCCCATGAGCAGGAAATCCAGTCGCGGCTGTTGATCAAACGGTTCGCCCGGCCTGAAGAAATTGC  >  NC_002947/3171412‑3171690

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