Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,845,705	+G	intergenic (‑72/+26)	PP_3394 ← / ← PP_3395	3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,845,703	1	.	G	100.0%	57.8 / NA	18	intergenic (‑70/+28)	PP_3394/PP_3395	3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (12/6);  total (12/6)

GATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCAT  >  NC_002947/3845571‑3845820                                                                                                                                      |                                                                                                                      gATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                         <  5:185392/148‑1 (MQ=255)            ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                               >  7:134988/1‑147 (MQ=255)                       ccATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGGCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCcgc                                                                                   <  1:45493/148‑1 (MQ=255)                                   ccGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACTATGGATCAGCCCGCGCCCTGCCAGGc                                                                       <  4:153721/148‑1 (MQ=255)                                       cTTAATTTATTGTTATTGTTGTCATCCGGGTCTCTCTGAGTCACCGGTTGCGCCTAACTAATTGCTGGATTCCGGGCACAAAATAGGTGTTTCGCGGGCTCACCCAGCGACAGCCACTTTGTATCAGCCCGCGCCCTGCCAGGCCGGt                                                                   <  6:86845/148‑1 (MQ=255)                                                   ttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACATCGTCAGCCCCGCTGGCTAAGCCCGCGCCCTGCCCGGCCGGTGCACgccgca                                                         >  5:117004/1‑146 (MQ=255)                                                      ttgttgTCAGCGGGGCCTCTCTGGGTGAAGCGTTGCGCCAGACTACGAGCTGGATGACGGGCAGAAAAGACGGTTTTCGCGGGCTCACCCAGCGTCCCCCACGTTGCATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                     <  2:128585/147‑1 (MQ=255)                                                           gtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtggt                                                >  1:37709/1‑147 (MQ=255)                                                              aTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCCACAGGTGGTCaa                                             >  2:74016/1‑147 (MQ=255)                                                                       tctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgcgc                                     >  5:135922/1‑146 (MQ=255)                                                                           tgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                >  3:76063/1‑147 (MQ=255)                                                                                  aGCGTTGCGCCAGACTAAGAGCTGGAGGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                         <  8:151812/147‑1 (MQ=255)                                                                                     gTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGCGCCGACAGCgaga                     >  2:135008/1‑148 (MQ=255)                                                                                          gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtg                >  3:112324/1‑148 (MQ=255)                                                                                           ccAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGCGCCGACAGCGAGAGCgtgt               >  1:7590/1‑148 (MQ=255)                                                                                            cAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCACCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGtgtg              >  1:9040/1‑148 (MQ=255)                                                                                                        tGGATGACGGGCAGAAAAGCCGTGTTTCGCGGGCTTATACATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGGCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCAt  >  4:36162/1‑148 (MQ=255)                                                                                                        tGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCAACGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGGGGTCACGCAGGCGCCGAGCCGACCGCGAGCGCGAGTGCCGGTCGCGCAt  >  4:33596/1‑148 (MQ=255)                                                                                                                                      |                                                                                                                      GATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCAT  >  NC_002947/3845571‑3845820

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