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%	197.1 / NA	54	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 (30/24);  total (30/24)

CGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATC  >  NC_000913/4192453‑4192608                                                                               |                                                                              cgccgcGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAgctgact                                                                             >  2:854064/1‑81 (MQ=255)  gccgcGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAgctgactg                                                                            <  1:854019/81‑1 (MQ=255)  gccgcGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAgctgactg                                                                            >  2:854083/1‑81 (MQ=255)    cgcgATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGcc                                                                          <  2:854020/81‑1 (MQ=255)       gATAAGCGGTTCGAGGATAGCGTCGTTGTGCAGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGa                                                                       <  2:854021/81‑1 (MQ=255)           aGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAgc                                                                    <  2:854022/80‑1 (MQ=255)            gCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAgcgc                                                                  >  1:854093/1‑81 (MQ=255)             cGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGcgcg                                                                 >  1:854084/1‑81 (MQ=255)             cGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGcgcg                                                                 >  1:854086/1‑81 (MQ=255)              ggTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCgg                                                                <  2:854023/81‑1 (MQ=255)                  cGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCg                                                            >  2:854098/1‑81 (MQ=255)                   gAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGc                                                           >  2:854088/1‑81 (MQ=255)                   gAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGc                                                           >  1:854091/1‑81 (MQ=255)                   gAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGc                                                           >  1:854104/1‑81 (MQ=255)                    aGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGcc                                                          <  1:854025/81‑1 (MQ=255)                    aGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGcc                                                          <  2:854024/81‑1 (MQ=255)                     ggATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCt                                                         <  2:854026/81‑1 (MQ=255)                      gATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTc                                                        <  2:854027/81‑1 (MQ=255)                      gATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTc                                                        >  2:854096/1‑81 (MQ=255)                       aTAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTc                                                        >  1:854106/1‑80 (MQ=255)                         aGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTccac                                                     >  2:854110/1‑81 (MQ=255)                          gCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTccac                                                     >  2:854124/1‑80 (MQ=255)                           cgtcgtTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTccacca                                                   <  1:854028/81‑1 (MQ=255)                            gtcgtTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCAt                                                  >  2:854115/1‑81 (MQ=255)                               gtTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGATGACTGCCGGAAGCGCGGATCGCCTCCACCATCAg                                               <  1:854029/81‑1 (MQ=255)                                ttGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGt                                              <  1:854030/81‑1 (MQ=255)                                  gtgCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTtg                                            >  1:854094/1‑81 (MQ=255)                                   tgCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTtgt                                           <  2:854031/81‑1 (MQ=255)                                        gcgcAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTgaag                                       <  2:854032/80‑1 (MQ=255)                                         cgcAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTgaag                                       >  1:854105/1‑79 (MQ=255)                                          gcAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGaagaag                                    >  1:854108/1‑81 (MQ=255)                                                cGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAtt                              >  1:854103/1‑81 (MQ=255)                                                cGACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAtt                              >  2:854099/1‑81 (MQ=255)                                                 gACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAttt                             <  1:854033/81‑1 (MQ=255)                                                 gACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAttt                             >  1:854133/1‑81 (MQ=255)                                                 gACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAttt                             >  1:854123/1‑81 (MQ=255)                                                 gACACGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAAttt                             >  1:854107/1‑81 (MQ=255)                                                  acacGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTc                            <  2:854034/81‑1 (MQ=255)                                                    acGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTccc                          >  1:854111/1‑81 (MQ=255)                                                     cGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCt                         >  1:854117/1‑81 (MQ=255)                                                     cGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCt                         >  1:854122/1‑81 (MQ=255)                                                     cGTTTCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCt                         <  2:854035/81‑1 (MQ=255)                                                       tttCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCtgt                       <  1:854036/81‑1 (MQ=255)                                                       tttCATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCtgt                       <  1:854037/81‑1 (MQ=255)                                                          cATCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGcc                    >  1:854127/1‑81 (MQ=255)                                                           aTCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGcc                    <  1:854038/80‑1 (MQ=255)                                                            tCGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCgg                  <  2:854039/81‑1 (MQ=255)                                                             cGCCAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCgg                  >  1:854092/1‑80 (MQ=255)                                                               ccAGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATCTCCCTGTGCCGGTa                >  2:854130/1‑80 (MQ=255)                                                                 aGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAAc             <  2:854041/81‑1 (MQ=255)                                                                 aGTGTCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAAc             <  2:854040/81‑1 (MQ=255)                                                                  gtgtCACCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGACTTTCCCTGTGCCGGTaaa              >  2:854114/1‑79 (MQ=255)                                                                   tgtCACCAGCTGACTGCCGGAAGCGCGGATCACCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAg           <  2:854042/81‑1 (MQ=255)                                                                       aCCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGAtgt       >  1:854126/1‑81 (MQ=255)                                                                       aCCAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGAtgt       >  2:854100/1‑81 (MQ=255)                                                                        ccAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGAtgt       >  2:854102/1‑80 (MQ=255)                                                                         cAGCTGACTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGa     <  2:854043/81‑1 (MQ=255)                                                                            ctgactgCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATc  <  2:854044/81‑1 (MQ=255)                                                                               |                                                                              CGCCGCGATAAGCGGTTCGAGGATAGCGTCGTTGTGCTGGCGCAAGTCGACACGTTTCATCGCCAGTGTCACCAGCTGGCTGCCGGAAGCGCGGATCGCCTCCACCATCAGTTGTGAAGAAGCGAATTTCCCTGTGCCGGTAAACAGATGTGAATC  >  NC_000913/4192453‑4192608

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