Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,662,350	C→T	P599S (CCC→TCC)	ynfF →	S‑ and N‑oxide reductase, A subunit, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,662,350	0	C	T	100.0%	108.6 / NA	33	P599S (CCC→TCC)	ynfF	S‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/15);  total (18/15)

ACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAA  >  NC_000913/1662222‑1662493                                                                                                                                 |                                                                                                                                                aCTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATg                                                                                                                                >  5:126650/1‑146 (MQ=255) aCTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGaa                                                                                                                              >  2:47220/1‑148 (MQ=255)    ggATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACgag                                                                                                                                      <  8:34532/137‑1 (MQ=255)    ggATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACgag                                                                                                                                      >  7:34532/1‑137 (MQ=255)     gATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGa                                                                                                                         >  7:73898/1‑149 (MQ=255)         cTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCg                                                                                                                      <  1:47220/148‑1 (MQ=255)                 aGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAGACGACCGGGATCTTTa                                                                                                             <  6:126650/149‑1 (MQ=255)                          aCGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGa                                                                                                                   >  1:63988/1‑134 (MQ=255)                          aCGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGa                                                                                                                   <  2:63988/134‑1 (MQ=255)                                aGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGccc                                                                                                <  8:73898/147‑1 (MQ=255)                                               tCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTc                                                                                <  7:88271/148‑1 (MQ=255)                                                                      gTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTg                                                             >  6:142614/1‑144 (MQ=255)                                                                      gTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGACCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTg                                                             >  6:140157/1‑144 (MQ=255)                                                                            aTGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCAtt                                                                  <  8:19319/133‑1 (MQ=255)                                                                            aTGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCAtt                                                                  >  7:19319/1‑133 (MQ=255)                                                                                         atatCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTg                                      >  8:2110/1‑148 (MQ=255)                                                                                             ctcCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAac                                 >  1:38295/1‑149 (MQ=255)                                                                                              tcCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGt                                        <  5:61405/141‑1 (MQ=255)                                                                                              tcCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGt                                        >  6:61405/1‑141 (MQ=255)                                                                                              tcCATGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAaca                                >  6:72618/1‑149 (MQ=255)                                                                                                 aTGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCaaa                                             <  8:101237/133‑1 (MQ=255)                                                                                                 aTGCGAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCaaa                                             >  7:101237/1‑133 (MQ=255)                                                                                                     gAAAACGAAGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCgg                         >  1:120724/1‑149 (MQ=255)                                                                                                             aGGAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATc                 >  6:101664/1‑149 (MQ=255)                                                                                                              ggAACGTAACCCTGAGATGTCCGACTACGAGGAGATTAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATcc                                                      <  2:168559/111‑1 (MQ=255)                                                                                                              ggAACGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATcc                                                      >  1:168559/1‑111 (MQ=255)                                                                                                                 aCGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGACCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGaaa             <  5:142614/149‑1 (MQ=255)                                                                                                                 aCGTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGACCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGaaa             <  5:140157/149‑1 (MQ=255)                                                                                                                   gTAACCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAAtt           >  4:140587/1‑149 (MQ=255)                                                                                                                     aaCCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGccc                                                                                                <  7:30642/62‑1 (MQ=255)                                                                                                                     aaCCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGccc                                                                                                >  8:30642/1‑62 (MQ=255)                                                                                                                      aCCCTGAGATGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTAt        <  2:38295/149‑1 (MQ=255)                                                                                                                            agaTGTCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGaa  <  5:72618/149‑1 (MQ=255)                                                                                                                                 |                                                                                                                                                ACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAA  >  NC_000913/1662222‑1662493

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