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

TGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAC  >  NC_002947/3845559‑3845837                                                                                                                                                  |                                                                                                                                       tgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑gggc                                                                                                                                    >  6:142356/1‑147 (MQ=255)          cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATc                                                                                                                            <  1:126405/149‑1 (MQ=255)              aTCATTTGCAGGTTCTGCATGCCTTCGCGAAGGCCTACTTCATTGATTTTAATTGTTGTCATCGGGTGCTCTCTGTCGTACGCGTTGCGCATGACTAATAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTc                                                                                                                         <  1:264978/148‑1 (MQ=255)                 aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTc                                                                                                                         <  2:299621/145‑1 (MQ=255)                 aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTc                                                                                                                         >  1:299621/1‑145 (MQ=255)                      cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                                                <  2:182412/149‑1 (MQ=255)                      cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                                                <  5:123493/149‑1 (MQ=255)                         gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCACCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                             <  1:8091/149‑1 (MQ=255)                             tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                 >  8:101051/1‑141 (MQ=255)                             tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                 <  7:101051/141‑1 (MQ=255)                                  gCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATcc                                                                                                                               <  6:89318/122‑1 (MQ=255)                                  gCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATcc                                                                                                                               >  5:89318/1‑122 (MQ=255)                                                   cTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTg                                                                                   <  5:142356/149‑1 (MQ=255)                                                          gattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgcag                                                                            >  2:22750/1‑149 (MQ=255)                                                           attgttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgcag                                                                            >  1:143615/1‑148 (MQ=255)                                                           attgttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgcag                                                                            >  3:286976/1‑148 (MQ=255)                                                             tgttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCa                                                                                 >  3:95360/1‑141 (MQ=255)                                                                      tgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCgg                                                                                     >  7:335756/1‑128 (MQ=255)                                                                             gcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAgg                                                         >  5:167021/1‑149 (MQ=255)                                                                              cgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGACCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTAAAGCAGgc                                                        >  1:129086/1‑149 (MQ=255)                                                                                     tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                 >  6:148969/1‑149 (MQ=255)                                                                                       tgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTGAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACATGTGGTCAAGCAGGCgccgagcc                                                >  3:127653/1‑148 (MQ=255)                                                                                             aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                          >  6:82193/1‑148 (MQ=255)                                                                                             aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCg                                         >  2:32342/1‑149 (MQ=255)                                                                                                  ttGCGCAACACTAAGAGCTGGATGACGGGCAGAAATGACGTGTTTCGCGGGCTTATCCTTCGGCCGCCAGTATGGATATGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGc                                    <  4:127653/149‑1 (MQ=255)                                                                                                             tatgagCTGGATGACGGGCTGAAAAGCCGTGTTTCGCGGGCTGATCCTTCGTCAGCCACTATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                         <  4:286976/146‑1 (MQ=255)                                                                                                               agagCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGTTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGTCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGCCAGCGAGAGCGTGTGCCGGTcgc                       >  7:309673/1‑149 (MQ=255)                                                                                                               agagCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACCATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCACGCGCCGAGCCGACAGCCAGAGCGTGTGCCGGTcgc                       >  7:34959/1‑149 (MQ=255)                                                                                                                gagCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGAGCCATGGTGAGCCACGATGGATCAGCCCGCGCCCAGGCAGGCCGGGGCAGGCCGCCACCGGGGGTCAAGCACGCGGCGGGGCGACCGGGGGAGGGGGGGCCGGTcggg                      >  8:306571/1‑147 (MQ=255)                                                                                                                 agCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGGCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCgcgc                     <  5:82193/149‑1 (MQ=255)                                                                                                                 agCTGGAGGCCGGGCAGAAACGACGTGTTCCGCGGGCTGATCCATCGGCAGCCACGAGGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCgcgc                     <  2:80650/149‑1 (MQ=255)                                                                                                                     ggATGACGGGCAGAAAAGACGTGTTTCGCGGGCAGATCCATCGTCAGCCACGATGCATCAGCCCGCCCCCTGCCAGGCCGGTGCACGCAGCAACACGTGGTCAAGCAGGCCCCGAGCCGCCAGCGAGAGCGTCTGCCGGTCGCGCATg                  >  8:14673/1‑148 (MQ=255)                                                                                                                      gATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtg                                                                   >  3:352832/1‑98 (MQ=255)                                                                                                                             ggCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCATGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc          >  8:279224/1‑148 (MQ=255)                                                                                                                              gCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc          <  6:339851/147‑1 (MQ=255)                                                                                                                              gCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc          >  5:339851/1‑147 (MQ=255)                                                                                                                               cAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGaa       >  2:123726/1‑149 (MQ=255)                                                                                                                                  aaaaGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGc                                    >  7:34875/1‑117 (MQ=255)                                                                                                                                    aaGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAc  <  1:22750/149‑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: ‑