Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,584,679	C→T	G603D (GGT→GAT)	ydeP ←	putative oxidoreductase YdeP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,584,679	0	C	T	100.0%	121.6 / NA	37	G603D (GGT→GAT)	ydeP	putative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/18);  total (19/18)

TCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCC  >  NC_000913/1584551‑1584807                                                                                                                                 |                                                                                                                                 tCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAatcacc                                                                                                                              <  2:41331/133‑1 (MQ=255)      aTGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGAt                                                                                                                         >  1:502103/1‑133 (MQ=255)              cTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTCATCAGGTGAAAATCACCGGGATGACGGATg                                                                                                                 >  2:128158/1‑133 (MQ=255)               ttACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATgc                                                                                                                <  1:341633/133‑1 (MQ=255)               ttACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATgc                                                                                                                <  1:386481/133‑1 (MQ=255)                tACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATgcg                                                                                                               <  2:477252/133‑1 (MQ=255)                        aaCGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTtat                                                                                                       >  2:30020/1‑133 (MQ=255)                             tGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCg                                                                                                  >  2:352566/1‑133 (MQ=255)                             tGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCg                                                                                                  >  1:274592/1‑133 (MQ=255)                              gAAGAGGGATGTTCTCACAGCCCTTTGCAGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACTGGGATGACGGATGCGCTGGTTATAGTCgg                                                                                                 <  2:173954/133‑1 (MQ=255)                                         ttctAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGgc                                                                                      >  1:494855/1‑133 (MQ=255)                                             aaCAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGgcagca                                                                                  >  1:101958/1‑133 (MQ=255)                                                  cccTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCATCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCt                                                                             >  2:469834/1‑133 (MQ=255)                                                  cccTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCt                                                                             >  2:57756/1‑133 (MQ=255)                                                  cccTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCt                                                                             >  2:23891/1‑133 (MQ=255)                                                     tttGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCa                                                                          >  2:258715/1‑133 (MQ=255)                                                           ggTAATGATATTAGCCTTCCCTTACGGCGTTATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCa                                                                    <  1:149315/133‑1 (MQ=255)                                                                   aaTTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAAt                                                            <  1:30020/133‑1 (MQ=255)                                                                      tAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCg                                                         >  1:386742/1‑133 (MQ=255)                                                                       aGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGa                                                        <  2:436379/133‑1 (MQ=255)                                                                           ttACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCAt                                                    <  1:57756/133‑1 (MQ=255)                                                                            tACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATa                                                   >  1:18748/1‑133 (MQ=255)                                                                            tACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATa                                                   >  1:267687/1‑133 (MQ=255)                                                                             aCCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACACTGGCAGCACAGCTACGATGTCATTGCGAAAGCGATCATaa                                                  >  2:226011/1‑133 (MQ=255)                                                                               cTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATc                                                <  1:258715/133‑1 (MQ=255)                                                                                 gACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCtt                                              <  1:23891/133‑1 (MQ=255)                                                                                         cATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGa                                      <  1:352566/133‑1 (MQ=255)                                                                                         cATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGa                                      <  2:502103/133‑1 (MQ=255)                                                                                                         gcagcTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCa                      <  2:274592/133‑1 (MQ=255)                                                                                                         gcagcTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCa                      <  1:128158/133‑1 (MQ=255)                                                                                                         gcagcTGCATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCa                      <  1:115537/133‑1 (MQ=255)                                                                                                              tGCATTTATCAGGTGAAAATCAACGGGATGCCGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGctc                 <  1:168287/133‑1 (MQ=255)                                                                                                                cATTTATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTg               >  1:12522/1‑133 (MQ=255)                                                                                                                  tttATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCa                      <  2:103343/124‑1 (MQ=255)                                                                                                                  tttATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCa                      >  1:103343/1‑124 (MQ=255)                                                                                                                  tttATCAGGTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTggg             <  1:503410/133‑1 (MQ=255)                                                                                                                          gTGAAAATCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGtgct     >  1:41798/1‑133 (MQ=255)                                                                                                                             aaaaTCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGctgcc  >  2:281594/1‑133 (MQ=255)                                                                                                                                 |                                                                                                                                 TCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCGTCATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCC  >  NC_000913/1584551‑1584807

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