Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,994,218	T→C	S123S (AGT→AGC)	yzcX →	protein YzcX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,994,218	0	T	C	100.0%	129.3 / NA	37	S123S (AGT→AGC)	yzcX	protein YzcX
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/22);  total (15/22)

TTCCGGCACGATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCACGCTGACGGAAC  >  NC_000913/3994141‑3994294                                                                              |                                                                             ttCCGGCACGATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTac                                                                           <  1:732800/81‑1 (MQ=255)    cGGCACGATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTacac                                                                         <  1:732801/80‑1 (MQ=255)    cGGCACGATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTAcacc                                                                        <  2:732801/81‑1 (MQ=255)     ggCACGATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTAcacca                                                                       <  2:732802/81‑1 (MQ=255)        aCGATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGt                                                                    >  2:732859/1‑81 (MQ=255)          gATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCt                                                                  <  1:732804/81‑1 (MQ=255)          gATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCt                                                                  >  2:732852/1‑81 (MQ=255)              cATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTcc                                                              >  2:732853/1‑81 (MQ=255)               aTCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTccc                                                             <  1:732805/81‑1 (MQ=255)                tCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCt                                                            <  1:732806/81‑1 (MQ=255)                  cGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGt                                                          <  1:732807/81‑1 (MQ=255)                   gTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGtt                                                         <  2:732807/81‑1 (MQ=255)                       ccGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTg                                                     >  1:732845/1‑81 (MQ=255)                          gAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTa                                                  <  2:732808/81‑1 (MQ=255)                             taatGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTaaa                                                >  2:732856/1‑80 (MQ=255)                             taatGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTaaa                                                >  1:732862/1‑80 (MQ=255)                             taatGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAAc                                               <  1:732810/81‑1 (MQ=255)                                tgatgTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACAtt                                            <  1:732811/81‑1 (MQ=255)                                 gatgTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTa                                           >  2:732854/1‑81 (MQ=255)                                 gatgTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTa                                           >  2:732829/1‑81 (MQ=255)                                      ctctGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGAt                                       >  1:732867/1‑80 (MQ=255)                                         tgtgTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATc                                      <  1:732812/78‑1 (MQ=255)                                           tgtAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGtt                                 <  2:732812/81‑1 (MQ=255)                                               gCGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTAc                              <  1:732814/80‑1 (MQ=255)                                                cGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTAcc                             >  1:732852/1‑80 (MQ=255)                                                cGAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACcg                            <  1:732815/81‑1 (MQ=255)                                                 gAAAGATTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACcga                           <  2:732815/81‑1 (MQ=255)                                                      aTTTGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACcgagcgag                      <  2:732816/81‑1 (MQ=255)                                                       tttGTCTCTTCATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACcgagcgag                      >  1:732842/1‑80 (MQ=255)                                                                 cATTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCAc             <  2:732817/79‑1 (MQ=255)                                                                  aTTAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCACGCt          <  2:732818/81‑1 (MQ=255)                                                                   ttAGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCACGCTg         <  2:732819/81‑1 (MQ=255)                                                                     aGGGCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCACGCTGAc       >  2:732833/1‑81 (MQ=255)                                                                      gggCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCACGCTGAc       <  1:732821/80‑1 (MQ=255)                                                                       ggCGCAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCACGCTTACg      >  1:732891/1‑80 (MQ=255)                                                                        gcgcAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCACGCTGACgg     >  2:732863/1‑80 (MQ=255)                                                                          gcAGCTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCACGCTGACGGAAc  >  1:732865/1‑81 (MQ=255)                                                                              |                                                                             TTCCGGCACGATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCATCCTCACGCTGACGGAAC  >  NC_000913/3994141‑3994294

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