Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,784,031	C→G	11.1%	A16P (GCG→CCG) ‡	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,784,031	0	C	G	11.1%	97.6 / 4.9	36	A16P (GCG→CCG) ‡	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (19/13);  new base G (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 = 7.79e-01

GCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCG  >  NC_002947/1783884‑1784169                                                                                                                                                    |                                                                                                                                             gCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGACTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGcc                                                                                                                                             >  7:237400/1‑149 (MQ=255)   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             >  4:175706/1‑149 (MQ=255)                                                                                                                               tCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAggg             >  3:226230/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)                                                                                                                                                    |                                                                                                                                             GCAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCTCGCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCAC‑CTGGGGCATCGCGTCACGTTCCTGTCG‑TGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCG  >  NC_002947/1783884‑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: ‑