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	94.6%	115.6 / 0.6	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 . (0/2);  new base G (24/11);  total (24/13)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑GGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGC  >  NC_002947/3845564‑3845844                                                                                                                                              |                                                                                                                                               tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑gggctga                                                                                                                                         <  1:370700/147‑6 (MQ=255)     cATTATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTGATTGATTGTTATTGTTGTCATCGCGGGCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATc                                                                                                                                    <  5:115013/149‑1 (MQ=255)         aTCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGt                                                                                                                                  >  1:16239/1‑147 (MQ=255)          tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑g                                                                                                                              >  2:52240/1‑149 (MQ=255)                 cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGAt                                                                                                                       >  7:9585/1‑149 (MQ=255)                  aGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCAGATACAACGTCA‑GCCaaccag                                                                                                                      >  5:60571/1‑144 (MQ=255)                    gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGa                                                                                                                    >  1:90717/1‑149 (MQ=255)                          cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATAGTCA‑GCCACGATGGATCAGcc                                                                                                              >  2:167456/1‑149 (MQ=255)                                             aCTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGt                                                                                           >  7:140992/1‑149 (MQ=255)                                                  attgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑TCGGGCTGATTCAAGGTAACCCCACAATGG‑TAATCCCGCGCCCGGCCAGGCCgggggaag                                                                                      >  5:253208/1‑143 (MQ=255)                                                  attgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCAc                                                                                       <  1:52240/148‑1 (MQ=255)                                                   ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCTGAAAAGACGTGTTTC‑GCGGGCTGAt                                                                                                                                         >  3:312204/1‑98 (MQ=255)                                                   ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCTGAAAAGACGTGTTTC‑GCGGGCTGAt                                                                                                                                         <  4:312204/98‑1 (MQ=255)                                                         gttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                                >  7:96362/1‑148 (MQ=255)                                                          ttattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                                >  8:272611/1‑147 (MQ=255)                                                               gttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtgg                                                                         >  2:226640/1‑149 (MQ=255)                                                                      tCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCa                                                                  <  6:91939/149‑1 (MQ=255)                                                                         cgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATAAATGGGGA‑GGCACCATGGATAAACCCCCGCCCAGCCACGCCGGTGCACCCCGCAACACGGGGTCAAGCAGgc                                                               >  1:119675/1‑149 (MQ=255)                                                                                tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                        >  5:113173/1‑149 (MQ=255)                                                                                tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                        >  2:35061/1‑149 (MQ=255)                                                                                tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCA‑GCCGGGGGACGCAGCCACACGGGGTCAAGCAGgcgccgcgcc                                                       >  4:243923/1‑149 (MQ=255)                                                                                       gAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                                 >  3:329405/1‑149 (MQ=255)                                                                                        aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                                 >  1:254312/1‑148 (MQ=255)                                                                                                 gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgag                                              <  5:352801/142‑1 (MQ=255)                                                                                                 gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgag                                              >  6:352801/1‑142 (MQ=255)                                                                                                   cataCTATGAGCTGGATTAGGGGCATAAAAGACGTGTTTCGGG‑GGGTGATCCATCGTCA‑GCCACGATGTATAAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                                     <  8:9585/146‑1 (MQ=255)                                                                                                       cTAAGAGCGGGAGGCCGGGCAGAAACGCCGGGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGt                                 <  6:99966/149‑1 (MQ=255)                                                                                                        tAAGAGCTGTATGACGGGCATAAAAGACGTTTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                                <  2:90717/149‑1 (MQ=255)                                                                                                        tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                        >  5:355137/1‑125 (MQ=255)                                                                                                        tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                        <  6:355137/125‑1 (MQ=255)                                                                                                          agagCTGGTTGCCGGGCATAAACGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGTATCTGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgc                              <  4:329405/149‑1 (MQ=255)                                                                                                          agagCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGCGCCGACAGCGAGAGCGTGTGCCGGTcgc                              >  5:175007/1‑149 (MQ=255)                                                                                                          agagCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgc                              >  7:178685/1‑149 (MQ=255)                                                                                                              ctgtttGAGGGGCAGAAAAGACGTGTTTCGGG‑GGCTGATCCTTCTTCA‑GCCACTATGGTTCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCa                           <  8:178685/143‑1 (MQ=255)                                                                                                                   tGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                                     >  3:347769/1‑133 (MQ=255)                                                                                                                                  aCGTGTTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCaa      <  1:35061/149‑1 (MQ=255)                                                                                                                                    gtgtTTC‑GCGGGCTGATCCATCGTCA‑GCCACGTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCATCAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc    <  7:272611/149‑1 (MQ=255)                                                                                                                                      gtTTC‑GCGGGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgcgc  >  2:239582/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                               TCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑GGCTGATCCATCGTCA‑GCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGC  >  NC_002947/3845564‑3845844

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