Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,593,973	C→T	8.7%	intergenic (+279/+20)	mccA → / ← PP_4068	methylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha/Cro/CI family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,593,973	0	C	T	8.7%	127.6 / 2.7	46	intergenic (+279/+20)	mccA/PP_4068	methylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha/Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (25/17);  new base T (2/2);  total (27/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 = 8.62e-01

CGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTCCTCGGCCAGCTCGCG  >  NC_002947/4593832‑4594118                                                                                                                                              |                                                                                                                                                  cGTATTCGCGGGT‑AACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCacta                                                                                                                                           >  4:44012/1‑149 (MQ=255)          gggTAAACCCGCTCGTACATCCGATCGCGCGGGTACACATGTGGGTGACTAGCCTGACGCCGCGGGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAAcc                                                                                                                                   <  7:304658/149‑1 (MQ=18)                aCCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTc                                                                                                                              <  3:44012/148‑1 (MQ=255)                         tACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCac                                                                                                                                             <  6:23468/124‑1 (MQ=255)                         tACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCac                                                                                                                                             >  5:23468/1‑124 (MQ=255)                            aTCCGAACGGGGGTGGAAACATGTTGCGGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTGCACGCGAAAGGGCCGGCCCTACTCGCCCCCTTCCACCCGCCGTTCACCGTTGCCGCCACTAGTCCAGCGCTTTCATAAAGCCCCATg                                                                                                                 >  6:73005/1‑149 (MQ=21)                                 aTCGGGGGTGGACACATGATGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCAccc                                                                                                            >  2:349441/1‑149 (MQ=255)                                           gACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTTCTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGaccac                                                                                                  >  5:330184/1‑149 (MQ=255)                                                aTGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCAcccc                                                                                               >  1:107930/1‑147 (MQ=255)                                                 tGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCAcccc                                                                                               >  2:43848/1‑146 (MQ=255)                                                      gTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGTCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAAtt                                                                                       >  8:78755/1‑149 (MQ=255)                                                           tAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCAGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCaa                                                                                  >  3:66633/1‑149 (MQ=255)                                                                      cgcgTGTAGGAGCGGGTTTACTCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGTCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGca                                                                       >  4:290443/1‑149 (MQ=255)                                                                      cgcgTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGca                                                                       >  6:249431/1‑149 (MQ=255)                                                                       gcgTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCATGACCACCCCAAGAAATTCGCAACCCTCTTCGca                                                                       >  3:332638/1‑148 (MQ=255)                                                                                aGCGGGTTTACTCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGTCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGcc                                                             <  3:290443/149‑1 (MQ=255)                                                                                aGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGcc                                                             <  6:330184/149‑1 (MQ=255)                                                                                 gCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCATGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCg                                                            <  4:332638/149‑1 (MQ=255)                                                                                 gCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCg                                                            >  3:203357/1‑149 (MQ=255)                                                                                 gCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCg                                                            >  7:17301/1‑149 (MQ=255)                                                                                   gggTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCggg                                                          >  7:41426/1‑149 (MQ=255)                                                                                     gTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCtt                                                                                                                                >  6:169177/1‑77 (MQ=255)                                                                                     gTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCtt                                                                                                                                <  5:169177/77‑1 (MQ=255)                                                                                     gTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTa                                                        <  5:249431/149‑1 (MQ=255)                                                                                     gTTTACCCACGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCggg                                                          >  7:156738/1‑147 (MQ=255)                                                                                            cgcgAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTg                                                 <  1:349441/149‑1 (MQ=255)                                                                                             gcgAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTACACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGt                                                >  4:75786/1‑149 (MQ=255)                                                                                               gAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGTCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTc                                              <  7:78755/149‑1 (MQ=255)                                                                                               gAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGtt                                               >  8:40734/1‑148 (MQ=255)                                                                                                 aaGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCaa                                            <  1:43848/149‑1 (MQ=255)                                                                                                        ggccCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTc                                     <  3:75786/149‑1 (MQ=255)                                                                                                            ctgctgGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGa                                  >  6:92708/1‑148 (MQ=255)                                                                                                            ctgctgGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGAt                                 <  2:107930/149‑1 (MQ=255)                                                                                                             tgctgGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCGCTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGAt                                 >  2:101059/1‑148 (MQ=255)                                                                                                               ctgGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATa                                <  1:390348/147‑1 (MQ=255)                                                                                                                        ttCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCAc                     >  4:27913/1‑149 (MQ=255)                                                                                                                        ttCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCAc                     >  1:160942/1‑149 (MQ=255)                                                                                                                          ccACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCAcc                    <  4:66633/148‑1 (MQ=255)                                                                                                                            aCCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTc                 >  8:49721/1‑149 (MQ=255)                                                                                                                            aCCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTc                 >  7:91239/1‑149 (MQ=255)                                                                                                                            aCCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTc                 <  8:91239/149‑1 (MQ=255)                                                                                                                             cccGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTcc                <  1:151301/149‑1 (MQ=255)                                                                                                                               cgccgTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTCCTc              <  5:92708/149‑1 (MQ=255)                                                                                                                                  cgTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTCCTCg             <  6:139338/147‑1 (MQ=255)                                                                                                                                         cGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTCCTCGGCCAGCTcg    >  4:312742/1‑149 (MQ=255)                                                                                                                                           ttGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACCGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTCCTCGGCCAGCTcgcg  >  7:273287/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                  CGTATTCGCGGGTAAACCCGCTCCTACATCCGATCGGGGGTGGACACATGTTGGTGATTAGGCGCAGGCCGCGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCACTACTCAAACCCTTTCAAAAAGCCCCATGCACCCGCACGACCACCCCAAGAAATTCGCAACCCTCTTCGCACAACAGCGCCGGGTAGTTACTGTTCAATGCCTTCAGATAACGCTGCCCACCTTCCTCGGCCAGCTCGCG  >  NC_002947/4593832‑4594118

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