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	88.9%	106.1 / 7.1	37	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/3);  new base G (17/15);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑C‑GG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  NC_002947/3845557‑3845842                                                                                                                                                      |                                                                                                                                                 gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑C‑gg                                                                                                                                                <  8:374936/149‑1 (MQ=255)      tgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGATGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CTGG‑CTGa                                                                                                                                           >  4:230196/1‑149 (MQ=255)          ggCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATcc                                                                                                                                        <  2:382556/148‑1 (MQ=255)          ggCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCa                                                                                                                                       >  5:185420/1‑149 (MQ=255)          ggCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑C‑GGTCTGATCCa                                                                                                                                       <  1:129577/149‑1 (MQ=255)            cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATc                                                                                                                                     >  5:222804/1‑149 (MQ=255)             aTGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCg                                                                                                                                    <  7:92619/149‑1 (MQ=255)               gATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGt                                                                                                                                   >  5:407189/1‑148 (MQ=255)                     ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCAGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCAc                                                                                                                            >  7:199036/1‑149 (MQ=255)                         aGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATg                                                                                                                        >  7:250495/1‑149 (MQ=255)                          ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATgg                                                                                                                       <  5:186105/149‑1 (MQ=255)                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGACGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGAt                                                                                                                     >  5:313111/1‑148 (MQ=255)                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCa                                                                                                                                 <  6:409263/136‑1 (MQ=255)                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCa                                                                                                                                 >  5:409263/1‑136 (MQ=255)                                 cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGc                                                                                                                 <  3:311903/148‑1 (MQ=255)                                                         attgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACg                                                                                    >  4:33360/1‑149 (MQ=255)                                                             attgttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACgcag                                                                                 >  7:173869/1‑148 (MQ=255)                                                              ttgttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACgcagc                                                                                >  7:318059/1‑148 (MQ=255)                                                                    ttgttgTCATCGCGGTCTCTCTGTCTTAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAACCCGTGTTTC‑G‑CGGG‑CTTATCCATCGTCAGCCACTATGTATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAgg                                                                          <  4:347875/148‑1 (MQ=255)                                                                      gttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCt                                                                                                    <  5:397059/120‑1 (MQ=255)                                                                      gttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCt                                                                                                    >  6:397059/1‑120 (MQ=255)                                                                         gtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCC‑GG‑CATAAGCATCGTGAGCCCAGATTGGACAAC‑‑‑‑CCGCCCCCAGGCAGGGCCGTGCACGCAGCAACAGGTGGTc                                                                     >  2:323826/1‑148 (MQ=255)                                                                                       tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                      >  3:73191/1‑149 (MQ=255)                                                                                           tctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagcc                                                     <  6:407189/146‑1 (MQ=255)                                                                                             tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAg                                                >  7:66301/1‑149 (MQ=255)                                                                                                  cGTTGCGCCATACTACGCGGTGGAGGAGGGGCTCAAACGAGGTGTTGC‑G‑CGGG‑CTCATCCATCGTCAGCCACGATGTATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgaga                                           <  6:313111/149‑1 (MQ=255)                                                                                                        gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt                                     >  6:165580/1‑149 (MQ=255)                                                                                                         ccAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGtgtg                                    >  7:306786/1‑149 (MQ=255)                                                                                                          cAGACTAAGCGATGCATGACGGGAACAAAAGACTTGTTTC‑G‑CGGG‑CTCATCCATCTCCTGCCACTATGTATCAGCCCCTCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgt                                       <  1:323826/149‑1 (MQ=255)                                                                                                          aaGACTAAGAGCTGTATGACGGGAAGAAAATACGTGTTTC‑G‑CGGG‑CTGATCCATCTGCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                                   <  8:199036/148‑1 (MQ=255)                                                                                                            gACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCg                                 >  1:28067/1‑149 (MQ=255)                                                                                                                 agagCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                                           <  7:272095/102‑1 (MQ=255)                                                                                                                 agagCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                                           >  8:272095/1‑102 (MQ=255)                                                                                                                      gggATGAGGGCCATAAACTCCGTGTTGC‑G‑CGGG‑ATTATCAATCGCCACCAACTTTGAAACAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                       <  3:230196/148‑1 (MQ=255)                                                                                                                            aCGGGCACAAAAGCCGTGTTTC‑G‑C‑GGCCGGAACCATCGTCACCCACGAGGTATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt                  <  7:86505/148‑1 (MQ=255)                                                                                                                             cGGGCATAAAAGAGGTTTTTC‑G‑CGGG‑CTTTTCCATCGTCAGCCACGTTGTATCTGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATg                 <  4:73191/148‑1 (MQ=255)                                                                                                                               ggCAGAAAAGACGTGTTTC‑G‑CGGG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACg              <  8:66301/149‑1 (MQ=255)                                                                                                                                           gtgtTTCGG‑G‑GG‑CTTATCCATCGTCAGCCACTTTGTATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc  <  2:28067/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑G‑C‑GG‑CTGATCCATCGTCAGCCACGATGGATCAGC‑‑‑‑CCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  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: ‑