Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,683,248	C→G	18.1%	intergenic (+57/‑174)	recJ → / → PP_1478	ssDNA‑specific exonuclease RecJ/xenobiotic reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,683,248	0	C	G	18.1%	85.8 / 10.7	38	intergenic (+57/‑174)	recJ/PP_1478	ssDNA‑specific exonuclease RecJ/xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base C (17/14);  new base G (2/5);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCAC  >  NC_002947/1683102‑1683395                                                                                                                                                     |                                                                                                                                                    tGGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGc‑‑cgcgcg                                                                                                                                                   >  8:152851/1‑149 (MQ=255)   ggTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCcgcg                                                                                                                                                 <  8:121336/149‑1 (MQ=255)          ttGGCGTTCAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCGCGGGCGgcgc                                                                                                                                           <  8:178791/148‑1 (MQ=255)              cGTACAAGCTGGATGTGAACGAGTTTCGCTGTAATTAGAGTGTGCATTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTTCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCGCGGGCGGCGCTCGAt                                                                                                                                      <  4:191013/149‑1 (MQ=38)               gTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCGCGGGCGGCGCTCGAtt                                                                                                                                     <  7:152851/149‑2 (MQ=255)                   aaGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTca                                                                                                                                  >  8:185269/1‑148 (MQ=255)                              gAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCGCGGGCGGCGCTCGAt                                                                                                                                      <  1:119896/133‑1 (MQ=255)                              gAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCGCGGGCGGCGCTCGAt                                                                                                                                      >  2:119896/1‑133 (MQ=255)                                    gTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCgcggcg                                                                                                                                            <  8:15249/121‑1 (MQ=255)                                    gTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCgcggcg                                                                                                                                            >  7:15249/1‑121 (MQ=255)                                     tttCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCg                                                                                                                                        >  5:32640/1‑124 (MQ=255)                                     tttCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCg                                                                                                                                        <  6:32640/124‑1 (MQ=255)                                                  agagTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTc                                                                                                               >  8:79598/1‑136 (MQ=255)                                                  agagTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTc                                                                                                               <  7:79598/136‑1 (MQ=255)                                                      tgtgCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAGATCGAGCGC‑‑CGCGGGGGCGGCGCTCGAtttgt                                                                                                                                  >  5:149641/1‑109 (MQ=255)                                                      tgtgCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAGATCGAGCGC‑‑CGCGCGGGCGGCGCTCGAtttgt                                                                                                                                  <  6:149641/113‑5 (MQ=38)                                                          cAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCca                                                                                          <  2:24392/149‑1 (MQ=255)                                                           aGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGaa                                                                                                         <  3:211528/133‑1 (MQ=255)                                                           aGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGaa                                                                                                         >  4:211528/1‑133 (MQ=255)                                                                     tGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGc                                                                                                            >  8:214954/1‑120 (MQ=255)                                                                              ggAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCGAAGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCCGCGAACCCCTCGCCATCCCACAGAACCTCCCCccacc                                                                         >  1:80509/1‑148 (MQ=255)                                                                              ggAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAAt                                                                      >  5:201719/1‑149 (MQ=255)                                                                                aGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCt                                                                    >  3:207887/1‑149 (MQ=255)                                                                                 gCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTg                                                                   <  7:185269/149‑1 (MQ=255)                                                                                   cgcgCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGt                                                                 <  4:207887/149‑1 (MQ=255)                                                                                    gcgcTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCg                                                                                                                                        >  8:196290/1‑77 (MQ=255)                                                                                    gcgcTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCg                                                                                                                                        <  7:196290/77‑1 (MQ=255)                                                                                                gCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTaa                                                    >  4:68546/1‑149 (MQ=255)                                                                                                gCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTaa                                                    >  2:103692/1‑149 (MQ=255)                                                                                                 cGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAAt                                                   <  3:68546/149‑1 (MQ=255)                                                                                                 cGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAAt                                                   <  1:103692/149‑1 (MQ=255)                                                                                                   tCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGGAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAAt                                                   >  1:84829/1‑147 (MQ=255)                                                                                                   tCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAAt                                                   <  2:84829/147‑1 (MQ=255)                                                                                                          cctgcctTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCAc                                          >  6:20306/1‑149 (MQ=255)                                                                                                                  ggATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCa                                   <  5:149073/148‑1 (MQ=255)                                                                                                                  ggATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAg                                  >  1:242445/1‑149 (MQ=255)                                                                                                                            cTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCAc                        <  2:242445/149‑1 (MQ=255)                                                                                                                                    caAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGCACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgtgt                >  8:75390/1‑149 (MQ=255)                                                                                                                                     aaaTCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCCCAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGCACCTCCCCCCACCAATCTGGGCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAgtggc               >  7:217234/1‑147 (MQ=255)                                                                                                                                                     ccgcGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCAc  >  7:56807/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGC‑‑CGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCAC  >  NC_002947/1683102‑1683395

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