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

TACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTA  >  NC_000913/4395773‑4395924                                                                               |                                                                          taCGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACt                                                                         <  1:906065/81‑1 (MQ=255)  aCGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTg                                                                        >  1:906108/1‑81 (MQ=255)   cGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTgg                                                                       <  2:906067/81‑1 (MQ=255)   cGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTgg                                                                       <  2:906066/81‑1 (MQ=255)   cGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTg                                                                        >  1:906113/1‑80 (MQ=255)        gTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGttt                                                                  <  2:906068/81‑1 (MQ=255)        gTCGAACCCTATACGATCGACAACTTAATGGTATATCACTTTGATTTGTACCGCATTGCCGATCACGAGGAACTGGAGttt                                                                  <  2:979818/81‑1 (MQ=255)         tCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTa                                                                 >  1:906110/1‑81 (MQ=255)         tCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTa                                                                 >  2:906121/1‑81 (MQ=255)             aCCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATggg                                                             <  1:906069/81‑1 (MQ=255)                   taCGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATcc                                                        <  1:906070/80‑1 (MQ=255)                     cGCTCGACAACTTAATGGTCTATCACTTTGGTTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCcgc                                                      <  2:906071/80‑1 (MQ=255)                     cGCTCAACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGCTCCCGAGGAACTGGAGTTTATGGGGATCcgc                                                      <  2:906072/80‑1 (MQ=255)                          gACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGatt                                                  <  2:906073/79‑1 (MQ=255)                               cTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGc                                           <  2:906074/81‑1 (MQ=255)                                 tAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGcc                                          >  1:906116/1‑80 (MQ=255)                                      gTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGAt                                    >  2:906140/1‑81 (MQ=255)                                         tATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGTGATCCGCGATTATTTTGCCAACGATGc                                  >  2:906117/1‑80 (MQ=255)                                          aTCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCa                                >  1:906118/1‑81 (MQ=255)                                          aTCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCa                                >  1:906124/1‑81 (MQ=255)                                            cACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATc                              >  2:906098/1‑81 (MQ=255)                                              cTTTGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTg                            <  1:906075/81‑1 (MQ=255)                                                 tGATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCt                         >  2:906132/1‑81 (MQ=255)                                                  gATTTGTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCtg                        >  2:906120/1‑81 (MQ=255)                                                       gTACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCtggtg                     >  1:906119/1‑79 (MQ=255)                                                        tACCGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGa                   <  2:906076/80‑1 (MQ=255)                                                          ccGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGt                 <  1:906078/80‑1 (MQ=255)                                                          ccGCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTg                <  2:906077/81‑1 (MQ=255)                                                            gCCTTGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGc              >  1:906131/1‑81 (MQ=255)                                                               ttGCCGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCca            >  2:906128/1‑80 (MQ=255)                                                                  ccGATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaac        <  1:906079/81‑1 (MQ=255)                                                                    gATCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCacaac        >  2:906125/1‑79 (MQ=255)                                                                      tCCCGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAgg    >  2:906136/1‑81 (MQ=255)                                                                        ccGAGGAACTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGt   <  1:906080/80‑1 (MQ=255)                                                                          gaggaaCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTa  >  1:906123/1‑79 (MQ=255)                                                                          gaggaaCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTa  >  1:906146/1‑79 (MQ=255)                                                                               |                                                                          TACGCTGGTCGAACCCTATACGCTCGACAACTTAATGGTCTATCACTTTGATTTGTACCGCCTTGCCGATCCCGAGGAGCTGGAGTTTATGGGGATCCGCGATTATTTTGCCAACGATGCCATCTGCCTGGTGGAGTGGCCACAACAAGGTA  >  NC_000913/4395773‑4395924

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