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%	99.5 / NA	30	H334Q (CAC→CAG)	relA	ATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (15/15);  total (15/15)

CGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAC  >  NC_002947/1851814‑1852093                                                                                                                                                 |                                                                                                                                        cgtgcgtgCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGc                                                                                                                                    <  1:888495/150‑1 (MQ=255)           gtgcgCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGc                                                                                                                          <  2:686642/150‑1 (MQ=255)           gtgcgCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGc                                                                                                                          >  1:686642/1‑150 (MQ=255)                cgTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGc                                                                                                                                    <  1:665402/135‑1 (MQ=255)                cgTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGc                                                                                                                                    >  2:665402/1‑135 (MQ=255)                                 tCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAgg                                                                                                    <  1:621760/150‑1 (MQ=255)                                       aCTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAGGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGa                                                                                              >  2:811997/1‑150 (MQ=255)                                         tGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATc                                                                                            <  2:380256/150‑1 (MQ=255)                                           cTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAgg                                                                                                    >  1:618173/1‑140 (MQ=255)                                           cTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAgg                                                                                                    <  2:618173/140‑1 (MQ=255)                                              cacCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTAc                                                                                       <  1:92770/150‑1 (MQ=255)                                                          cATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCa                                                                           >  2:673891/1‑150 (MQ=255)                                                           ctCGTGCACACCATGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGGCCCGAGAGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCAc                                                                          <  1:441898/149‑1 (MQ=255)                                                              gTGCACACCATGTGGCGGCATATTCCCAAGGAGTTCGAGGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGAGCTGTGCCGACGGCATGCACGag                                                                       >  2:441898/1‑149 (MQ=255)                                                                          tgGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACtt                                                           >  2:527253/1‑150 (MQ=255)                                                                             ctgcATATGCCCAAGGAGTTCGAGGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGc                                                        <  1:416313/148‑1 (MQ=255)                                                                             cggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGc                                                        >  1:525451/1‑150 (MQ=255)                                                                              ggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCg                                                       <  2:888071/150‑1 (MQ=255)                                                                                cATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTa                                                     <  2:525451/150‑1 (MQ=255)                                                                                   atTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATgc                                                  <  2:85865/150‑1 (MQ=255)                                                                                     tCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATgcgc                                                >  1:753624/1‑150 (MQ=255)                                                                                         aaGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCAc                                            >  2:60793/1‑150 (MQ=255)                                                                                              gTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGgcg                                       >  2:547778/1‑150 (MQ=255)                                                                                                  gacgacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTAc                                   >  1:830120/1‑150 (MQ=255)                                                                                                      acTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAgg                               <  1:130950/150‑1 (MQ=255)                                                                                                               ccAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGAGCACGTCCTGACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACg                      >  2:74143/1‑150 (MQ=255)                                                                                                                 aaCCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGTTCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTc                    <  1:60793/150‑1 (MQ=255)                                                                                                                   ccccAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGAAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCACAGACTGTCGCTACAAGGGCACCGACGTCaa                  >  2:214187/1‑150 (MQ=255)                                                                                                                                ggCTACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAAc     <  2:830120/150‑1 (MQ=255)                                                                                                                                   tACCGCTCGTTGCAGACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTNGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAc  >  1:72052/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                        CGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAC  >  NC_002947/1851814‑1852093

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