Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,784,033	T→C	11.5%	Q15R (CAG→CGG) ‡	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,784,033	0	T	C	11.5%	85.8 / 3.4	36	Q15R (CAG→CGG) ‡	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (18/13);  new base C (2/2);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

AGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCG  >  NC_002947/1783886‑1784169                                                                                                                                                    |                                                                                                                                           aGGATATTGTCGCCGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTg                                                                                                                                           <  5:152046/149‑1 (MQ=255)   gATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTgg                                                                                                                                          >  5:236321/1‑148 (MQ=255)         gTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGGGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTgg                                                                                                                                    <  6:50568/148‑1 (MQ=255)          tCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGc                                                                                                                                   >  2:176683/1‑148 (MQ=255)          tCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGc                                                                                                                                   >  4:261769/1‑148 (MQ=255)             acgaGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTGGc                                                                                                                                   <  6:143464/145‑1 (MQ=255)             acgaGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTGGc                                                                                                                                   >  5:143464/1‑145 (MQ=255)              cgaGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCAGGACCTGGCCACg                                                                                                                               >  7:39109/1‑148 (MQ=255)                   cAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGaa                                                                                                                         >  3:255282/1‑149 (MQ=255)                      gcgcGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACa                                                                                                                       <  8:181982/148‑1 (MQ=255)                          gggCCTCGATCAGCGGTCGGATGTCGCTGCCGGTGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCg                                                                                                                  <  1:268516/149‑1 (MQ=255)                              cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACa                                                                                                                       >  3:229746/1‑140 (MQ=255)                              cTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACa                                                                                                                       <  4:229746/140‑1 (MQ=255)                                       cGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTggg                                                                                                    <  6:236321/149‑1 (MQ=255)                                               tGTCTCTGACGGCGCGAAACGGTTTGGCGAGCAGCTCGCGGGCCAGGCGGATGACTTTCTTGAAGGCGGCGATGGGGCTGGCATTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATcgc                                                                                             <  8:237400/148‑1 (MQ=255)                                                gTCGCTGCCGGCGCGAAACCGTTTGTCGAGCACCTCGCCTGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAAc                                                                                                                        >  2:25798/1‑121 (MQ=255)                                                gTCGCTGCCGGCGCGAAACCGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAAc                                                                                                                        <  1:25798/121‑1 (MQ=255)                                                gTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGt                                                                                           >  5:126843/1‑149 (MQ=255)                                                                aaaCGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑Tgg                                                                          >  4:126503/1‑149 (MQ=255)                                                                       tGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGc                                                                    <  4:113632/148‑1 (MQ=255)                                                                         tCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTg                                                                 >  2:190204/1‑149 (MQ=255)                                                                            aGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCg                                                              <  3:261769/149‑1 (MQ=255)                                                                                 cTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTg                                                         >  2:80761/1‑149 (MQ=255)                                                                                   cGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCg                                                       <  3:126503/149‑1 (MQ=255)                                                                                             gCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCg                                              >  6:82172/1‑148 (MQ=255)                                                                                                      tttCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCgg                                    >  7:264/1‑149 (MQ=255)                                                                                                        tCTTGCAGGCGGCGATGGGGCGCGCGTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGGGCGGTCAGGCCGAAGTGCGCggg                                   <  5:270133/148‑1 (MQ=255)                                                                                                           tGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATgg                               >  8:225224/1‑149 (MQ=255)                                                                                                               ggcggcGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtc                           >  6:69920/1‑149 (MQ=255)                                                                                                                         gggCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTTAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTgcgc                 <  1:176683/149‑1 (MQ=255)                                                                                                                             tCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAggg             >  3:226230/1‑149 (MQ=255)                                                                                                                             tCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAggg             >  4:175706/1‑149 (MQ=255)                                                                                                                               gCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCGCTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAgg              <  1:80761/147‑1 (MQ=255)                                                                                                                                   ccgtggCCAGTTCCAGGCGGAACTGGCCCCGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTg           >  6:127161/5‑145 (MQ=255)                                                                                                                                       ggCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtctc   >  2:269916/1‑149 (MQ=255)                                                                                                                                         ccAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACGCTGGGGCAGCGCGTCACGTTCGTGGCG‑GGGGGTGGCGTGGCGGGGTGCGGGCAGGCCGAAGTGCGCGGGATGGTGGCGGCGTTGCGCAGGGTGAAGATCTCg  >  1:195526/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                           AGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCG  >  NC_002947/1783886‑1784169

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