Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,571,110	T→C	12.7%	intergenic (‑28/+16)	PP_2254 ← / ← dctA‑III	GntR family transcriptional regulator/C4‑dicarboxylate transport protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,571,110	0	T	C	12.7%	115.3 / 6.5	48	intergenic (‑28/+16)	PP_2254/dctA‑III	GntR family transcriptional regulator/C4‑dicarboxylate transport protein
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base C (3/3);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CAGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCACGGTGG  >  NC_002947/2570963‑2571253                                                                                                                                                       |                                                                                                                                                cAGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGGTGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑Tg                                                                                                                                                >  3:276629/1‑149 (MQ=255)    ttgttCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCTCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCg                                                                                                                                             <  6:230452/149‑1 (MQ=255)           gCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGACGGCGTCA‑TGCCGGCCACGt                                                                                                                                      >  8:130687/1‑149 (MQ=255)            ccAGTTCGCGCTGGCCTGGCAGCAGGGCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCATGA‑CGCCGGcccttcc                                                                                                                                     <  7:348769/149‑6 (MQ=255)             cAGTTGGCGGTGGCGTGGCAGCATGTGGGGGGTGGGCGATTGGCGGGGGGCTTGGGCTTGGGGGAGTTTTTCCACCATTGTG‑TTGA‑CGAGGGTTTATGAGCTGATCAGGTTTCGCTCTTGGAAGGGCGGGCGTCA‑TGCCGGCCACGTTc                                                                                                                                    <  8:290444/149‑1 (MQ=255)             cAGGTCGCGCGGGCCTGCCAGCATGGCGCCGGTGCGCCATGGGCCGCGGGCCTGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTCATCACGGTTCGCTCCTGGAAGGGCCGGCATGA‑CGCCGGcccttcca                                                                                                                                    <  6:231849/149‑7 (MQ=255)              aGTTCGCGCGGGCCGGGCCGCATGTCGCGGGTGCGCCATTGGCCGCGGGCCTGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTCA‑CGCCGTTTGTTGCGCTGATCCCGGTGCGCCCCTGTACGGGCCGGCGTCA‑TGCCGGCCACGTTCg                                                                                                                                   <  6:127739/149‑1 (MQ=255)                tgcgcgCTGGCCGGGCAGCATGGGGCCGGTGCGCCTTTGGCCGCGGGCATGGGCCGGGCGTAGTTTTTCCACCCCTGAG‑TTGC‑CGACGGTTGATGCGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCt                                                                                                                                 <  8:183328/147‑1 (MQ=255)                           ctggcAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTg                                                                                                                      >  7:22348/1‑149 (MQ=255)                            tggcAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGt                                                                                                                     >  1:342368/1‑149 (MQ=255)                                agcaTGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGAt                                                                                                                 >  2:345537/1‑149 (MQ=255)                                   aTGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATggg                                                                                                              <  1:231041/149‑1 (MQ=255)                                        gCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATGCCTGTTGATGGGCGCTg                                                                                                         >  7:191969/1‑149 (MQ=255)                                                      ggCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgtg                                                                                            <  3:81992/148‑1 (MQ=255)                                                        ccGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgtggcg                                                                                         >  1:281739/1‑149 (MQ=255)                                                        ccGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTgcgtggcg                                                                                         >  7:107547/1‑149 (MQ=255)                                                           cgGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGcgtggcgtgg                                                                                      <  5:295441/149‑1 (MQ=255)                                                           cgGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGcgtggcgtgg                                                                                      <  1:109917/149‑1 (MQ=255)                                                            gggCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTggg                                                                                     >  6:119519/1‑149 (MQ=255)                                                             ggCCAGGGCCTGGCGGAGTTTTTCCACCACTGAGTTTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTggg                                                                                     >  1:241452/1‑149 (MQ=255)                                                                  gggCCGGGCGTTGTTTTTCCACCCCTGA‑‑TTGACCGACGGTGCCTGAGGTTATCACGGTTCGCTCGTGCATGGGCCGGCGTCA‑TGCGGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCaaa                                                                               <  7:130687/149‑1 (MQ=255)                                                                       tGGCGGAGTTTTTCCCCCCCTGCG‑TTGC‑CGCCGGTGCATGCGCTGATCCCGGTGCGCTCCTGGACGGGCCGGCGTCT‑TGCGGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAg                                                                          <  1:138513/149‑1 (MQ=255)                                                                       tGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAg                                                                          >  5:351683/1‑149 (MQ=255)                                                                         gCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTGCTGTTGATGGGCGCTGGGTTCGTTGCGTGGCTTGGGGGCAAAGCCAGGc                                                                        >  5:217224/1‑149 (MQ=255)                                                                         gCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGc                                                                        <  2:241452/149‑1 (MQ=255)                                                                             aGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTTATGAGCTGATCACGTTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGc                                                                    <  2:209018/149‑1 (MQ=255)                                                                                 tttCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTc                                                                >  7:26980/1‑149 (MQ=255)                                                                                        cacTGAG‑TTTA‑CGACGGTTTATTAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACg                                                         <  2:281739/149‑1 (MQ=255)                                                                                             cg‑TTGA‑CGACGGTTGCTGAGCTGATCACGGTTCGCTCCTGGAAGGGCGGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCt                                                     <  6:346826/147‑1 (MQ=255)                                                                                                  ga‑cgacGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAg                                                                                                                              >  3:24602/1‑71 (MQ=255)                                                                                                  ga‑cgacGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGttgctt                                                                                                 >  1:366158/1‑100 (MQ=255)                                                                                                        cgGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGttg                                           >  5:49145/1‑148 (MQ=255)                                                                                                               tgaGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTc                                     >  6:213995/1‑148 (MQ=255)                                                                                                                 aGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAAt                                  <  5:213995/149‑1 (MQ=255)                                                                                                                  gCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATg                                 >  1:53824/1‑149 (MQ=255)                                                                                                                     gctcccGGTGCGCCCCTGGACGGG‑CGGCGGCACTGCGGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTcc                              <  2:342368/144‑1 (MQ=255)                                                                                                                     gATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTcc                              >  5:49217/1‑149 (MQ=255)                                                                                                                      aTCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑AGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGTTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCt                             >  2:75349/1‑149 (MQ=255)                                                                                                                         aCGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTtct                          >  4:95644/1‑149 (MQ=255)                                                                                                                           ggTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTcc                        >  7:326188/1‑149 (MQ=255)                                                                                                                                 cTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACgg                  <  8:22348/149‑1 (MQ=255)                                                                                                                                      aagggggCCGGCATGA‑CGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTtct                          <  1:148762/132‑1 (MQ=255)                                                                                                                                        aaGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATgg           <  8:107547/149‑1 (MQ=255)                                                                                                                                          tggCCGGCATGA‑CGCCGGCCACGTTCGCTCAGAATACCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCGCCCCACGGGCAATGGc          >  1:205661/2‑148 (MQ=255)                                                                                                                                           ggCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCa        >  5:350371/1‑149 (MQ=255)                                                                                                                                           ggCCGGCATGA‑CGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCa        >  6:383639/1‑149 (MQ=255)                                                                                                                                           ggCCGGCATGA‑CGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCa        >  6:226413/1‑149 (MQ=255)                                                                                                                                                 cttCA‑TGCCGGCCACGTTCGGTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCATAGCCAGTCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTGCTTCTCCCAACGGGCAATGGCCACGGTgg  <  6:217224/147‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                CAGTTGTTCGGCCAGTTCGCGCTGGCCTGGCAGCATGTCGCCGGTGCGCCATTGGCCGCGGGCCAGGGCCTGGCGGAGTTTTTCCACCACTGAG‑TTGA‑CGACGGTTGATGAGCTGATCACGGTTCGCTCCTGGAAGGGCCGGCGTCA‑TGCCGGCCACGTTCGCTCAGAATTCCTGTTGATGGGCGCTGGGTTGCTTGCGTGGCGTGGGGGCAAAGCCAGGCTTGCCGTCGAGCACGTTCTGCGCACGTTCCAGGTCAATGTCCTTCTCCCAACGGGCAATGGCCACGGTGG  >  NC_002947/2570963‑2571253

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