Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,579,597	A→G	H35R (CAT→CGT)	yhhZ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,579,597	0	A	G	100.0%	154.1 / NA	48	H35R (CAT→CGT)	yhhZ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/24);  total (24/24)

TACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAAT  >  NC_000913/3579509‑3579686                                                                                         |                                                                                          tACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAg                                                                                       >  1:2776776/1‑93 (MQ=255) tACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAg                                                                                       >  2:2816631/1‑93 (MQ=255) tACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAg                                                                                       <  2:3922555/93‑1 (MQ=255) tACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAg                                                                                       >  1:113904/1‑93 (MQ=255)     tGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATg                                                                                   >  1:2455112/1‑93 (MQ=255)       aCAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGaa                                                                                 >  1:2186330/1‑93 (MQ=255)        cAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGaaa                                                                                >  2:772625/1‑93 (MQ=255)         aGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAat                                                                               >  1:2342273/1‑93 (MQ=255)          gTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAata                                                                              <  2:4419492/93‑1 (MQ=255)             aCGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATAttt                                                                           >  2:661051/1‑93 (MQ=255)               gggAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTAc                                                                         >  2:3439263/1‑93 (MQ=255)               gggAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGCGAGGATGAAATATTTAc                                                                         >  1:244027/1‑93 (MQ=255)                ggAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACa                                                                        >  2:4072939/1‑93 (MQ=255)                  agaACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACAtt                                                                      >  1:4591529/1‑93 (MQ=255)                  agaACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACAtt                                                                      >  1:2934353/1‑93 (MQ=255)                      caacaaGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCa                                                                  <  2:1669162/93‑1 (MQ=255)                           aGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCtt                                                             >  2:1656392/1‑93 (MQ=255)                              aaGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTaaa                                                          >  1:3515432/1‑93 (MQ=255)                                   tctcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAataata                                                     <  2:1950568/93‑1 (MQ=255)                                     tcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATAtt                                                   >  1:3936858/1‑93 (MQ=255)                                        gCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaat                                                >  2:1878809/1‑93 (MQ=255)                                          aGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaataa                                              >  1:3017484/1‑93 (MQ=255)                                            gTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAataata                                            >  1:4422800/1‑93 (MQ=255)                                             ttGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATAc                                           >  2:547119/1‑93 (MQ=255)                                                tgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACggg                                        >  1:46178/1‑93 (MQ=255)                                                tgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACggg                                        <  2:2997051/93‑1 (MQ=255)                                                tgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACggg                                        <  2:2997053/93‑1 (MQ=255)                                                  ggACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGc                                      <  1:3758344/93‑1 (MQ=255)                                                  ggACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGc                                      <  2:2704470/93‑1 (MQ=255)                                                  ggACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGc                                      <  2:4599714/93‑1 (MQ=255)                                                    aCTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCtt                                    >  2:303601/1‑93 (MQ=255)                                                     cTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATCCGGGGCTTg                                   <  2:4086724/93‑1 (MQ=255)                                                            gTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTcaca                            <  2:902060/93‑1 (MQ=255)                                                            gTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTcaca                            <  2:902061/93‑1 (MQ=255)                                                              cTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGt                          >  1:4316611/1‑93 (MQ=255)                                                              cTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGt                          <  1:2160679/93‑1 (MQ=255)                                                                aCAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTc                        <  2:1284422/93‑1 (MQ=255)                                                                     tAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATgg                   <  1:2095828/93‑1 (MQ=255)                                                                      aaTCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGt                  <  2:1295642/93‑1 (MQ=255)                                                                       aTCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGta                 <  1:1737419/93‑1 (MQ=255)                                                                        tCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGtat                <  1:1813062/93‑1 (MQ=255)                                                                              gCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACAtt          <  2:448792/93‑1 (MQ=255)                                                                              gCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACAtt          >  1:2561781/1‑93 (MQ=255)                                                                              gCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACAtt          <  1:1392607/93‑1 (MQ=255)                                                                                agagCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACAtttt        <  1:3931323/93‑1 (MQ=255)                                                                                   gCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTa     <  1:274957/93‑1 (MQ=255)                                                                                   gCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTa     <  1:2639072/93‑1 (MQ=255)                                                                                   gCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTa     <  2:4560543/93‑1 (MQ=255)                                                                                   gCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTa     <  1:2639071/93‑1 (MQ=255)                                                                                      ggCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAaat  >  2:2304508/1‑93 (MQ=255)                                                                                         |                                                                                          TACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTAAAT  >  NC_000913/3579509‑3579686

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