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%	229.4 / NA	72	H35R (CAT→CGT)	yhhZ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/38);  total (34/38)

TTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTA  >  NC_000913/3579507‑3579683                                                                                           |                                                                                       tttACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTg                                                                                      >  2:3134789/1‑93 (MQ=255)   tACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAg                                                                                    <  2:887241/93‑1 (MQ=255)    aCCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAgg                                                                                   >  2:5662093/1‑93 (MQ=255)     ccTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGa                                                                                  <  1:3776724/93‑1 (MQ=255)      cTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGAt                                                                                 >  2:1434124/1‑93 (MQ=255)       tGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGAGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATg                                                                                >  1:4086089/1‑93 (MQ=255)        gACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGa                                                                               <  2:5794697/93‑1 (MQ=255)           aGTAACGGGAGAACAACAAGGAAGCATCTCCGCGGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAat                                                                            >  2:5247131/1‑93 (MQ=255)           aGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAat                                                                            >  2:4127485/1‑93 (MQ=255)            gTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAata                                                                           >  2:3625227/1‑93 (MQ=255)            gTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAata                                                                           >  1:3125850/1‑93 (MQ=255)               aCGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATAttt                                                                        >  1:2709619/1‑93 (MQ=255)                cGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTa                                                                       >  2:423510/1‑93 (MQ=255)                cGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTa                                                                       <  2:3249802/93‑1 (MQ=255)                  ggAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACa                                                                     >  1:2625711/1‑93 (MQ=255)                  ggAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACa                                                                     <  2:1673912/93‑1 (MQ=255)                   gagaACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACAt                                                                    <  2:1829490/93‑1 (MQ=255)                     gaACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATtc                                                                  >  1:4646245/1‑93 (MQ=255)                     gaACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATtc                                                                  <  1:1570877/93‑1 (MQ=255)                       acaacaAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATtctc                                                                >  2:3727605/1‑93 (MQ=255)                        caacaaGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCa                                                               <  2:4707926/93‑1 (MQ=255)                        caacaaGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCa                                                               <  1:84890/93‑1 (MQ=255)                          acaaGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCAct                                                             >  1:1673850/1‑93 (MQ=255)                            aaGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCActct                                                           >  2:2152631/1‑93 (MQ=255)                                aaGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTaaa                                                       <  1:6011878/93‑1 (MQ=255)                                  gCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCATTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAaata                                                     >  1:1579267/1‑93 (MQ=255)                                   cATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAaataa                                                    <  1:2424093/93‑1 (MQ=255)                                   cATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAaataa                                                    <  2:5056989/93‑1 (MQ=255)                                    aTCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAaataat                                                   <  1:86823/93‑1 (MQ=255)                                     tctcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAataata                                                  >  1:5183588/1‑93 (MQ=255)                                     tctcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAataata                                                  <  2:4192910/93‑1 (MQ=255)                                     tctcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAataata                                                  >  2:5366199/1‑93 (MQ=255)                                     tctcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAataata                                                  >  2:4279743/1‑93 (MQ=255)                                     tctcCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAataata                                                  <  2:4742359/93‑1 (MQ=255)                                          gCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaat                                             >  1:6024882/1‑93 (MQ=255)                                          gCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaat                                             <  1:5071368/93‑1 (MQ=255)                                           cAGGTTGTGTGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaata                                            <  1:2430836/93‑1 (MQ=255)                                           cAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATtaata                                            <  2:2016000/93‑1 (MQ=255)                                              gTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAataata                                         >  1:5386260/1‑93 (MQ=255)                                              gTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAataata                                         >  2:1320165/1‑93 (MQ=255)                                               ttGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATAc                                        >  2:1999954/1‑93 (MQ=255)                                               ttGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATAc                                        <  2:3611851/93‑1 (MQ=255)                                                 gtgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACgg                                      >  1:3943824/1‑93 (MQ=255)                                                 gtgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACgg                                      <  2:3248184/93‑1 (MQ=255)                                                  tgGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACggg                                     <  1:4283022/93‑1 (MQ=255)                                                     gACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCt                                  <  2:3094308/93‑1 (MQ=255)                                                         tCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGt                              <  1:5508342/93‑1 (MQ=255)                                                         tCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGt                              <  2:739957/93‑1 (MQ=255)                                                            gAGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTca                           >  2:5675177/1‑93 (MQ=255)                                                             aGTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTcac                          <  2:4981057/93‑1 (MQ=255)                                                              gTCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTcaca                         <  1:3867737/93‑1 (MQ=255)                                                               tCTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAg                        >  2:345408/1‑93 (MQ=255)                                                                cTACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGt                       <  1:5518772/93‑1 (MQ=255)                                                                 tACAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGtt                      <  2:6201627/93‑1 (MQ=255)                                                                   cAGGTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTcc                    <  2:2047373/93‑1 (MQ=255)                                                                     ggTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCAt                  >  1:2448802/1‑93 (MQ=255)                                                                     ggTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCAt                  <  1:5047185/93‑1 (MQ=255)                                                                     ggTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCAt                  <  2:4421785/93‑1 (MQ=255)                                                                      gTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATg                 <  1:904706/93‑1 (MQ=255)                                                                      gTAATCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATg                 <  1:1763775/93‑1 (MQ=255)                                                                        aaTCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGt               >  2:5738246/1‑93 (MQ=255)                                                                         aTCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGta              >  1:4155059/1‑93 (MQ=255)                                                                         aTCGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGta              >  2:1250955/1‑93 (MQ=255)                                                                          tCGTTGGCAGAGCGGGCGTGAGGATGAAGTATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGtat             <  2:5233534/93‑1 (MQ=255)                                                                           cGTTGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGtata            <  1:2161775/93‑1 (MQ=255)                                                                             ttGGCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAAc          >  1:4594219/1‑93 (MQ=255)                                                                                gCAGAGCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACAtt       >  1:6155047/1‑93 (MQ=255)                                                                                  agagCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACAtttt     <  1:994653/93‑1 (MQ=255)                                                                                   gagCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTg    >  2:2166870/1‑93 (MQ=255)                                                                                   gagCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCAATGGTATAACATTTTg    >  2:2929819/1‑93 (MQ=255)                                                                                     gCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTa  <  2:5537464/93‑1 (MQ=255)                                                                                     gCGGGCGTGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTa  <  2:4814512/93‑1 (MQ=255)                                                                                           |                                                                                       TTTACCTGACAGTAACGGGAGAACAACAAGGAAGCATCTCCGCAGGTTGTGGGACTTCTGAGTCTACAGGTAATCGTTGGCAGAGCGGGCATGAGGATGAAATATTTACATTCTCACTCTTAAATAATATTAATAATACGGGGCTTGGTTCACAGTTCCATGGTATAACATTTTGTA  >  NC_000913/3579507‑3579683

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