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	95.9%	157.3 / ‑0.3	49	intergenic (‑70/+28)	PP_3394/PP_3395	3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (24/23);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑C‑GG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAT  >  NC_002947/3845557‑3845833                                                                                                                                                     |                                                                                                                                    gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTcg‑c‑g                                                                                                                                    <  1:154572/148‑1 (MQ=255) gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTcg‑c‑g                                                                                                                                    <  3:92316/148‑1 (MQ=255)  ctgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGTAGAAAAGACGTGTTTCG‑C‑ggg                                                                                                                                  <  6:18102/149‑2 (MQ=255)           gCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCAt                                                                                                                         >  1:2891/1‑149 (MQ=255)            cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATc                                                                                                                        <  6:92953/149‑1 (MQ=255)             aTGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCg                                                                                                                       <  1:345534/149‑1 (MQ=255)                aTCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCa                                                                                                                    <  3:233555/149‑1 (MQ=255)                 tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAg                                                                                                                   >  4:97443/1‑149 (MQ=255)                   aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATc                                                                                                                            >  7:4330/1‑138 (MQ=255)                   aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATc                                                                                                                            <  8:4330/138‑1 (MQ=255)                     ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGAG‑GG‑CTGATCCATCGTCAGCCAc                                                                                                               >  3:75328/1‑149 (MQ=255)                            ctctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGAt                                                                                                        <  8:75723/149‑1 (MQ=255)                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGtggctga                                                                                                                <  6:123026/140‑7 (MQ=255)                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGtggctga                                                                                                                >  5:123026/1‑134 (MQ=255)                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGAt                                                                                                        >  4:293667/1‑148 (MQ=255)                                 cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCa                                                                                                      <  3:97443/146‑1 (MQ=255)                                   aGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGTGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATAAGCCcg                                                                                                 >  5:279569/1‑149 (MQ=255)                                    gCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCcgc                                                                                                >  6:111541/1‑149 (MQ=255)                                       aTCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGcc                                                                                             >  4:352411/1‑149 (MQ=255)                                            gcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGc                                                                                         >  7:94390/1‑148 (MQ=255)                                                  cGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCg                                                                                  >  3:318438/1‑149 (MQ=255)                                                   gACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCgg                                                                                 >  8:102572/1‑149 (MQ=255)                                                    aCTTCATTTTTTGTTTTTGTTGTCATCGCGGTCTCTCTGTCTTAAGCGTTGCGCCATACTAAGAGCTGGATGACGGGCAGAAAAGCCGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGt                                                                                <  8:94390/149‑1 (MQ=255)                                                     cTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTg                                                                               >  2:108524/1‑149 (MQ=255)                                                          ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgc                                                                          <  3:293667/149‑1 (MQ=255)                                                                 ttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                     >  7:108830/1‑147 (MQ=255)                                                                  tattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                                  >  1:256674/1‑149 (MQ=255)                                                                    ttgttgTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggt                                                                <  1:297422/149‑1 (MQ=255)                                                                      gttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGTAGAAAAGACGTGTTTCG‑CGGG‑CTGATACATCGTCAGCCACGATGGATCAGCCCCCGCCCTGCCCGGCCGGTGCACGCAGCAACAggtgg                                                              >  1:186512/1‑149 (MQ=255)                                                                        tgtCATCGCGGTCTCTCTGTCTCATGCGTTCCGCCAGACTAAGAGCTGGATGACGGGCATAAAACACTTGTTTCG‑CGGG‑CTTATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtggt                                                             <  2:256674/148‑1 (MQ=255)                                                                                gggTGTCTCTGGCTGAAGCGTTGCGCCTTATTAATTGGTGGATGACTGGCATAAATGACGTGTTTCG‑CGGG‑CTGATCCTTGGTCAGCCTCTTTTGTTCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgc                                                    <  6:279569/148‑1 (MQ=255)                                                                                       tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGc                                                                                                    <  8:223788/94‑1 (MQ=255)                                                                                       tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGc                                                                                                    >  7:223788/1‑94 (MQ=255)                                                                                         tgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACTATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                             <  4:75328/147‑1 (MQ=255)                                                                                         tgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGCTGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagac                                            >  5:68964/1‑146 (MQ=255)                                                                                           tctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagccga                                          >  1:204372/1‑148 (MQ=255)                                                                                             tgtaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAg                                       >  7:76833/4‑149 (MQ=255)                                                                                             tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgcg                                                   >  5:395735/1‑137 (MQ=255)                                                                                             tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgcg                                                   <  6:395735/137‑1 (MQ=255)                                                                                               aaGCGTTGCGCCAGACTAATCGCTGGAGGCCGGGCACAAAAGACGTGTTGCG‑CGGG‑CTGAGCCATCGTCCGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACa                                        <  7:104545/146‑1 (MQ=255)                                                                                                  gtttGCGCCATACTAAGAGCTGGTTGACGGGCATAAAAGACGTGTTTCG‑CGGG‑CTGATCCTTCGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgaga                                  <  8:108830/147‑1 (MQ=255)                                                                                                      gccccAGAATAAGCGCTGGATGAGGGGCAGAAAACCCGTGTTGCG‑CGGG‑CTCATCCAGCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgt                              <  3:352411/146‑1 (MQ=255)                                                                                                        gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACa                                        >  5:13034/1‑137 (MQ=255)                                                                                                           aGACTAAGAGCTGGATCCCGGGCACAAAAGCCGTGTTTCG‑CGGG‑CGGACCCAGCGCCAGCCACGAGGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGcc                         <  5:111541/149‑1 (MQ=255)                                                                                                                 tgagCTGGATTACGGTCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGAGGGATCAGCCGGCGCCCTGCCAGGCCGGTGCACTCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcg                    <  6:68964/147‑1 (MQ=255)                                                                                                                 agagCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgc                   <  7:102572/149‑1 (MQ=255)                                                                                                                  gagCTGGATGACGGGCACAAAACCCGTGTTTCG‑C‑GGCCTGATCCATCGTCAGCCACGATGCATAAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgcg                  <  4:29046/149‑1 (MQ=255)                                                                                                                   atcTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCgcgc                 <  3:247230/147‑1 (MQ=255)                                                                                                                    gCTGGATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCgcgc                 <  4:318438/148‑1 (MQ=255)                                                                                                                        gATGACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCgcgc                 >  4:247230/1‑144 (MQ=255)                                                                                                                           gACGGGCAGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt         >  6:89497/1‑149 (MQ=255)                                                                                                                                  aGAAAAGACGTGTTTCG‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAt  >  7:51221/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                    GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCG‑C‑GG‑CTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAT  >  NC_002947/3845557‑3845833

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