Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	648,289	A→T	16.9%	intergenic (+58/+33)	acoR → / ← accC	acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	648,289	0	A	T	16.9%	69.0 / 7.8	35	intergenic (+58/+33)	acoR/accC	acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (12/15);  new base T (2/4);  total (16/19)
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

TGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTCGA  >  NC_002947/648144‑648437                                                                                                                                                   |                                                                                                                                                     tGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATTAGCACCGCATCTACATGCCGCTGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAgg                                                                                                                                                    <  7:35707/148‑1 (MQ=255)       aCGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGGCCCCGCTTCAGGCGCGGCAAAAAAAAAACCC‑GCACCAGGCCgggg                                                                                                                                             >  3:33592/1‑149 (MQ=255)       aCGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAGAACCCCGGCCCTGTGCCgggg                                                                                                                                             <  8:106872/149‑1 (MQ=255)             tttGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGTTGAAGGGGTATAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTgg                                                                                                                                       <  8:222296/149‑1 (MQ=255)               tGGGCATGTCGCGGGCGACGCTTTATCCGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAACACCCCGGCCCTGTGCCGGGGtttttta                                                                                                                                      <  5:5733/148‑6 (MQ=255)                gggCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCCCTGTGCCGGGGttttttat                                                                                                                                     <  7:242400/148‑7 (MQ=255)                 ggCATTTCGCGGGCGACGCTTTATCGGCGAGTGCAGGAGCACCGAATCTACATGCCGCTGATGAAGGTGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACa                                                                                                                                   <  4:33592/149‑1 (MQ=255)                           gggCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTaca                                                                                                                         <  4:342031/149‑1 (MQ=255)                                      tATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTcc                                                                                                               <  5:333240/148‑1 (MQ=255)                                                        cACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCCCCGCCCCACAATCTGGC‑AGGCAAAAAAAAAACCCCGCACACGGCCGGGGGTGTGGTAAATAAAACAAAACGATAACACAGGGGTGGCCACGTTGTGCt                                                                                            >  6:342245/1‑149 (MQ=255)                                                           cgaaTCGACATGCCGCGGATGATGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                         <  2:124021/146‑1 (MQ=255)                                                             cATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAgg                                                                                       <  1:359117/149‑1 (MQ=255)                                                             cATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAgg                                                                                       <  5:342245/149‑1 (MQ=255)                                                                 gACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAg                                                                                    <  8:151386/148‑1 (MQ=255)                                                                                aaGGGGTAGAGCCCTGGGGGCGCGCCGCGGCCCCAGATTATGGG‑CGGCGATAAAAAAACCCCGCACCGGGGCGGGGGTGGGGTTCAGAACACAACACCAGCACGCCGGGGTGGGCAGGTTGTGGTCGAGGGAGGGCATGTTGAcgccgc                                                                    >  8:103739/1‑149 (MQ=255)                                                                                             cTGGGGCCGCGTCGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGcc                                                       >  1:77902/1‑149 (MQ=255)                                                                                                               cccAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTc                                     >  6:309339/1‑149 (MQ=255)                                                                                                                ccAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                    >  5:19387/1‑133 (MQ=255)                                                                                                                ccAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                    <  6:19387/133‑1 (MQ=255)                                                                                                                    aTTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGt                                >  5:307131/1‑149 (MQ=255)                                                                                                                                   tAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTg                  <  5:309339/149‑1 (MQ=255)                                                                                                                                   tAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACCCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCggg                        <  1:217581/143‑1 (MQ=255)                                                                                                                                   tAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACCCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCggg                        >  2:217581/1‑143 (MQ=255)                                                                                                                                    ccagaaCCCCGGCCCTGTGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTgg                 >  4:266117/5‑149 (MQ=255)                                                                                                                                      aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGc                           >  2:349218/1‑137 (MQ=255)                                                                                                                                      aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGc                           <  1:349218/137‑1 (MQ=255)                                                                                                                                       gaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCt              >  7:11444/2‑149 (MQ=255)                                                                                                                                       aaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCt              <  7:103739/149‑1 (MQ=255)                                                                                                                                           cccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGAt          >  8:5158/1‑149 (MQ=255)                                                                                                                                            ccGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATg         <  8:11444/149‑1 (MQ=255)                                                                                                                                            ccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa           <  2:73009/147‑1 (MQ=255)                                                                                                                                            ccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa           >  1:73009/1‑147 (MQ=255)                                                                                                                                            ccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATg         >  8:52630/1‑149 (MQ=255)                                                                                                                                            ccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATg         >  7:69469/1‑149 (MQ=255)                                                                                                                                               gCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCg      >  4:112662/1‑149 (MQ=255)                                                                                                                                                cacaGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGt     <  7:5158/149‑1 (MQ=255)                                                                                                                                                   aGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCAcc                                          <  7:81669/109‑1 (MQ=255)                                                                                                                                                   aGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCAcc                                          >  8:81669/1‑109 (MQ=255)                                                                                                                                                   aGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTcga  <  7:52630/149‑1 (MQ=255)                                                                                                                                                   aGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGATGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGAt                      <  2:150810/129‑1 (MQ=255)                                                                                                                                                   aGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGATGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGAt                      >  1:150810/1‑129 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGC‑AGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTCGA  >  NC_002947/648144‑648437

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