Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,853,211	G→C	7.5%	intergenic (+56/‑68)	relA → / → mazG	ATP:GTP 3'‑pyrophosphotransferase/modified nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,853,211	0	G	C	7.5%	164.9 / 3.9	55	intergenic (+56/‑68)	relA/mazG	ATP:GTP 3'‑pyrophosphotransferase/modified nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base C (2/2);  total (26/29)
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.71e-01

GATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACCT  >  NC_002947/1853065‑1853355                                                                                                                                                    |                                                                                                                                                     gaTTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGa                                                                                                                                                     >  5:337965/1‑148 (MQ=255)     ccAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAAt                                                                                                                                                  <  1:80108/148‑1 (MQ=255)     ccAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCt                                                                                                                                                >  7:301621/1‑149 (MQ=255)         gCCTGGACGCGCTGGGGCGCCTGCAGGGGCGGATCTCGCAGTTGCCGAACATCATCCAGACGCGGCGTAATCGTACACCTTGAGCCCGGGGCGTATG‑CTTAGCGGGTAAACCCGCTGCTACAGGATGGAAGCAAGTTTGAAATCTGCAc                                                                                                                                            >  4:18242/1‑149 (MQ=255)             ggACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTa                                                                                                                                        >  1:143222/1‑149 (MQ=255)                  cgcTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCtgt                                                                                                                                   <  6:89701/149‑1 (MQ=255)                    cTGGGGCGCCTGCTGGGGCGGATCTCGGAGTTGCGGAACATCATCGAGACGCGGCGTAATCGTAGCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAAGCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTgg                                                                                                                                 <  7:4104/149‑1 (MQ=255)                       gggCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAg                                                                                                                              >  5:179223/1‑149 (MQ=255)                        ggCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGc                                                                                                                             <  5:264345/149‑1 (MQ=255)                           gcCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCggg                                                                                                                          <  6:337965/149‑1 (MQ=255)                                   gggCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCAcc                                                                                                                                           >  5:24226/1‑124 (MQ=255)                                   gggCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCAcc                                                                                                                                           <  6:24226/124‑1 (MQ=255)                                       ggATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCcgcg                                                                                                                >  6:104314/1‑147 (MQ=255)                                       ggATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTA‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCcgcg                                                                                                                >  8:202862/1‑147 (MQ=255)                                            tcGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGc                                                                                                                                              <  1:329032/112‑1 (MQ=255)                                              gCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGcca                                                                                                       >  4:110157/1‑149 (MQ=255)                                              gCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTCTGGGAGCGGGCAAGCCCGCGAAGAAGcc                                                                                                        <  1:200636/148‑1 (MQ=255)                                                      cGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGTAAGCCCGCGAAGAAGCCACCA‑CCg                                                                                                >  2:257564/1‑147 (MQ=255)                                                                      gcggcgTAATCGTAACCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGa                                                                            <  5:104314/149‑1 (MQ=255)                                                                          cgTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGTAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTa                                                                        <  1:257564/149‑1 (MQ=255)                                                                          cgTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTa                                                                        <  6:179223/149‑1 (MQ=255)                                                                              aTCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACAcc                                                                    <  4:304184/149‑1 (MQ=255)                                                                                    ccccTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCTCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTgac                                                              <  1:262985/149‑1 (MQ=255)                                                                                       cTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTgacga‑c                                                          >  3:152518/1‑149 (MQ=255)                                                                                          gagaCCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑cc                                                         <  8:301621/147‑1 (MQ=255)                                                                                           agaCCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑Cctgc                                                      >  2:290557/1‑149 (MQ=255)                                                                                             aCCGGGGCGTCTG‑CTTAGCGGGTTAACCCGCACCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑Cctgctg                                                    >  2:351394/1‑149 (MQ=255)                                                                                               cGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGGATGACCTACACCCTTGACGA‑CCTGCTGCa                                                  >  4:144585/1‑149 (MQ=255)                                                                                                 gggCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGAcc‑tgc                                                                                     >  3:312051/1‑114 (MQ=255)                                                                                                 gggCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGAcc‑tgc                                                                                     <  4:312051/114‑1 (MQ=255)                                                                                                         g‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGcc                                        >  1:124441/1‑149 (MQ=255)                                                                                                           cTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGccc                                       <  3:110157/149‑1 (MQ=255)                                                                                                                 gggTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTgc                                 <  7:202862/149‑1 (MQ=255)                                                                                                                      aaCCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGAc                            >  2:268616/1‑149 (MQ=255)                                                                                                                       aCCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGAcc                           >  1:175994/1‑149 (MQ=255)                                                                                                                         ccGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑Cctgc                                                      <  6:363537/120‑1 (MQ=255)                                                                                                                         ccGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑Cctgc                                                      >  5:363537/1‑120 (MQ=255)                                                                                                                          cGCTCCCACAGGTACGGTGCAGATTTCAAACCTGCATCCATCCTGTAGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCg                         >  5:7749/1‑148 (MQ=21)                                                                                                                          cGCTCCCACAGGTACGGTGCAGATTTCAAACCTGCATCCATCCTGTAGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCg                         <  6:7749/148‑1 (MQ=21)                                                                                                                                  cAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                 <  2:124441/148‑1 (MQ=255)                                                                                                                                  cAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                 <  1:198402/148‑1 (MQ=255)                                                                                                                                  cAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACg                  <  1:140188/147‑1 (MQ=255)                                                                                                                                      atggatgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTgcgc                               >  1:322032/1‑130 (MQ=255)                                                                                                                                      atggatgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTgcgc                               <  2:322032/130‑1 (MQ=255)                                                                                                                                         ggtgCAGATTTCAAACCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGAc                            <  1:269937/128‑1 (MQ=255)                                                                                                                                         ggtgCAGATTTCAAACCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGAc                            >  2:269937/3‑130 (MQ=255)                                                                                                                                          atgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGATCCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGt         <  1:268616/149‑1 (MQ=255)                                                                                                                                          atgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTg        >  1:48873/1‑149 (MQ=255)                                                                                                                                            gCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTgg       >  1:1020/1‑148 (MQ=255)                                                                                                                                             cAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGa     <  1:351394/149‑1 (MQ=255)                                                                                                                                             acggTTTTAAATATGCACCGTACCTGTGGTTGGGGGCAAGCCCGCGAATAATCCACCATCCCAATCCCCGGGCCACCTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGt         <  1:290557/146‑1 (MQ=255)                                                                                                                                              aGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                 <  5:94663/136‑1 (MQ=255)                                                                                                                                              aGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                 >  6:94663/1‑136 (MQ=255)                                                                                                                                               ggTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGAcc   <  2:175994/149‑1 (MQ=255)                                                                                                                                                 tttGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACCt  >  2:183815/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACCT  >  NC_002947/1853065‑1853355

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