Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,395,850	A→G	E89G (GAG→GGG)	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,850	0	A	G	100.0%	167.2 / NA	49	E89G (GAG→GGG)	tsaE	N(6)‑L‑threonylcarbamoyladenine synthase, TsaE subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/25);  total (24/25)

CTTTAGCCGGGGCTTTTTACAGGCTCTGGGTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGCACAAGGCCGT  >  NC_000913/4395713‑4395983                                                                                                                                          |                                                                                                                                      cTTTAGCCGGGGCTTTTTACAGGCTCTGGGTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGc                                                                                                                                     >  2:475323/1‑140 (MQ=255)   ttAGCCGGGGCTTTTTACAGGCTCTGGGTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTg                                                                                                                                   <  1:272652/140‑1 (MQ=255)     aGCCGGGGCTTTTTACAGGCTCTGGGTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGa                                                                                                                                 <  1:400743/140‑1 (MQ=255)                            gggTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGattatt                                                                                                          >  1:373665/1‑140 (MQ=255)                                      ggTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGa                                                                                                <  1:120599/140‑1 (MQ=255)                                       gTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGAt                                                                                               >  2:16716/1‑140 (MQ=255)                                          aTGTCAAAAGCCCCACTTATTCGCTGGTCGAACCCGATACGCGCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCTGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGcc                                                                                            <  2:722219/140‑1 (MQ=255)                                           tGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCa                                                                                           >  1:427594/1‑140 (MQ=255)                                            gTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCAt                                                                                          <  2:113262/140‑1 (MQ=255)                                                aaaGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGc                                                                                      >  1:300172/1‑140 (MQ=255)                                                aaaGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGc                                                                                      >  2:375313/1‑140 (MQ=255)                                                aaaGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGAAAACTTAATGGGCTATCACTTTGATTTGTACCCCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGc                                                                                      >  1:722219/1‑140 (MQ=255)                                                        aCTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGa                                                                              <  1:375313/140‑1 (MQ=255)                                                            ataCGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTgg                                                                          >  1:264069/1‑140 (MQ=255)                                                             taCGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGc                                                                         <  1:16716/140‑1 (MQ=255)                                                              aCGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGcc                                                                        >  2:659012/1‑140 (MQ=255)                                                                  tGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGACCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaa                                                                    <  2:478330/140‑1 (MQ=255)                                                                   ggTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaac                                                                   <  1:192874/140‑1 (MQ=255)                                                                    gTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCTacaaca                                                                  >  1:573628/1‑140 (MQ=255)                                                                    gTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaaca                                                                  >  2:774097/1‑140 (MQ=255)                                                                    gTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaaca                                                                  >  2:551488/1‑140 (MQ=255)                                                                    gTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaaca                                                                  >  2:201767/1‑140 (MQ=255)                                                                    gTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaaca                                                                  >  2:115461/1‑140 (MQ=255)                                                                    gTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaaca                                                                  >  1:478349/1‑140 (MQ=255)                                                                    gTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaaca                                                                  >  1:478730/1‑140 (MQ=255)                                                                           ccTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTACCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACa                                                           <  1:264568/140‑1 (MQ=255)                                                                                     cGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGttcttc                                                 <  2:264069/140‑1 (MQ=255)                                                                                       aCAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCt                                               >  1:234849/1‑140 (MQ=255)                                                                                       aCAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCt                                               >  1:602394/1‑140 (MQ=255)                                                                                          aCTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGAc                                            <  2:358700/140‑1 (MQ=255)                                                                                              aaTGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCgg                                        >  2:584173/1‑140 (MQ=255)                                                                                              aaTGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCgg                                        <  2:300172/140‑1 (MQ=255)                                                                                                     tATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGa                                 >  2:338765/1‑140 (MQ=255)                                                                                                       tCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGTCCCGGATGTCGaaa                               <  1:201767/140‑1 (MQ=255)                                                                                                         aCTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATa                             <  1:774097/140‑1 (MQ=255)                                                                                                         aCTTTGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATa                             <  1:584173/140‑1 (MQ=255)                                                                                                           tttGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATaca                           <  2:602394/140‑1 (MQ=255)                                                                                                             tGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATAcaca                         <  2:478349/140‑1 (MQ=255)                                                                                                             tGATTTGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATAcaca                         >  1:139555/1‑140 (MQ=255)                                                                                                                 ttGTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGa                     >  1:224510/1‑140 (MQ=255)                                                                                                                   gTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACACGGATGTCGAAATa                             <  1:513300/130‑1 (MQ=255)                                                                                                                   gTACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACACGGATGTCGAAATa                             >  2:513300/1‑130 (MQ=255)                                                                                                                    tACCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTa                  <  2:224510/140‑1 (MQ=255)                                                                                                                     aCCGCCTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGTCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTAt                 <  2:611034/140‑1 (MQ=255)                                                                                                                         ccTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATc                >  2:688225/1‑137 (MQ=255)                                                                                                                         ccTTGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATc                <  1:688225/137‑1 (MQ=255)                                                                                                                           ttGCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGca           <  2:234849/140‑1 (MQ=255)                                                                                                                             gCCGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGcaca         <  1:115461/140‑1 (MQ=255)                                                                                                                              ccGATCCCGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGCACaa        >  2:360231/1‑140 (MQ=255)                                                                                                                                    ccGAGGGGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGCACAAGGCCGt  <  2:373665/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      CTTTAGCCGGGGCTTTTTACAGGCTCTGGGTCATCAGGGTAATGTCAAAAGCCCCACTTATACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTACAGGTGTTCTTCCTGACCCGGATGTCGAAATACACATTGATTATCAGGCACAAGGCCGT  >  NC_000913/4395713‑4395983

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