Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,395,743	T→C	G53G (GGT→GGC)	tsaE →	N(6)‑L‑threonylcarbamoyladenine synthase, TsaE subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,395,743	0	T	C	100.0%	112.8 / NA	34	G53G (GGT→GGC)	tsaE	N(6)‑L‑threonylcarbamoyladenine synthase, TsaE subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/16);  total (18/16)

CGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGCACAAGGCCG  >  NC_000913/4395500‑4395982                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                cgTTTTGGAACGCGCGCGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGCATTGTTAACACGGAAGTGACTGATAACAACCATGATGAATCGAGTAATTCCGCTCCCTGCTGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCAGGGCCATCAgg                                                                                                                                                                                                                                          >  2:150550/1‑251 (MQ=255)          aCGCGCTGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGAGGAGCAGGCAACATTAGACCTGGGCTAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGTTAAAACCACCGTTAGCCGTGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCa                                                                                                                                                                                                                                 <  2:300308/251‑1 (MQ=255)                          cAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATAc                                                                                                                                                                                                                 >  1:305339/1‑251 (MQ=255)                                          cACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGAAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCtat                                                                                                                                                                                                 >  1:273970/1‑251 (MQ=255)                                                aCAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTc                                                                                                                                                                                           >  1:195677/1‑251 (MQ=255)                                                     gTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACa                                                                                                                                                                                       >  2:105471/1‑250 (MQ=255)                                                             aaCCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCACAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATgg                                                                                                                                                                              >  1:279359/1‑251 (MQ=255)                                                                   gAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTTGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATTGTCTATc                                                                                                                                                                        <  2:32310/251‑1 (MQ=255)                                                                       tGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACtt                                                                                                                                                                    >  1:31824/1‑251 (MQ=255)                                                                         aCTGATAAAAACCATGATGACTCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGCCAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTAGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTc                                                                                                                                                                                                          >  2:139768/1‑211 (MQ=255)                                                                         aCTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTCTACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGACTATCACTTTg                                                                                                                                                                  >  1:295111/1‑251 (MQ=255)                                                                         aCTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACAGTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTAGCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTGTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGGCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTg                                                                                                                                                                  <  1:202682/251‑1 (MQ=255)                                                                         aCTGATAAAAACCATGATGAATCGAGTAATGCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTc                                                                                                                                                                                                          <  1:139768/211‑1 (MQ=255)                                                                            gATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGAtt                                                                                                                                                               <  1:295728/251‑1 (MQ=255)                                                                            gATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGAtt                                                                                                                                                               <  1:211481/251‑1 (MQ=255)                                                                                                    aaTTCCTCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCTCTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCgagg                                                                                                                                       <  2:85932/251‑1 (MQ=255)                                                                                                    aaTTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCgagg                                                                                                                                       >  1:31122/1‑251 (MQ=255)                                                                                                    aaTTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCgag                                                                                                                                        <  1:27861/250‑1 (MQ=255)                                                                                                                 gatgaGCAGGCAACAGTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACGTTGATTTGTACCGCCTTGCCGATCCCGATGAGCTGGAGTTTa                                                                                                                           <  2:145325/250‑1 (MQ=255)                                                                                                                 gatgaGCAGGCAACAGTAGACCGGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTAGCTGGATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTTTGGGCCATCAGTGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCCCTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTa                                                                                                                           <  1:37/250‑1 (MQ=255)                                                                                                                    gaGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGTTTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATggg                                                                                                                       >  2:54094/1‑251 (MQ=255)                                                                                                                              cATTAGACCTGGGCGAGCGGGTAGCGAAAGTCTGCGATTGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCTGGGCTTTTTACAGGCTCTGGGCCATCATGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGAGTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGat                                                                                                             <  1:53519/251‑1 (MQ=255)                                                                                                                                                gggTAGCGAAAGCCTGCGATGGCGCAACCTTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCTCTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGcc                                                                                           <  1:300102/251‑1 (MQ=255)                                                                                                                                                           gCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTGTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCtggt                                                                                <  1:253788/251‑1 (MQ=255)                                                                                                                                                                   tGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGCGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGACTGGTGGAGTGGc                                                                        >  1:252621/1‑251 (MQ=255)                                                                                                                                                                            cGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTACCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAACGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaaca                                                                 >  1:181040/1‑249 (MQ=255)                                                                                                                                                                             gTAATATATCTGTATGGCGCTTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGCGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGACATCTGCCTGGTGGCGTGGCCACAACAAgg                                                              >  2:185854/1‑251 (MQ=255)                                                                                                                                                                                    agctgcacggCGATTTCGGCGCAAGTACAACCACATTTAGCCGGGGCTTTTTAGAGGCTCTCAGCCATCCGGGTAATGTAAACAGCCCCACTTCTCCGCTGGTAGACCCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGGACCGACTTGCCCACCCCGAGGAGCTGGCGTTTATGGGCATCAGCGATTCTTTTGCCACTGATGCCATCTGCCTGGTGGAGTGGCCACACCAAGGTACGGgt                                                       >  2:126735/9‑251 (MQ=255)                                                                                                                                                                                     tctGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGgtg                                                      >  2:187777/1‑251 (MQ=255)                                                                                                                                                                                        ggttGGCGAGTTAGGAGCAGGTAACACCGCCTTTAGCTGTGTCTTTGTTCAGGCTCTGGGCCATCAGGGTAATGTCAAAATCCCCATTTATACTCTGTTCGAACCCTATACGCTCGACAACTTAATGGTCTATCTCTTGGATTTGTACCGCCTTGCCTATCCCGAGGAGCTGAAGTTTATGTGGATCCTCGATTATTTTGCCGACGATGCCATGTGCCTGGTGGTGTGGCCAGAACAAGGTACATGTGtt                                                    <  2:302134/247‑1 (MQ=255)                                                                                                                                                                                                   tAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGAccc                                         <  1:100759/250‑1 (MQ=255)                                                                                                                                                                                                                aaCCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGTCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATac                           <  1:187777/251‑1 (MQ=255)                                                                                                                                                                                                                       tttAGCCGGGGCTTTTTACAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCCCTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGGGTTCTTCCTGACCCGGATGTCGAAATACACATTg                     >  2:255498/1‑250 (MQ=255)                                                                                                                                                                                                                                         cAGGCTCTGGGCCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGTAATACACATTGATTATCAGGCACAAGGCCg  >  2:293786/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                CGTTTTGGAACGCGCGGGATGCTGGCAACCGATCTCTTTTCCACGCTACAGCGTATTGTTAACCCGGAAGTGACTGATAAAAACCATGATGAATCGAGTAATTCCGCTCCCTGATGAGCAGGCAACATTAGACCTGGGCGAGCGGGTAGCGAAAGCCTGCGATGGCGCAACCGTAATCTATCTGTATGGCGATTTAGGCGCAGGTAAAACCACCTTTAGCCGGGGCTTTTTACAGGCTCTGGGTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGCACAAGGCCG  >  NC_000913/4395500‑4395982

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