Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	648,287	A→C	9.7%	intergenic (+56/+35)	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,287	0	A	C	9.7%	111.5 / 2.9	39	intergenic (+56/+35)	acoR/accC	acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (14/21);  new base C (2/2);  total (16/23)
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

TGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCC  >  NC_002947/648144‑648432                                                                                                                                                |                                                                                                                                                  tGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGGGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCCACGCGGCCCCCCATTATGGCAGGCAATAAAAAAACCCCGCCCAcgg                                                                                                                                              >  4:267391/1‑149 (MQ=255)    ggCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCg                                                                                                                                           <  7:243415/149‑1 (MQ=255)      cACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCggg                                                                                                                                         <  8:101939/149‑1 (MQ=255)          ggCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACTTGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTc                                                                                                                                     <  1:77321/149‑1 (MQ=255)           gCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCTGTGCCGGGGtttt                                                                                                                                    <  4:307341/149‑3 (MQ=255)            cTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTg                                                                                                                                   <  8:143002/149‑1 (MQ=255)               tGGGCATTTCGCGGGCGACGCTTTATCGCCGCGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGt                                                                                                                                 <  5:13915/148‑1 (MQ=255)                 ggCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTAc                                                                                                                               <  8:188833/148‑1 (MQ=255)                                            cGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTgg                                                                                                   <  1:229435/149‑1 (MQ=255)                                                   attaGCACCGCATCGACATGCCGCGGATTAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGttt                                                                                            <  5:71120/146‑1 (MQ=255)                                                        cACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCt                                                                                       <  3:53354/149‑1 (MQ=255)                                                           cGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                    <  3:267391/149‑1 (MQ=255)                                                           cGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                    <  1:24150/149‑1 (MQ=255)                                                                  aCATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTg                                                                             <  4:135709/149‑1 (MQ=255)                                                                        cgcgGATGAAGGGGTAGCGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGtt                                                                       <  7:261311/149‑1 (MQ=255)                                                                           ggATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAc                                                                    <  2:28523/149‑1 (MQ=255)                                                                                                ggcccGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt                                                <  5:221376/145‑1 (MQ=255)                                                                                                    cgcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtcatca                                           >  6:242209/1‑149 (MQ=255)                                                                                                     gcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCg                                                                                     >  4:136940/1‑106 (MQ=255)                                                                                                     gcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCg                                                                                     <  3:136940/106‑1 (MQ=255)                                                                                                              ccccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTg                                                          >  3:68849/1‑124 (MQ=255)                                                                                                              ccccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTg                                                          <  4:68849/124‑1 (MQ=255)                                                                                                                ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcg                               >  2:101487/1‑149 (MQ=255)                                                                                                                        tGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAg                       >  7:69266/1‑149 (MQ=255)                                                                                                                            aggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCgg                    >  4:203710/1‑148 (MQ=255)                                                                                                                                aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt               >  8:647/1‑149 (MQ=255)                                                                                                                                aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt               >  7:149261/1‑149 (MQ=255)                                                                                                                                aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt               >  2:910/1‑149 (MQ=255)                                                                                                                                  tAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTg             <  3:203710/149‑1 (MQ=255)                                                                                                                                    cagaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtcatc                                            >  5:45649/4‑116 (MQ=255)                                                                                                                                    cagaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtcatc                                            <  6:45649/113‑1 (MQ=255)                                                                                                                                    aaaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTGAATCACGCCCCAGGTCGCGGTGCAGAGGGATGTTGGt           >  1:110674/1‑149 (MQ=255)                                                                                                                                    aaaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGAGCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGt           >  1:114226/1‑149 (MQ=255)                                                                                                                                     aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTc          >  8:12920/1‑149 (MQ=255)                                                                                                                                         ccccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa      <  7:51260/149‑1 (MQ=255)                                                                                                                                          cccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCGGAAGCCTTCATCACGCACCAgg                                  <  4:294444/120‑1 (MQ=255)                                                                                                                                          cccGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCGGAAGCCTTCATCACGCACCAgg                                  >  3:294444/1‑120 (MQ=255)                                                                                                                                            cGGCACAGGGACGGGGTTCTGGTACAACACACTACACGATCACTCCTGCTTGGCCACTTTGAGCTCGAGGTAGTGAATGTTCACCCCGCATTTGCAGAAGCGTACATCACGCACCAGCTCGAGGTCCAGCGGGATGATGGTCTAGATGc   >  4:226250/1‑149 (MQ=255)                                                                                                                                             ggCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGcc  <  5:242209/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  TGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCC  >  NC_002947/648144‑648432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑