Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,097,674	G→T	22.2%	L16M (CTG→ATG)	PP_3605 ←	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,097,674	0	G	T	22.2%	33.6 / 4.8	18	L16M (CTG→ATG)	PP_3605	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (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 = 9.69e-01

GCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTCCAGGCATTCATGGGCCTTGGC  >  NC_002947/4097564‑4097821                                                                                                               |                                                                                                                                                    gCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCGAGTGCATCAGCGCGACGACTGTCGAAGGt                                                                                                              >  1:335220/1‑150 (MQ=255)         ccAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGATAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATgg                                                                                                      >  2:432619/1‑150 (MQ=255)         ccAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGATAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATgg                                                                                                      <  1:432619/150‑1 (MQ=255)         ccAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATgg                                                                                                      >  2:389330/1‑150 (MQ=255)                       aaaGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCa                                                                                        >  1:193659/1‑150 (MQ=255)                                                 cGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCGATGAAATCAGGGGGCGGATCa                                                              >  2:94059/1‑150 (MQ=255)                                                 cGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCa                                                              <  2:193659/150‑1 (MQ=255)                                                   tAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAgg                                                            <  1:389330/150‑1 (MQ=255)                                                      ggTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAg                                                         <  2:159922/150‑1 (MQ=255)                                                            gcTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGAACAGGCAGTCAGCa                                                   >  1:91475/1‑150 (MQ=255)                                                             cTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCAAGCAGCGTGGCAAGCAAAAAAAAAAGACTGGCGTCGGACCTGTCCTGGGCAATGCTGGGGTTAGGGCAGTTAAATCAGGGGGCGGATCAGGCAGTCAGCAc                                                  <  2:384652/150‑1 (MQ=255)                                                                                         gAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTcc                      >  1:379101/1‑150 (MQ=255)                                                                                                      aCAAGGCCATCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTCCAGGCATTCATg           >  1:6012/1‑148 (MQ=255)                                                                                                      aCAAGGCCATCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTCCAGGCATTCATg           >  1:6006/1‑148 (MQ=255)                                                                                                      aCAAGGCCATCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTCCAGGCATTCATg           <  2:6006/148‑1 (MQ=255)                                                                                                      aCAAGGCCATCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTCCAGGCATTCATg           <  2:6012/148‑1 (MQ=255)                                                                                                       cAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTCCAGGCATTCATGGGc        <  2:329192/150‑1 (MQ=255)                                                                                                         aGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTCCAGGCATTCATGGGCCt      >  2:391815/1‑150 (MQ=255)                                                                                                             cagcagGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTCCAGGCATTCATGGGCCTTGgc  >  2:162151/1‑150 (MQ=255)                                                                                                               |                                                                                                                                                    GCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCTCCGACCGCCCTTGCTCCAGGCATTCATGGGCCTTGGC  >  NC_002947/4097564‑4097821

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