Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,683,251	C→G	15.8%	intergenic (+60/‑171)	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,251	0	C	G	15.8%	93.7 / 11.4	38	intergenic (+60/‑171)	recJ/PP_1478	ssDNA‑specific exonuclease RecJ/xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base C (15/17);  new base G (2/4);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

GTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCA  >  NC_002947/1683110‑1683394                                                                                                                                              |                                                                                                                                                 gTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                           >  5:346839/1‑149 (MQ=255)    ggCGTACAAGATGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTGGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                         <  4:152/148‑1 (MQ=255)     gCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGAGCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTC‑Ga                                                                                                                                      <  5:155197/149‑1 (MQ=255)         aCAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGAGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GAtct                                                                                                                                   >  4:163311/1‑148 (MQ=255)          cAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTC‑GAtttgt                                                                                                                                 <  1:67166/149‑5 (MQ=255)                  aTGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCt                                                                                                                         <  5:84162/149‑1 (MQ=255)                        aCGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAAc                                                                                                                   <  3:232733/149‑1 (MQ=255)                          gAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCt                                                                                                                 <  4:171552/149‑1 (MQ=255)                            gTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCttt                                                                                                               >  2:208390/1‑149 (MQ=255)                                       aTGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGACCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAAccc                                                                                                     <  5:236276/148‑1 (MQ=255)                                        tGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTC‑g                                                                                                                                       <  3:238495/113‑1 (MQ=255)                                        tGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTC‑g                                                                                                                                       >  4:238495/1‑113 (MQ=255)                                           gagTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCgcgg                                                                                                                                              >  2:182868/1‑104 (MQ=255)                                           gagTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGc                                                                                                <  3:163311/149‑1 (MQ=255)                                           gagTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGCGCACAAATCGAGCGCCGCCCGCgcgg                                                                                                                                              <  1:177371/104‑1 (MQ=255)                                           gagTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGCGCACAAATCGAGCGCCGCCCGCgcgg                                                                                                                                              <  1:182868/104‑1 (MQ=255)                                           gagTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGACTGCCTTGGATATTGCGCTGGGGCACAAATCGAGCGCCGCCCGCgcgg                                                                                                                                              >  2:177371/1‑104 (MQ=255)                                                 gCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATcc                                                                                           >  7:272642/1‑148 (MQ=255)                                                  cAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCca                                                                                         <  8:272642/149‑1 (MQ=255)                                                  cAGTTGATGATAGCGCATCTGGAGCCGCGCTGAGTTTAGCGTCAGGAGCCTGCCGTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCca                                                                                         <  5:182361/149‑1 (MQ=255)                                                                  atatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCt                                                                                 <  3:57147/141‑1 (MQ=255)                                                                  atatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCt                                                                                 >  4:57147/1‑141 (MQ=255)                                                                   tatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATcc                                                                                           >  3:79828/1‑130 (MQ=255)                                                                   tatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATcc                                                                                           <  4:79828/130‑1 (MQ=255)                                                                         gCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTg                                                                  <  2:141777/149‑1 (MQ=255)                                                                          ccGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCca                                                                                         <  4:101087/125‑1 (MQ=255)                                                                          ccGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCca                                                                                         >  3:101087/1‑125 (MQ=255)                                                                                       aGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCAACCCCCACCAATCTGGTCAAGTCTCa                                                       >  5:214152/1‑146 (MQ=255)                                                                                                        ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGc                                    >  4:281296/1‑148 (MQ=255)                                                                                                           gATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACGGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAgg                                >  2:323553/1‑149 (MQ=255)                                                                                                                     tGGGGCACAGATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACt                      >  8:144909/1‑149 (MQ=255)                                                                                                                     tGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCCCAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGCCCTGCACCGGGCCAGGGATGTGCAc                       >  1:227398/1‑149 (MQ=39)                                                                                                                       gggCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGACTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACtt                     <  6:346839/148‑1 (MQ=255)                                                                                                                        ggCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa                   <  3:207678/149‑1 (MQ=255)                                                                                                                                tCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGAc                                               >  2:289904/1‑113 (MQ=255)                                                                                                                                tCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGAc                                               <  1:289904/113‑1 (MQ=255)                                                                                                                                tCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCgg           <  1:66982/149‑1 (MQ=255)                                                                                                                                  gAGCGCCGCCCGCGCGGCGCTC‑GATATCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTCCACTTGAGAGTCCGTGc         >  2:73737/1‑149 (MQ=255)                                                                                                                                         gcccgcGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCCCAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCa  >  1:331835/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                 GTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGCCGGATCA  >  NC_002947/1683110‑1683394

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