Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,850,158	C→T	10.6%	intergenic (+75/‑112)	hemN → / → fnrA	oxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,850,158	0	C	T	10.6%	122.9 / 5.4	47	intergenic (+75/‑112)	hemN/fnrA	oxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base T (3/2);  total (26/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 = 9.87e-01

TGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACA  >  NC_002947/4850021‑4850298                                                                                                                                          |                                                                                                                                             tGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTGCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCgg                                                                                                                                    <  3:91286/148‑1 (MQ=255)                      cTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACcat                                                                                                               >  1:322132/1‑148 (MQ=255)                               cacaACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGca                                                                                                      >  3:13401/1‑148 (MQ=255)                                acaACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGACAGCGAGCACCATCATCGCGcaca                                                                                                    >  7:24667/1‑149 (MQ=255)                                   aCCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTg                                                                                                 <  2:322132/149‑1 (MQ=255)                                            gTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACAcc                                                                                        <  6:116977/149‑1 (MQ=255)                                                   gcgAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                  >  4:95189/1‑99 (MQ=255)                                                   gcgAGTGATCTGAAAAACAAGGCGTGCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCtgt                                                                                    <  6:165225/146‑1 (MQ=255)                                                    cgAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCt                                                                                      >  5:200033/1‑143 (MQ=255)                                                    cgAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCt                                                                                      <  6:200033/143‑1 (MQ=255)                                                          aTCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGt                                                                           >  8:180435/1‑148 (MQ=255)                                                             tGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAcc                                                                       >  4:279060/1‑149 (MQ=255)                                                             tGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAc                                                                        >  2:219575/1‑148 (MQ=255)                                                              gAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCgg                                                                                                                                    >  8:190792/1‑87 (MQ=255)                                                              gAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCgg                                                                                                                                    <  7:190792/87‑1 (MQ=255)                                                               aaaaaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAgggtaa                                                                       <  4:13401/147‑6 (MQ=255)                                                               aaaaaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAccc                                                                      <  1:201701/148‑1 (MQ=255)                                                               aaaaaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAc                                                                        <  1:219575/146‑1 (MQ=255)                                                                       ggCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACt                                                              >  1:233500/1‑148 (MQ=255)                                                                        gCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTa                                                            <  2:233500/149‑1 (MQ=255)                                                                          gTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtg                                                          >  4:233101/1‑149 (MQ=255)                                                                              aaCCGCATGGACAAGCTTCCTTGTCGGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtggtag                                                      >  4:84709/1‑145 (MQ=255)                                                                              aaCCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgtgt                                                       <  3:233101/148‑1 (MQ=255)                                                                                      ggACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTAGACCCTGTGTCGTGAGTTACCGTTCCGAGTGATGTGTGGTTTCCca                                              >  8:232333/1‑149 (MQ=255)                                                                                          aaGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTCCACCCTGTGTCGTg                                                                              <  7:301072/113‑1 (MQ=255)                                                                                          aaGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTg                                                                              >  8:301072/1‑113 (MQ=255)                                                                                             cTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCatcatc                                                                                                           <  2:86835/81‑1 (MQ=255)                                                                                              ttCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATc                                      <  8:24667/149‑1 (MQ=255)                                                                                                ccTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGa                                    >  7:145192/1‑149 (MQ=255)                                                                                                  ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                             <  7:162503/106‑1 (MQ=255)                                                                                                  ttGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                             >  8:162503/1‑106 (MQ=255)                                                                                                    gTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                 >  8:112868/1‑148 (MQ=255)                                                                                                     tCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCa                                <  8:145192/148‑1 (MQ=255)                                                                                                         gCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGt                           >  3:166143/1‑149 (MQ=255)                                                                                                             acTATGAGGACACTAGGACTCCCGGGAGGCCTCCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGa                       <  1:150874/149‑1 (MQ=255)                                                                                                             acTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGa                       <  7:112868/149‑1 (MQ=255)                                                                                                                      acacaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCa               >  2:274433/1‑148 (MQ=255)                                                                                                                        acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgt                                                         <  1:12359/104‑1 (MQ=255)                                                                                                                        acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgt                                                         >  2:12359/1‑104 (MQ=255)                                                                                                                        acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGc                     >  2:3937/1‑140 (MQ=255)                                                                                                                        acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGc                     <  1:3937/140‑1 (MQ=255)                                                                                                                        acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAAc            >  1:113651/1‑149 (MQ=255)                                                                                                                               accccGGCAGGTGGCCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGcc      >  1:122108/3‑148 (MQ=255)                                                                                                                                tCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCca    >  6:277385/1‑149 (MQ=255)                                                                                                                                  ccGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCcaca  >  7:242283/1‑149 (MQ=255)                                                                                                                                    ggCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCcac                                             >  4:229220/1‑104 (MQ=255)                                                                                                                                    ggCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCcac                                             <  3:229220/104‑1 (MQ=255)                                                                                                                                          |                                                                                                                                             TGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACA  >  NC_002947/4850021‑4850298

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