Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,683,290	C→G	22.2%	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,290	0	C	G	22.2%	50.5 / 10.7	27	intergenic (+57/‑174)	recJ/PP_1478	ssDNA‑specific exonuclease RecJ/xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base C (11/10);  new base G (2/4);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GCTGGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGC  >  NC_002947/1683142‑1683429                                                                                                                                                     |                                                                                                                                            gcTGGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGccgccc                                                                                                                                            >  2:653264/1‑150 (MQ=255)   tGGGTGGAGGTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCATGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGccgcccgc                                                                                                                                          >  1:313757/1‑150 (MQ=255)           gTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                  <  2:96679/150‑1 (MQ=255)            ttGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                 <  1:653264/150‑1 (MQ=255)              ggCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                               <  1:346842/150‑1 (MQ=255)               gCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCgcc                                                                                                                              <  1:713759/150‑3 (MQ=255)                       gCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCGCAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAgg                                                                                                                      >  1:606063/1‑150 (MQ=255)                          ggATGTGAACGAGTTTCGCGGTAATGAGGCCGGGTAGGTGGTGATTGCGCATATGGAGCCGCTCTGAGTTTAGCGGCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGgcgc                                                                                                                   <  2:606063/150‑1 (MQ=38)                                           gcgGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGACGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGa                                                                                                  >  2:377185/1‑150 (MQ=255)                                                             aGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCcaca                                                                                >  1:404217/1‑150 (MQ=255)                                                                          gcATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCcca                                                                   >  1:283787/1‑150 (MQ=255)                                                                             tatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCCCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGaa                                                                             <  1:437222/137‑1 (MQ=255)                                                                             tatGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCCCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGaa                                                                             >  2:437222/1‑137 (MQ=255)                                                                                  aGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTg                                                           >  2:114004/1‑150 (MQ=255)                                                                                        gcTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAg                                                     >  2:207888/1‑150 (MQ=255)                                                                                          tGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGaa                                                                                                 >  1:445702/1‑104 (MQ=255)                                                                                          tGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGaa                                                                                                 <  2:445702/104‑1 (MQ=255)                                                                                                  gCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAATCCCTCGCCATCCCACAGAACCTCCCCCCACCATTCTGGTCTAGTCTCATTAAt                                           <  2:283787/150‑1 (MQ=255)                                                                                                       aGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCg                                      <  2:404217/150‑1 (MQ=255)                                                                                                                  ttGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCa                           >  1:270464/1‑150 (MQ=255)                                                                                                                      aTCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGa                       >  1:317909/1‑150 (MQ=255)                                                                                                                          tGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCGCCCCGCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAtgtg                   <  2:313757/150‑1 (MQ=255)                                                                                                                                 gggCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa            <  1:207888/150‑1 (MQ=255)                                                                                                                                          tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAt                      >  1:632711/1‑131 (MQ=255)                                                                                                                                          tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGAt                      <  2:632711/131‑1 (MQ=255)                                                                                                                                          tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCggg   >  1:593275/1‑150 (MQ=255)                                                                                                                                           cGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCCAACCTTTAATCAAGCCCCTCGCCATCCCACCGACCCGCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGc  <  1:377185/150‑1 (MQ=38)                                                                                                                                             aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa            <  2:315440/138‑1 (MQ=255)                                                                                                                                             aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGa            >  1:315440/1‑138 (MQ=255)                                                                                                                                                     |                                                                                                                                            GCTGGGTGGAGTTGGCGTACAAGCTGGATGTGAACGAGTTTCGCGGTAATGAGAGTGTGCAGTTGATGATTGCGCATATGGAGCCGCGCTGAGTTTAGCGTCAGGTGCCTGCCTTGGATCTTGCGCTGGGGCACAAATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGCAAACCTTTCGGCGAACCCCTCGCCATCCCACAGAACCTCCCCCCACCAATCTGGTCTAGTCTCATTAATGACCGGCACCGGGCCAGGGATGTGCACTTGAGTGTCCGGGC  >  NC_002947/1683142‑1683429

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