Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,139,164	Δ1 bp	coding (438/867 nt)	yghU →	disulfide reductase/organic hydroperoxide reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,139,163	0	G	.	100.0%	176.8 / NA	42	coding (437/867 nt)	yghU	disulfide reductase/organic hydroperoxide reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (18/24);  total (18/24)

CCTGCTTTATCTGGCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACTGGCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCATTACGCACCGGTA  >  NC_000913/3139089‑3139236                                                                           |                                                                          ccTGCTTTATCTGGCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTc                                                                    >  2:580354/1‑81 (MQ=255)         aTCTGGCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTa                                                              >  1:580370/1‑79 (MQ=255)         aTCTGGCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTa                                                              >  1:580358/1‑79 (MQ=255)           cTGGCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAgg                                                          >  2:580356/1‑81 (MQ=255)              gCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAgg                                                          <  1:580301/78‑1 (MQ=255)              gCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAgg                                                          <  2:580302/78‑1 (MQ=255)              gCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGgcg                                                       <  2:580300/81‑1 (MQ=255)              gCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGgcg                                                       >  1:580366/1‑81 (MQ=255)              gCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGgc                                                        <  1:580299/80‑1 (MQ=255)                ggAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGcgcg                                                     >  2:580369/1‑81 (MQ=255)                    aaaTTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCAc                                                 <  1:580303/81‑1 (MQ=255)                        ttGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGtt                                             >  2:580400/1‑81 (MQ=255)                        ttGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGtt                                             >  1:580373/1‑81 (MQ=255)                        ttGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGtt                                             >  1:580352/1‑81 (MQ=255)                         tGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGtt                                             <  1:580304/80‑1 (MQ=255)                         tGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTc                                            >  1:580365/1‑81 (MQ=255)                          ggCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTcc                                           >  1:580351/1‑81 (MQ=255)                           gCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTcc                                           <  2:580305/80‑1 (MQ=255)                            cTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTc                                         <  1:580306/81‑1 (MQ=255)                            cTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTc                                         <  2:580307/81‑1 (MQ=255)                               cTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTcggc                                      >  1:580355/1‑81 (MQ=255)                               cTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTcg                                        <  1:580308/79‑1 (MQ=255)                                   cTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCggt                                   <  1:580309/80‑1 (MQ=255)                                   cTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCggt                                   <  2:580310/80‑1 (MQ=255)                                      ccGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGtt                               <  2:580311/81‑1 (MQ=255)                                       cGCAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGttt                              <  1:580312/81‑1 (MQ=255)                                         cAGGATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTg                            <  2:580313/81‑1 (MQ=255)                                            gATTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTc                         <  2:580315/81‑1 (MQ=255)                                            gATTTGGCAAAGCGTACTGAAACGATGAACT‑CCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTc                         <  2:580314/81‑1 (MQ=255)                                             aTTTGGCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCa                        >  1:580357/1‑81 (MQ=255)                                                 ggCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACttt                    >  1:580377/1‑81 (MQ=255)                                                  gCAAAGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACttt                    >  2:580376/1‑80 (MQ=255)                                                     aaGCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTAcc                <  2:580316/81‑1 (MQ=255)                                                       gCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCAt              <  1:580317/81‑1 (MQ=255)                                                       gCGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCAt              >  1:580363/1‑81 (MQ=255)                                                        cGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCAtt             <  2:580319/81‑1 (MQ=255)                                                        cGTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCAtt             <  2:580318/81‑1 (MQ=255)                                                         gTACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCATTa            >  2:580375/1‑81 (MQ=255)                                                          tACTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCATTAc           <  1:580320/81‑1 (MQ=255)                                                            cTGAAACGATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCATTACg          <  1:580321/80‑1 (MQ=255)                                                                   gATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCATTACGCACCGGt   >  2:580342/1‑80 (MQ=255)                                                                   gATGAACT‑GCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCATTACGCACCGGTa  <  2:580322/81‑1 (MQ=255)                                                                           |                                                                          CCTGCTTTATCTGGCGGAGAAATTTGGCTACTTCCTGCCGCAGGATTTGGCAAAGCGTACTGAAACGATGAACTGGCTGTTCTGGTTACAGGGCGCGGCACCGTTCCTCGGCGGTGGTTTTGGTCACTTTTACCATTACGCACCGGTA  >  NC_000913/3139089‑3139236

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