Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,503,790	A→G	12.8%	intergenic (+39/‑59)	xanP → / → tusD	xanthine permease/sulfur transfer protein complex subunit TusD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,503,790	0	A	G	12.8%	102.4 / 8.2	49	intergenic (+39/‑59)	xanP/tusD	xanthine permease/sulfur transfer protein complex subunit TusD
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base G (3/3);  total (27/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 = 2.95e-01

CGCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGGG  >  NC_002947/4503642‑4503938                                                                                                                                                      |                                                                                                                                                       cgcACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGATTGAGAGTAGCCCGCGCCCGCTTGca                                                                                                                                                        <  5:96790/149‑1 (MQ=255)         ccGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgg                                                                                                                                                 <  3:178742/148‑1 (MQ=255)             cggcgATTCGCTTTGTGCGGGAGTCGGGTGTGCGTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGATTGAGAGTAGCCCGCGCCCGCTTGCACAAgcgggcgc‑g                                                                                                                                           <  6:343013/149‑1 (MQ=255)              ggcgATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgc‑gg                                                                                                                                          <  5:287296/149‑1 (MQ=255)                 ccatCGCAATGTGCTGGAGTCGGGGGTGGCTAGCGGGGGGATGAGTCCGCTGGGGGTGAACGGGTTCTTGCCGGAGGGGAAGCAGC‑CGCCTGAACCGGTGAGCCGGAATGAGAGTAGCCCGCGCCGGCGGGGACGAgcgggcgc‑gggcg                                                                                                                                       >  8:236373/4‑148 (MQ=255)                      gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCtttttt                                                                                                                                  >  4:66690/1‑149 (MQ=255)                      gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCtttttt                                                                                                                                  >  6:361317/1‑149 (MQ=255)                      gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCtttttt                                                                                                                                  >  3:72344/1‑149 (MQ=255)                        aaTGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGc                                                                                                                                >  3:253891/1‑149 (MQ=255)                        aaTGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGc                                                                                                                                >  7:175187/1‑149 (MQ=255)                          tgtgCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGTGGGCGC‑GGGCTTTTTTGCtt                                                                                                                              <  1:45369/149‑1 (MQ=255)                               tGGAGTCGGGTGTGGCTACCGGGGGGATTTCTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCa                                                                                                                         >  1:32913/1‑149 (MQ=255)                                 gAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCa                                                                                                                         <  8:117339/147‑1 (MQ=255)                                  aGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATgg                                                                                                                      <  7:286076/149‑1 (MQ=255)                                      gggTGTGGCTACCGAGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCGTTGTTGCGTTATCATGGCAc                                                                                                                   >  2:262106/1‑148 (MQ=255)                                            ggCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑t                                                                                                           >  4:75685/1‑149 (MQ=255)                                              cTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑ttt                                                                                                         <  4:72344/149‑1 (MQ=255)                                               tACCGGGGGGATTTGTGCGCTGATGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCATGCACAAGCGGGCGC‑GGGCtatt                                                                                                                                    >  7:234032/1‑122 (MQ=255)                                                  cGGGGGGATTTGTTCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCa                                                                                                      <  3:225242/148‑1 (MQ=255)                                                  cGGGGGGATTTGTGCGCTGGTGCTGAACTGGTGCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGACTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCAt                                                                                                     <  7:60031/149‑1 (MQ=255)                                                    gggggATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGa                                                                                                   <  5:153396/149‑1 (MQ=255)                                                         aTTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAg                                                                                              <  3:50729/149‑1 (MQ=255)                                                              tgCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCAt                                                                                          <  3:66690/148‑1 (MQ=255)                                                                             cTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATcgcg                                                                          <  4:253891/149‑1 (MQ=255)                                                                                gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgtg                                                                       <  4:113208/149‑1 (MQ=255)                                                                                gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCGTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgtg                                                                       >  5:334452/1‑149 (MQ=255)                                                                                gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCGTTTTTGCGTTATAAAGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTAGCGATCGCGgtg                                                                       >  6:111826/1‑149 (MQ=255)                                                                                  tCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCCCCAGCGGGCGC‑GGGCTTTTTTGCTTTATCAAGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGgtgt                                                                      >  8:80664/1‑148 (MQ=255)                                                                                  tCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGtt                                                                     >  1:95376/1‑149 (MQ=255)                                                                                       ccGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTcc                                                                 <  8:175187/148‑1 (MQ=255)                                                                                                 aGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATgc                                                      <  2:95376/149‑1 (MQ=255)                                                                                                         gCTAGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCGCtt                                              >  2:78856/1‑149 (MQ=255)                                                                                                           ttGAACC‑TTGAGTCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGa                                                                                                   <  3:239030/94‑1 (MQ=255)                                                                                                           ttGAACC‑TTGAGTCTGAATGAGAGTAGCCCGCGCCCCCATGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGTGACCATTCC‑ATTGCATGa                                                                                                   >  4:239030/1‑94 (MQ=255)                                                                                                           ttGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCCCGTGGCCAAGCGGGGGCGGGGGTTTTTTGGTTT‑TGAAGGGACCCGTGC‑GTTGGATGGGGTATGCATGAAAATAGCGTTAGCGGGGGTGTCCCCCGCCCCTGCGcccgcgtgg                                            >  2:353680/1‑143 (MQ=255)                                                                                                              aaCC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATAATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGa                                         >  3:45681/1‑149 (MQ=255)                                                                                                                        aCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTtc                                                                    <  4:221670/113‑1 (MQ=255)                                                                                                                        aCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTtc                                                                    >  3:221670/1‑113 (MQ=255)                                                                                                                           gAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGGGC‑GGGCTTTTTTGGTTTATCATGGGACCATTCC‑GTTGGATGAGTTAGCCATGGAATTAGCTATCGCGGTGTTGTCCCCGGCACATGCGCCCTCTTAGCGAAGCGCCCTgcgc                             >  2:25761/1‑149 (MQ=255)                                                                                                                                      ccccgcgcccgcTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTGATCATGGCACCATTCC‑TTTGCATGAGTTAGCCCTGAAATTCGCTCTCGCGGTGGTCTCCCCGGCCCCTGCGCCCTCTTCGCGAAGCGCAGGGCGCGTCGcgaagt                   <  7:236373/147‑6 (MQ=255)                                                                                                                                          gcgcccgcTTGTGCAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg              >  1:143576/1‑149 (MQ=255)                                                                                                                                          gcgcccgcTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg              >  4:105486/1‑149 (MQ=255)                                                                                                                                           cgcccgcTTGTGCAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGc             >  1:342255/1‑149 (MQ=255)                                                                                                                                            gcccgcTTGTGCAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg              <  7:20168/147‑1 (MQ=255)                                                                                                                                            gcccgcTTGTGCAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg              >  8:20168/1‑147 (MQ=255)                                                                                                                                              ccgcTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCCGTTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTg           >  1:276889/1‑149 (MQ=255)                                                                                                                                              ccgcTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTgg          <  4:45681/149‑1 (MQ=255)                                                                                                                                                  ttGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcgg      >  3:40511/1‑149 (MQ=255)                                                                                                                                                  ttGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcgg      >  1:234407/1‑149 (MQ=255)                                                                                                                                                    gCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggc     <  1:25761/148‑1 (MQ=255)                                                                                                                                                     cacaAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcgg   <  5:361317/149‑1 (MQ=255)                                                                                                                                                      acaAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCggg  >  1:62138/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CGCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACC‑TTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGC‑GGGCTTTTTTGCTTTATCATGGCACCATTCC‑TTTGCATGAGTTAGCCATGAAATTCGCTATCGCGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGGG  >  NC_002947/4503642‑4503938

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