Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,597,085	A→C	10.5%	intergenic (‑80/+16)	nemA ← / ← PP_3174	FMN‑linked N‑ethylmaleimide reductase/oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,597,085	0	A	C	10.5%	124.1 / 6.0	46	intergenic (‑80/+16)	nemA/PP_3174	FMN‑linked N‑ethylmaleimide reductase/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base C (3/2);  total (24/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 = 6.63e-01

AGGGGGGCCATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCC  >  NC_002947/3596938‑3597221                                                                                                                                                    |                                                                                                                                         cggggggCCATGAATCCGCGGTTGGGCAGGGTGAGTGGGCCGATTGGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAACGGGGGGACCGGGAGCCCGGCTCAg                                                                                                                                         <  8:96741/148‑1 (MQ=255)     gggCCATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCCGGGCt                                                                                                                                     >  2:86549/1‑149 (MQ=255)        ccATGATTAGGGGGTGGGGCTGGTGTAGTGTGCCGTTTTGCATGGGTTGAAAGATTTTCATTGCAGCTCCAGAGGTGGATATCCAATGTTGGGAACGATACCGGGGTGATTGGGGAATGGGGGGACCGGGAGCCCTGAGCCGGGCTccc                                                                                                                                  <  6:222292/149‑1 (MQ=255)         cttcaaTCCGCGGTTGGGCTGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTcccc                                                                                                                                 <  6:50151/145‑1 (MQ=255)         cATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTAAGGGCTccc                                                                                                                                  >  7:303241/1‑148 (MQ=255)               tACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGc                                                                                                                            >  2:144702/1‑148 (MQ=255)                tcgcgGGCGGGCCGGGTGAGTGGGCCTATTTGCCGCGGTTGCATCAGTTTCATTGCAGCTCCATCGGTCTAAAACCAATGGTGCGCACGCTTCCGGGGTTATTGGGGAACGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGt                                                                                                                          <  4:137963/148‑1 (MQ=255)                  gcgGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGt                                                                                                                          >  3:29361/1‑147 (MQ=255)                         ggCAGCGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAAAGCCTCCTCTGcat                                                                                                                 >  4:193859/1‑147 (MQ=255)                                    ggggCCGTTTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGGGGAAAACCTATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAgccgcc                                                                                                       <  4:212985/147‑1 (MQ=255)                                            tttGCAGCGGTGCCAAGACTTTCATTGCAGCGCCAGAGGTGCAAAACCAATGGTCCGCACGATCCCGGGGTGATTGGGCAAAGGGCGCACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTcac                                                                                                <  7:31650/147‑1 (MQ=255)                                               gccgCGGTTGCAAGAGTTCCATTGCAGCCCCAGAGGTGGAAACCCAATGGTGCGCACGATACCGGGGGGATTGGGGAAAGGGGGCACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATaa                                                                                            <  5:370917/145‑1 (MQ=255)                                               gCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAAt                                                                                           >  3:86921/1‑149 (MQ=255)                                               gCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGACGCCCACTCACATAAt                                                                                           >  5:68677/1‑149 (MQ=255)                                                        gCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGt                                                                                  >  1:223170/1‑149 (MQ=255)                                                          aaGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACGCACATAATGCGCGACGTCa                                                                                >  8:268265/1‑149 (MQ=255)                                                                   aTTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGccc                                                                       >  8:56149/1‑149 (MQ=255)                                                                    ttGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGccc                                                                       >  7:238286/1‑148 (MQ=255)                                                                         gCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTg                                                                 >  8:118529/1‑149 (MQ=255)                                                                             cAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGa                                                               <  6:100897/147‑1 (MQ=255)                                                                              agagGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACg                                                             >  4:162779/1‑148 (MQ=255)                                                                                   tGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGcc                                                       >  8:110597/1‑149 (MQ=255)                                                                                   tGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGc                                                        <  8:303241/148‑1 (MQ=255)                                                                                       aaaCCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTg                                                    <  1:144702/148‑1 (MQ=255)                                                                                        aaCCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGc                                                  >  3:239937/1‑149 (MQ=255)                                                                                           cAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGtt                                               <  4:86921/149‑1 (MQ=255)                                                                                                gTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGt                                                >  6:288496/1‑143 (MQ=255)                                                                                                gTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGt                                                <  5:288496/143‑1 (MQ=255)                                                                                                     cACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAg                                                                          >  4:114359/1‑112 (MQ=255)                                                                                                     cACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAg                                                                          <  3:114359/112‑1 (MQ=255)                                                                                                           aCCGGGGTTATTGGGGAATGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATg                                <  3:193859/148‑1 (MQ=255)                                                                                                             cGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAAt                             <  8:238286/149‑1 (MQ=255)                                                                                                             cGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCAGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAAt                             <  2:252600/149‑1 (MQ=255)                                                                                                             cGGGGTCATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAAt                             <  7:110312/149‑1 (MQ=255)                                                                                                              ggggTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAAt                             >  8:110312/1‑148 (MQ=255)                                                                                                                ggTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTg                                                    <  4:284315/123‑1 (MQ=255)                                                                                                                ggTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTg                                                    >  3:284315/1‑123 (MQ=255)                                                                                                                    aTTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCa                       <  4:239937/148‑1 (MQ=255)                                                                                                                    aTTGGGGAAACGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGCGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCa                       <  4:29361/148‑1 (MQ=255)                                                                                                                                ggggACCGGGAGCCCGGCTCAGGGCTCCCCAATGAGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTaa          >  4:254729/1‑149 (MQ=255)                                                                                                                                 gggACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAAt         <  7:110597/149‑1 (MQ=255)                                                                                                                                 gcattggggAGCCCTGAGCCGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTaa          >  4:228768/7‑148 (MQ=255)                                                                                                                                      ggggAGCCCTGAGCCGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAg                                                                          <  4:258325/78‑1 (MQ=38)                                                                                                                                      ggggAGCCCTGAGCCGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAg                                                                          >  3:258325/2‑79 (MQ=38)                                                                                                                                        ggAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTcgcc  >  7:273976/1‑149 (MQ=255)                                                                                                                                        ggAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTcgcc  >  1:322467/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                         AGGGGGGCCATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCC  >  NC_002947/3596938‑3597221

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