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%	83.5 / NA	25	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 (18/7);  total (18/7)

TGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAC  >  NC_002947/3845559‑3845837                                                                                                                                                  |                                                                                                                                       tgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑ggg                                                                                                                                     <  4:132669/148‑2 (MQ=255)           aTGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATc                                                                                                                            <  1:121347/148‑1 (MQ=255)              atggcTTGCAGAGCCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCg                                                                                                                           >  2:136402/5‑146 (MQ=255)                   ttGCACGCTCTGCAAGCCATCGCGCAGGCCGCCTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCCTCGTCAGCCa                                                                                                                    >  1:142827/1‑148 (MQ=255)                        ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATg                                                                                                               <  5:69509/148‑1 (MQ=255)                           tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGGCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGGTTCGCGGGCTGATCCATCGTCAGCCACGCTGGa                                                                                                             >  6:197293/1‑147 (MQ=255)                                aaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTGGCTGACGATGGATCa                                                                                                          <  7:56727/145‑1 (MQ=255)                                  gCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCcg                                                                                                     >  8:115451/1‑148 (MQ=255)                                               ccGTCTACATTTATTGTTATTGTTGTCTTCGCGGTCTGCCGGTCTGAACCTTTGCGCCAGACTAACTGCTGGATGAGGGGCCTAAAAGCCGTGTTTCGCGGGCTGATCCATCGCCAGCCACGATGGATCAGCCCGCGCCCTGCCAg                                                                                          <  6:54824/146‑1 (MQ=255)                                                    ttCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGt                                                                                    >  6:176841/1‑147 (MQ=255)                                                       attgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCAc                                                                                >  8:6541/1‑148 (MQ=255)                                                              gttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACCTCGTCAGCCACCATGGATAAGCCCCCGCCCTGCCAGGCCCGTGCAcgccgcc                                                                          >  8:69860/1‑146 (MQ=255)                                                                 attgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCCGGCCGGTGCACGCAGCAACAg                                                                      >  7:60044/1‑148 (MQ=255)                                                                          aTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCaa                                                              >  7:38426/1‑147 (MQ=255)                                                                          aTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCaa                                                              >  7:149094/1‑147 (MQ=255)                                                                                   tctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccg                                                    >  8:91530/1‑148 (MQ=255)                                                                                     tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCCCGCAGCAACAGGTGGTCAAGCAGGCgccgag                                                  >  5:141075/1‑148 (MQ=255)                                                                                                      gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtg                                 >  7:157077/1‑148 (MQ=255)                                                                                                        cAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGCGCGtgtg                               >  3:106979/1‑148 (MQ=255)                                                                                                         aGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                              >  2:52647/1‑148 (MQ=255)                                                                                                           aCTAAGAGCGGGATCCCGGGCACAAACCCCGGGTTGCGCGGGCTCACCCAGCGCCCCCCACTAGGGATCCGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCg                            <  7:115451/148‑1 (MQ=255)                                                                                                             tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGt                          >  2:33470/1‑148 (MQ=255)                                                                                                                        tGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATAAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACACCGAGAGCCTGTGCCGGGCGCGCAAGCa                >  7:129090/1‑147 (MQ=255)                                                                                                                          aCGGGCATAAATGAGGTGTTTCGCGGGGTGACCCTTCGTCAGCAAGTTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGa              <  6:141075/147‑1 (MQ=255)                                                                                                                            gggCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCCGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGCGCCGCCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGa           >  3:160042/1‑148 (MQ=255)                                                                                                                                     aGACGTGTTTCGCGGGCTTATCAATCGTCAGCCACTATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAc  <  8:149094/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                       TGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAC  >  NC_002947/3845559‑3845837

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