Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	648,289	A→T	10.6%	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	10.6%	93.9 / 3.7	38	intergenic (+58/+33)	acoR/accC	acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (15/19);  new base T (2/2);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

CTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTCGA  >  NC_002947/648155‑648437                                                                                                                                       |                                                                                                                                                      cTTTGGGCATTTCGCTGCCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTg                                                                                                                                         <  1:120407/149‑1 (MQ=255)    tGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGCGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAGAACCCCGGCCCTGTGCCGGGGtttttt                                                                                                                                        <  3:137267/147‑5 (MQ=255)     gggCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAG‑CCTGGGGCCGCGACCCGGCCCCACAAAATGGCAGGCAATAAAAAAACCCCCCACAAGGCCGGGGTTCTGGttc                                                                                                                                     >  2:306926/1‑146 (MQ=255)            tCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATcaca                                                                                                                               <  3:41167/148‑1 (MQ=255)                 ggCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGGCCCAGCTTCTGGCAGGGCATAAAAAAACCCCGCACAGGGGCCGGGTTGTGGTACAGAacacaaaac                                                                                                                         >  6:96124/1‑144 (MQ=255)                    caCGTTTTCTCGGCGTGTGCATTAGCACCGCATCGACATGCCGCGGATCACTGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGAt                                                                                                                     <  5:45715/148‑1 (MQ=255)                        cTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGGCCCAGATTGGGGGAGGGAATAAAAAAACCCCGCCCAGGGCCGGGGGTCGGGTTCATAACACTAACACGaccct                                                                                                                  >  6:169106/1‑145 (MQ=255)                                       cATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGtt                                                                                                  >  5:353116/1‑149 (MQ=255)                                        aTGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGttt                                                                                                 <  6:353116/149‑1 (MQ=255)                                                cGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                         <  5:96124/149‑1 (MQ=255)                                                             cgcgGATGAAGGGGTAGAG‑CCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTg                                                                           <  1:306926/149‑1 (MQ=255)                                                                ggATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAc                                                                         <  8:238455/149‑1 (MQ=255)                                                                 gATGAAGGGGTAGAGCCCTGGGGCCGCGACCCGGCCCCCGAATATGGCAGGCAATAAAAAAAACCCGCCCAGGGGCCGGGGTGTGGTACAACACACA‑AAACCATCACTCCTGGGTGGCCAGTTTTTGGTCGAGGTAGTGTATGtttaag                                                                        >  2:68149/1‑145 (MQ=255)                                                                      aGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGt                                                                                                           <  1:305444/109‑1 (MQ=255)                                                                             gagCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTtgt                                                                                               <  1:118129/114‑1 (MQ=255)                                                                                 tcTGGGGCCGCGTCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTa                                                                                     >  8:89789/2‑120 (MQ=255)                                                                                      ggCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAc                                                                         >  2:161058/1‑127 (MQ=255)                                                                                          gcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                    >  5:189628/1‑144 (MQ=255)                                                                                          gcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtc                                                    <  6:189628/144‑1 (MQ=255)                                                                                              cgcgGCCCCAGATGCTGGCAGGCTATAAAAAAGCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCAc                                           <  4:223676/149‑1 (MQ=255)                                                                                               gcgGCCCCAGATTCTTGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtca                                                   <  2:215273/140‑1 (MQ=255)                                                                                               gcgGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGa                                                                                 <  6:308905/110‑1 (MQ=255)                                                                                               gcgGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGa                                                                                 >  5:308905/1‑110 (MQ=255)                                                                                               gcgGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtca                                                   >  1:215273/1‑140 (MQ=255)                                                                                                 ggCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGc                                             >  8:11809/1‑144 (MQ=255)                                                                                                 ggCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGc                                             <  7:11809/144‑1 (MQ=255)                                                                                                     ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcg                                    >  2:153402/1‑149 (MQ=255)                                                                                                          ttCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTg                               >  5:266597/1‑149 (MQ=255)                                                                                                             tGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAg                            >  5:272488/1‑149 (MQ=255)                                                                                                                     aaTAACAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                    <  7:186924/149‑1 (MQ=255)                                                                                                                     aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                    >  6:326336/1‑149 (MQ=255)                                                                                                                       tAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAgg                                                                                       <  2:286091/80‑1 (MQ=255)                                                                                                                          agaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcgccgc                                                                    <  2:86347/94‑1 (MQ=38)                                                                                                                          agaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcgccgc                                                                    >  1:86347/3‑96 (MQ=38)                                                                                                                          agaaCCCCGGCCCTGTGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTc               >  3:142209/3‑149 (MQ=255)                                                                                                                              ccccGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGATGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa           <  7:328701/149‑1 (MQ=255)                                                                                                                               cccGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGAt          <  6:272488/149‑1 (MQ=255)                                                                                                                                  ggCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGGCTTGATGcc       <  6:266597/149‑1 (MQ=255)                                                                                                                                     acaGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTc    <  1:68149/149‑1 (MQ=255)                                                                                                                                       aGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                                                          <  5:232346/93‑1 (MQ=255)                                                                                                                                       aGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                                                          >  6:232346/1‑93 (MQ=255)                                                                                                                                       aGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTcga  >  7:297981/1‑149 (MQ=255)                                                                                                                                       |                                                                                                                                                      CTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACT‑ACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTCGA  >  NC_002947/648155‑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: ‑