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%	132.9 / NA	40	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 (23/17);  total (23/17)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATT  >  NC_002947/3845557‑3845834                                                                                                                                                    |                                                                                                                                     gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑gg                                                                                                                                    >  4:199932/1‑149 (MQ=255)     ctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑gggctg                                                                                                                                <  4:290700/149‑5 (MQ=255)      tgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑gggctga                                                                                                                               >  8:139883/1‑144 (MQ=255)        tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGAt                                                                                                                               <  3:171630/148‑1 (MQ=255)           gCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCAt                                                                                                                           <  3:199932/149‑1 (MQ=255)           gCATGATCACTTGCAGGCTCTGCAAGCCATCGCACAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACAt                                                                                                                           >  6:44896/1‑149 (MQ=255)              tGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGAGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                                        <  2:111351/149‑1 (MQ=255)                      tGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCAc                                                                                                                 >  4:5042/1‑148 (MQ=255)                          ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATgg                                                                                                            <  3:5042/149‑1 (MQ=255)                                gCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGc                                                                                                     >  7:228917/1‑149 (MQ=255)                                 cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTGGCTGACGATGGA‑TCa                                                                                                       <  5:298149/146‑1 (MQ=255)                                    gCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATAGA‑TCAGCCcgc                                                                                                 <  1:182068/149‑1 (MQ=255)                                          gcgcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTg                                                                                           >  6:70971/1‑149 (MQ=255)                                           cgcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGc                                                                                          <  7:139883/149‑1 (MQ=255)                                                gCCGACTTCATTGATTGTTATTGTTGTAATCGCGGTCTCTCTGTCTGAAGAGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCATATAAATCGTAAGACACCATGGA‑TAAAACCGCGCCCGGCCAGGc                                                                                     >  5:37793/1‑149 (MQ=255)                                                                gttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                      >  7:271855/1‑148 (MQ=255)                                                                 ttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                                              >  2:230278/1‑108 (MQ=255)                                                                 ttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                                              <  1:230278/108‑1 (MQ=255)                                                                 ttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                      >  2:53435/1‑147 (MQ=255)                                                                  tattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                                   >  4:221770/1‑149 (MQ=255)                                                                             tCGCGTTCTCTCTGTCTTAAGGGTTGCGCCATAATAAGTGGTTGTTGAGGGGCATAAAAGACGTGTTTCGCGGGCTTATCCATCGTCAGCCACTTTGTA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCa                                                        <  8:271855/149‑1 (MQ=255)                                                                             tCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACATAGTCAGCCACGCTGGG‑TCACCCCCCGCCCGGCCAGGCCCGTGCACGCAGCAACACGTGGTCAAGCa                                                        >  1:154287/1‑149 (MQ=255)                                                                                cgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgc                                                     <  5:256875/149‑1 (MQ=255)                                                                                      ctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGc                                                                                                     <  4:41767/95‑1 (MQ=255)                                                                                      ctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGc                                                                                                     >  3:41767/1‑95 (MQ=255)                                                                                      ctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgag                                               >  8:225840/1‑149 (MQ=255)                                                                                          gtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagccga                                           >  5:223824/1‑149 (MQ=255)                                                                                           tctgAAGGGTTTCGCTAGACTAAGAGCTGTATGAGGGGCATAAACGACGTGTTTCGCGGGCTGTTCCATGGTCAGCCAGTTTGGT‑TCAGCCCGCGCCCTGCCAGGCGGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagccg                                            <  8:228917/147‑1 (MQ=255)                                                                                                    ttGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGc                                 <  6:223824/149‑1 (MQ=255)                                                                                                        gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt                             >  7:1608/1‑149 (MQ=255)                                                                                                         ccAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGCGCGCGAGAGCGtgtg                            >  1:106449/1‑149 (MQ=255)                                                                                                          cAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                           >  6:217127/1‑149 (MQ=255)                                                                                                           aGACTAAGAGCTGGATGAGGGGCAGAAACGACGGGTTGCGCGGGCTGATCCAGCGTCCGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGcc                          <  6:134193/149‑1 (MQ=255)                                                                                                           aGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGcc                          >  3:268042/1‑149 (MQ=255)                                                                                                               tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                      >  3:104938/1‑149 (MQ=255)                                                                                                                aaGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCAACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcg                     >  3:113245/1‑149 (MQ=255)                                                                                                                          tGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGCCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGa           <  1:53435/149‑1 (MQ=255)                                                                                                                             cGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGa        <  7:225840/149‑1 (MQ=255)                                                                                                                               ggCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGAGTCAGCCCGCGCCCTGCAAGGCCGGTGAACGCAGCAACAGGTGGTAAAGCAGCCCCCCAGCCGACAGCCACACCTTGTGCCGGTcg                     >  1:71263/1‑135 (MQ=255)                                                                                                                               ggCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc       >  1:293612/1‑148 (MQ=255)                                                                                                                                 cAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGAGCACGCAGCAACAGGTGGTCAAGCAGGAGCCGAGCCGTCAGCGAGAGCGTGTGCCGGTCGCGCAAGCAGATGACGaa    >  6:280995/1‑149 (MQ=255)                                                                                                                                              tttCGCGGGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAtt  >  4:49227/1‑138 (MQ=255)                                                                                                                                              tttCGCGGGCTGATCCAGCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAtt  <  3:49227/138‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                     GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGA‑TCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATT  >  NC_002947/3845557‑3845834

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: ‑