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	97.6%	139.0 / ‑3.4	42	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/0);  new base G (24/17);  total (25/17)
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

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑G‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  NC_002947/3845557‑3845842                                                                                                                                                    |                                                                                                                                             gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑g‑g                                                                                                                                           <  7:135625/149‑1 (MQ=255)  ctgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑g‑gg                                                                                                                                          >  5:41066/1‑148 (MQ=255)        tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATc                                                                                                                                     <  2:97055/149‑1 (MQ=255)               gATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGt                                                                                                                               >  5:117468/1‑148 (MQ=255)                     ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTAATTGATTGTTATTGTTGTCATCGCCGTCTCTCTGTCTGAAGCGTTGCGCCAGAATAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GCGCTGATCCATCGTCAGCCAc                                                                                                                        >  2:45830/1‑149 (MQ=255)                           gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCa                                                                                                                             >  7:184078/1‑138 (MQ=255)                           gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGa                                                                                                                  <  4:207139/149‑1 (MQ=255)                            ctctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGAt                                                                                                                 <  6:41066/149‑1 (MQ=255)                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATc                                                                                                                                 <  4:343045/132‑1 (MQ=255)                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATc                                                                                                                                 >  3:343045/1‑132 (MQ=255)                                   aGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATg                                                                                                                    >  1:44763/1‑139 (MQ=255)                                   aGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATg                                                                                                                    <  2:44763/139‑1 (MQ=255)                                         cgcgCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGAATAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCt                                                                                                    >  7:269495/1‑149 (MQ=255)                                                      ttCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGc                                                                                                  >  6:91067/1‑138 (MQ=255)                                                      ttCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGc                                                                                                  <  5:91067/138‑1 (MQ=255)                                                        cATTGATTGTTTTTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGGGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGGGTTTCGCGG‑GCTGATCCTTCGTCAGCCACGATGGAGCAGCCCGCGCCCTGCCAGGCCGGTGc                                                                                       <  1:45830/147‑1 (MQ=255)                                                                 ttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATgg                                                                                                                   <  7:283427/110‑1 (MQ=255)                                                                   attgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGGCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                                           <  6:167649/148‑1 (MQ=255)                                                                    ttgttgTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGccc                                                                                                           <  5:357841/115‑1 (MQ=255)                                                                    ttgttgTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGccc                                                                                                           >  6:357841/1‑115 (MQ=255)                                                                    ttgttgTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGAAGGATCAGACCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggt                                                                         >  3:258356/1‑149 (MQ=255)                                                                            aTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAg                                                                  >  6:161382/1‑148 (MQ=255)                                                                             tCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGa                                                                                                                      >  7:348648/1‑95 (MQ=255)                                                                              cgcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAg                                                               >  2:279269/1‑149 (MQ=255)                                                                                      ctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgag                                                       >  6:365391/1‑149 (MQ=255)                                                                                                  cGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCga                                             >  2:64178/1‑147 (MQ=255)                                                                                                    ttGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGc                                         >  5:121308/1‑149 (MQ=255)                                                                                                       cgcAAGACTAAGAGCTGGATGACGGGCAGAAATGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgt                                       <  6:87702/148‑1 (MQ=255)                                                                                                                 agagCTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgc                            <  1:279269/149‑1 (MQ=255)                                                                                                                     cTGGATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCCTCAGCCACCAAGGATCACCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACACGTGGTCAAGCAAGCCCCGAGCCGCCAGCGAGACGGTGTGCCGGAcgccc                          >  3:201145/1‑145 (MQ=255)                                                                                                                       ggATGACGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                                           >  1:337585/1‑96 (MQ=255)                                                                                                                        gctggcGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCGACGTTTGATCCTCCCGCGCCCTGCCAGGCCGGGGGGCGCAGCAACAGGTGGTCCAGCAGGCGCCGAGCCGACAGCGAGAGCGGGTGCCGGTCGCGGATg                       >  7:72997/6‑147 (MQ=255)                                                                                                                            aCGGGCAGAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt                  <  5:161382/148‑1 (MQ=255)                                                                                                                                   gAAAAGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGCCAGCGAGAGCGTGGGCCGGTCGCGCATGCAGATGACGCAtt          >  5:128848/1‑149 (MQ=255)                                                                                                                                    aaaaGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAtt          >  2:56721/1‑148 (MQ=255)                                                                                                                                    aaaaGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTc         >  5:286686/1‑149 (MQ=255)                                                                                                                                     aaaGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCa        <  6:121308/149‑1 (MQ=255)                                                                                                                                      aaGACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAc       <  6:117468/149‑1 (MQ=255)                                                                                                                                        gACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAcc      <  6:340609/148‑1 (MQ=255)                                                                                                                                        gACGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAcc      >  5:340609/1‑148 (MQ=255)                                                                                                                                         cggtgtTTGGCGG‑GCTGACCAAGCGCCCGCCCCGATGCATCCCCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCaa    <  7:199325/147‑1 (MQ=255)                                                                                                                                         aCGTGTTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt                                     >  1:135982/1‑116 (MQ=255)                                                                                                                                            tgtTTCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCga                                             >  7:358974/1‑105 (MQ=255)                                                                                                                                               ttCGCGG‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGCGCCGGCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc  >  1:304838/1‑145 (MQ=255)                                                                                                                                                    |                                                                                                                                             GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑G‑GCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  NC_002947/3845557‑3845842

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