Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,302,066	G→A	R230C (CGT→TGT)	napA ←	nitrate reductase, periplasmic, large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,302,066	0	G	A	100.0%	111.0 / NA	35	R230C (CGT→TGT)	napA	nitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/18);  total (17/18)

AGTTGGCGATGTAGTTCAGGATCACCAGGTCAGATTGCGGCGTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAA  >  NC_000913/2301934‑2302195                                                                                                                                     |                                                                                                                                  aGTTGGCGATGTAGTTCAGGATCACCAGGTCAGATTGCGGCGTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGt                                                                                                                             <  1:727524/139‑1 (MQ=255)   ttGGCGATGTAGTTCAGGATCACCAGGTCAGATTGCGGCGTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGa                                                                                                                           <  2:635954/139‑1 (MQ=255)     ggCGATGTAGTTCAGGATCACCAGGTCAGATTGCGGCGTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATg                                                                                                                         >  2:122061/1‑139 (MQ=255)           gTAGTTCAGGATCACCAGGTCAGATTGCGGCGTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGAc                                                                                                                   >  2:444633/1‑139 (MQ=255)                 cAGGATCACCAGGTCAGATTGCGGCGTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAgg                                                                                                             <  1:444633/139‑1 (MQ=255)                                      cggcgTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATg                                                                                                                         <  2:462432/106‑1 (MQ=255)                                      cggcgTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATg                                                                                                                         >  1:462432/1‑106 (MQ=255)                                       ggcgTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATgc                                                                                                                        >  1:736445/1‑106 (MQ=255)                                       ggcgTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATgc                                                                                                                        <  2:736445/106‑1 (MQ=255)                                                 tgatgCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCca                                                                             <  1:335162/139‑1 (MQ=255)                                                   atgCCGTTATTCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCcaca                                                                           >  1:407407/1‑139 (MQ=255)                                                       cGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGcac                                                                       >  2:533414/1‑139 (MQ=255)                                                                            tACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGAt                                                  <  1:122061/139‑1 (MQ=255)                                                                                     ggtaggtaGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACaa                                                                     >  1:490506/1‑111 (MQ=255)                                                                                     ggtaggtaGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACaa                                                                     <  2:490506/111‑1 (MQ=255)                                                                                      gtaggtagAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGAt                                                  <  2:554937/129‑1 (MQ=255)                                                                                      gtaggtagAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGAt                                                  >  1:554937/1‑129 (MQ=255)                                                                                           tagAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATgc                                                                                                                        <  1:324273/54‑1 (MQ=255)                                                                                           tagAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATgc                                                                                                                        >  2:324273/1‑54 (MQ=255)                                                                                             gAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCa                                                                                                   <  1:242477/73‑1 (MQ=255)                                                                                             gAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCa                                                                                                   >  2:242477/1‑73 (MQ=255)                                                                                             gAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCa                                 <  1:533414/139‑1 (MQ=255)                                                                                                aGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATgcgc                                                                                                                      >  1:360888/1‑51 (MQ=255)                                                                                                aGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATgcgc                                                                                                                      <  2:360888/51‑1 (MQ=255)                                                                                                aGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGc                                                                         >  2:285217/1‑96 (MQ=255)                                                                                                aGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGc                                                                         <  1:285217/96‑1 (MQ=255)                                                                                                 gCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCgg                             <  2:407407/139‑1 (MQ=255)                                                                                                       ccACGGTGACGTTCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATc                       >  2:178738/1‑139 (MQ=255)                                                                                                                  ttCTGGTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCa                                                                                                   >  2:510303/1‑52 (MQ=255)                                                                                                                  ttCTGGTTAGACAGACGACAGTTAGTGATGCGCGACCTGAGGATCGGGTGCa                                                                                                   <  1:510303/52‑1 (MQ=255)                                                                                                                       gTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCa                                 <  2:286233/113‑1 (MQ=255)                                                                                                                       gTTAGAGAGACGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCa                                 >  1:286233/1‑113 (MQ=255)                                                                                                                            agagaCGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCgg                             >  2:420011/1‑112 (MQ=255)                                                                                                                            agagaCGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCgg                             <  1:420011/112‑1 (MQ=255)                                                                                                                            agagaCGACAGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATaaa  >  1:726147/1‑139 (MQ=255)                                                                                                                                 cgacaGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTc                          >  1:597359/1‑110 (MQ=255)                                                                                                                                 cgacaGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTc                          <  2:597359/110‑1 (MQ=255)                                                                                                                                     |                                                                                                                                  AGTTGGCGATGTAGTTCAGGATCACCAGGTCAGATTGCGGCGTAAAGATGATGCCGTTATCCGCCAGCTCGAAGCTACGATGCTGGTAGGTAGAAAGCACCGCCACGGTGACGTTCTGGTTAGAGAGACGACGGTTAGTGATGCGCGACCAGAGGATCGGGTGCATTTCCGCCATGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAA  >  NC_000913/2301934‑2302195

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