Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,851,958	C→G	H334Q (CAC→CAG)	relA →	ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,851,958	0	C	G	100.0%	120.8 / NA	37	H334Q (CAC→CAG)	relA	ATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (19/18);  total (19/18)

TGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAA  >  NC_002947/1851831‑1852103                                                                                                                                |                                                                                                                                                  tGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAgg                                                                                                                             <  1:926220/150‑1 (MQ=255) tGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAgg                                                                                                                             <  1:738602/150‑1 (MQ=255)      gTGCCGGAAATCCGCGACTGCTACACTGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCTAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAg                                                                                                                        <  2:143847/150‑1 (MQ=255)          cGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGc                                                                                                                    >  2:206208/1‑150 (MQ=255)           ggAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGc                                                                                                                                    >  1:157143/1‑133 (MQ=255)           ggAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGc                                                                                                                                    <  2:157143/133‑1 (MQ=255)                  cgcgACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTg                                                                                                            >  1:922484/1‑150 (MQ=255)                   gcgACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGc                                                                                                           >  2:449091/1‑150 (MQ=255)                              acCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAATGTGCTGGAGGTGCAGATCAGTAcc                                                                                                >  2:177240/1‑150 (MQ=255)                                cgcgCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCa                                                                                              >  1:343222/1‑150 (MQ=255)                                   gcTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACgg                                                                                           <  2:343222/150‑1 (MQ=255)                                      gggCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCAt                                                                                        <  1:437040/150‑1 (MQ=255)                                      gggCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCAt                                                                                        <  2:922484/150‑1 (MQ=255)                                         cATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCa                                                                                     >  1:793868/1‑150 (MQ=255)                                          aTCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCAc                                                                                    >  1:558411/1‑150 (MQ=255)                                          aTCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCAc                                                                                    <  2:1007535/150‑1 (MQ=255)                                           tCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCAt                                                                                        <  1:207070/145‑1 (MQ=255)                                           tCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCAt                                                                                        >  2:207070/1‑145 (MQ=255)                                                 acacCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAgg                                                                             >  1:115440/1‑150 (MQ=255)                                                    ccTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCg                                                                          <  2:115440/150‑1 (MQ=255)                                                          ggcggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTg                                                                    <  2:348550/150‑1 (MQ=255)                                                          ggcggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTg                                                                    <  2:793868/150‑1 (MQ=255)                                                             ggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCg                                                                 >  2:422151/1‑150 (MQ=255)                                                                atatTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTAt                                                              >  1:193661/1‑150 (MQ=255)                                                                       cAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGGACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCa                                                       >  1:453363/1‑150 (MQ=255)                                                                                        aCATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCa                                      <  1:422151/150‑1 (MQ=255)                                                                                         cATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCAc                                     <  1:206208/150‑1 (MQ=255)                                                                                          aTCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCAcc                                    >  2:562687/1‑150 (MQ=255)                                                                                           tCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCg                                   >  1:521165/1‑150 (MQ=255)                                                                                             gCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGAc                                 <  1:177240/150‑1 (MQ=255)                                                                                               cAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAAc                                                                       <  1:1025856/110‑1 (MQ=255)                                                                                               cAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAAc                                                                       >  2:1025856/1‑110 (MQ=255)                                                                                                              cGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCGCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCaa                >  1:769702/1‑150 (MQ=255)                                                                                                                  tACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCa                                                                                     <  2:943914/77‑1 (MQ=255)                                                                                                                  tACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCa                                                                                     >  1:943914/1‑77 (MQ=255)                                                                                                                           ttGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGAGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACgaagaa   >  2:274667/1‑150 (MQ=255)                                                                                                                            tGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGaaa  <  2:558411/150‑1 (MQ=255)                                                                                                                                |                                                                                                                                                  TGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAA  >  NC_002947/1851831‑1852103

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