Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,891,040	G→A	G76G (GGC→GGT)	tsaB ←	N(6)‑L‑threonylcarbamoyladenine synthase, TsaB subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,891,040	0	G	A	100.0%	51.1 / NA	17	G76G (GGC→GGT)	tsaB	N(6)‑L‑threonylcarbamoyladenine synthase, TsaB subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/10);  total (7/10)

CAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTC  >  NC_000913/1890922‑1891158                                                                                                                       |                                                                                                                       cAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATgcg                                                                                                          <  2:300117/133‑1 (MQ=255)   aTGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCa                                                                                                        <  1:294594/133‑1 (MQ=255)   aTGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCa                                                                                                        <  1:775/133‑1 (MQ=255)                               ttttGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGccc                                                                            >  1:98502/1‑133 (MQ=255)                                           ccccTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGc                                                                <  2:9795/133‑1 (MQ=255)                                                  cgcCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAata                                                         <  1:61774/133‑1 (MQ=255)                                                                 cGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGtt                                          >  1:41279/1‑133 (MQ=255)                                                                          gCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCtggtg                                 >  2:494757/1‑133 (MQ=255)                                                                           ccAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGGTAATATCAGTCTGGAAAGTTCCGCtggtgg                                >  1:203192/1‑133 (MQ=255)                                                                                   gggTAATGCCGCGCCAAGCCCAAGCGCTTGCGCGATACCGTTACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGGCAGGGAAGTTCCGCTGGTGGTCAGGata                        <  2:472451/132‑1 (MQ=255)                                                                                   cGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGata                        <  2:98502/133‑1 (MQ=255)                                                                                      tAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCATTCAGGGAAGTTCCGCTTGTGGACAGGATATcc                     >  2:170234/1‑133 (MQ=255)                                                                                      tAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGCTCCGCTGGCGGTCAGGATATcc                     >  1:472451/1‑133 (MQ=255)                                                                                            cgcgCCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCAcc               <  1:494757/133‑1 (MQ=255)                                                                                             gcgcCAAGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAATTCCGCTTGTGGTCAGGATATGCTGCACCa              >  1:194728/1‑133 (MQ=255)                                                                                                  aaGCGCCAGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGt         <  1:498363/133‑1 (MQ=255)                                                                                                     cgcCAGCCCTTGCGCGATACCAATCCCACTGCACACGCCAGGAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAg      <  2:194728/133‑1 (MQ=255)                                                                                                         aGCCCTTGCGCGATACCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTc  <  1:69202/133‑1 (MQ=255)                                                                                                                       |                                                                                                                       CAATGGCTGCCAGCACGCGGGTTGCGCCGTTTTTGCGCCACGCCCCTTGCGCCATCGTCATTAGCGTGGAGACGCCAATCATCGGTAATTCCGCGCCAAGCGCCAGCCCTTGCGCGATGCCAATACCAATGCGCACGCCAGTAAAGCTACCGGGGCCGCGCCCGTAAGCCAGAGCGTTAATATCAGTCAGGGAAGTTCCGCTGGTGGTCAGGATATCCTGCACCATCGGTAAGATTC  >  NC_000913/1890922‑1891158

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