Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,845,703:1	+G	100%	intergenic (‑70/+28)	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%	106.3 / NA	32	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 (17/15);  total (17/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCA  >  NC_002947/3845559‑3845839                                                                                                                                                   |                                                                                                                                         tgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTcgc‑g                                                                                                                                         <  6:178714/146‑1 (MQ=255) tgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTcgc‑g                                                                                                                                         >  5:178714/1‑146 (MQ=255) tgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGC‑gggc                                                                                                                                      >  2:308756/1‑147 (MQ=255)   ctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGC‑gggc                                                                                                                                      <  6:305161/147‑3 (MQ=255)          cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGC‑gggc                                                                                                                                      <  8:348931/140‑3 (MQ=255)          cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGC‑gggc                                                                                                                                      >  7:348931/1‑138 (MQ=255)            tGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                                            >  2:106998/1‑149 (MQ=255)              aTCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTc                                                                                                                           <  8:69197/148‑1 (MQ=255)               tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAg                                                                                                                         <  7:32036/149‑1 (MQ=255)                cACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGc                                                                                                                        >  2:386695/1‑149 (MQ=255)                cACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGc                                                                                                                        >  3:60974/1‑149 (MQ=255)                   ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGtggctgac                                                                                                                     >  2:215392/1‑142 (MQ=255)                      cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                   <  8:411873/148‑1 (MQ=255)                         gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCAc                                                                                                                     <  1:231095/143‑1 (MQ=255)                         gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCAc                                                                                                                     >  2:231095/1‑143 (MQ=255)                           tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGAt                                                                                                              >  7:241774/1‑148 (MQ=255)                              gCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGc                                                                                                          >  5:197001/1‑149 (MQ=255)                               cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                                                  <  3:303119/140‑1 (MQ=255)                               cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                                                  >  4:303119/1‑140 (MQ=255)                               cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCa                                                                                                            <  7:73036/146‑1 (MQ=255)                               cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGc                                                                                                          <  8:200881/148‑1 (MQ=255)                                                cGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCg                                                                                        <  5:422114/149‑1 (MQ=255)                                                cGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCg                                                                                        >  7:18123/1‑149 (MQ=255)                                                                 attgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAgg                                                                       >  5:25928/1‑149 (MQ=255)                                                                              cgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAgg                                                           >  1:90861/1‑149 (MQ=255)                                                                                       tgtctgAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                   <  7:391469/147‑1 (MQ=255)                                                                                           tgAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCAAGCAGCAACAGGTGGTAAAGCAGGCGCCGAGCCGACAg                                             >  3:18141/1‑149 (MQ=255)                                                                                            gAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                            >  2:54463/1‑149 (MQ=255)                                                                                              aGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCga                                          <  1:106998/149‑1 (MQ=255)                                                                                                        cAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                                >  7:62563/1‑149 (MQ=255)                                                                                                             tAAGAGCTGTAT‑GCGGGGCACAAACGCCGTGTTGCGCGGGCTTACCCAGCGTCCCCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                           <  1:386695/149‑1 (MQ=255)                                                                                                                gagCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgcg                        <  6:197001/149‑1 (MQ=255)                                                                                                                  gCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCgcgc                       <  2:90861/148‑1 (MQ=255)                                                                                                                           aCGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt               <  1:54463/148‑1 (MQ=255)                                                                                                                            cGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGa             >  3:160035/1‑149 (MQ=255)                                                                                                                                 aGAAAAGACGTGTTTCGCGGGCTGATCAATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAt        >  6:332908/1‑149 (MQ=255)                                                                                                                                      agaggggtTTCGCGGGCTGATCCATCGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAcc   <  8:18123/143‑1 (MQ=255)                                                                                                                                        aCGTGTTTCGCGGGCTGATACAACGGCAGGCACGATGGATAAGCCCGCGCCCAGCCAGGCCCGTGCACGCAGCAACAGGTGGTAAAGCCCGCGCCGAGCCGCCAGCGGCAGAGTGTGCCGGTCGCGCGTGCAGATGACGGAttaccca  >  4:8215/1‑143 (MQ=255)                                                                                                                                                   |                                                                                                                                         TGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGAT‑GACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCA  >  NC_002947/3845559‑3845839

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