Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,192,531	G→A	S35S (AGC→AGT)	thiG ←	1‑deoxy‑D‑xylulose 5‑phosphate:thiol sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,192,531	0	G	A	100.0%	112.0 / NA	31	S35S (AGC→AGT)	thiG	1‑deoxy‑D‑xylulose 5‑phosphate:thiol sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/16);  total (15/16)

CGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAAT  >  NC_000913/4192453‑4192607                                                                               |                                                                             cgccgcGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAgctgact                                                                            <  1:786333/81‑1 (MQ=255)    cgcgATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGcc                                                                         <  2:786333/81‑1 (MQ=255)    cgcgATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGcc                                                                         >  1:786389/1‑81 (MQ=255)     gcgATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCg                                                                        >  1:786383/1‑81 (MQ=255)       gATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCgg                                                                       >  1:786387/1‑80 (MQ=255)       gATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGa                                                                      <  1:786335/81‑1 (MQ=255)                    aGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGcc                                                         >  1:786393/1‑81 (MQ=255)                       aTAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTc                                                       >  2:786377/1‑80 (MQ=255)                          gCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTccacc                                                   <  2:786335/81‑1 (MQ=255)                            gtcgtTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCAt                                                 >  2:786397/1‑81 (MQ=255)                               gtTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAg                                              >  1:786436/1‑81 (MQ=255)                                  gtgCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGtt                                            <  2:786336/80‑1 (MQ=255)                                     cTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTga                                        <  2:786337/81‑1 (MQ=255)                                      tGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTtgtg                                         >  2:786393/1‑79 (MQ=255)                                       ggCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTtgtg                                         >  2:786398/1‑78 (MQ=255)                                        gcgcAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTgaaga                                     <  1:786339/81‑1 (MQ=255)                                           cAAGTCGACACGTTTCATCGCCAGTGTCACCGGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGc                                  <  2:786339/81‑1 (MQ=255)                                               tCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAt                              <  1:786341/81‑1 (MQ=255)                                                  acacGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTc                           <  2:786341/81‑1 (MQ=255)                                                    acGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTccc                         >  2:786409/1‑81 (MQ=255)                                                           aTCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGcc                   <  1:786343/80‑1 (MQ=255)                                                             cGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCgg                 >  2:786394/1‑80 (MQ=255)                                                              gCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGt                <  1:786344/80‑1 (MQ=255)                                                              gCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTa               <  1:786345/81‑1 (MQ=255)                                                                 aGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAAc            <  2:786345/81‑1 (MQ=255)                                                                  gtgtCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTaaa             >  2:786407/1‑79 (MQ=255)                                                                  gtgtCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACa           <  2:786346/81‑1 (MQ=255)                                                                    gtCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTaaa             >  2:786391/1‑77 (MQ=255)                                                                     tCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGAt        <  1:786348/81‑1 (MQ=255)                                                                       aCCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGAtgt      >  1:786401/1‑81 (MQ=255)                                                                       aCCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGAtgt      >  2:786415/1‑81 (MQ=255)                                                                       aCCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGAtg       <  1:786349/80‑1 (MQ=255)                                                                           gctgactgCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAAt  <  2:786349/81‑1 (MQ=255)                                                                               |                                                                             CGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAAT  >  NC_000913/4192453‑4192607

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