Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,538,930	A→G	F1221L (TTC→CTC)	narZ ←	nitrate reductase Z subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,538,930	0	A	G	100.0%	132.6 / NA	38	F1221L (TTC→CTC)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/14);  total (24/14)

GATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCC  >  NC_000913/1538803‑1539059                                                                                                                                |                                                                                                                                  gATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGtt                                                                                                                       <  2:291838/140‑1 (MQ=255) gATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGtt                                                                                                                       <  2:727784/140‑1 (MQ=255) gATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGtt                                                                                                                       <  2:453393/140‑1 (MQ=255)               ggTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGc                                                                                                         <  2:913163/140‑1 (MQ=255)                gTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGcc                                                                                                        >  2:31948/1‑140 (MQ=255)                     gCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCgtagt                                                                                                   >  2:286503/1‑140 (MQ=255)                     gCACCATGCCGACTTGTGAACGGATTTTCA‑TTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCgtagta                                                                                                  >  1:555539/1‑140 (MQ=255)                         cATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATATTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGtt                                                                                               >  2:526229/1‑140 (MQ=255)                           tGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCgtag                                                                                                    <  1:396284/133‑1 (MQ=255)                           tGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCgtag                                                                                                    >  2:396284/1‑133 (MQ=255)                           tGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTaa                                                                                             >  1:283074/1‑140 (MQ=255)                           tGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTaa                                                                                             <  1:490405/140‑1 (MQ=255)                              cGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGc                                                                                          >  1:54306/1‑140 (MQ=255)                                       aaCGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGccc                                                                                        <  2:653202/133‑1 (MQ=255)                                       aaCGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGccc                                                                                        >  1:653202/1‑133 (MQ=255)                                                      tttCGCCTCCTGTACCTGATCGCGCCCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGc                                                                  <  2:147521/140‑1 (MQ=255)                                                      tttCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGACCAGCTGCGCGTAACCGc                                                                  >  1:355429/1‑140 (MQ=255)                                                      tttCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGACCAGCTGCGCGTAACCGc                                                                  >  1:644377/1‑140 (MQ=255)                                                       ttCGCCTCCTGTACCTGATCGCGACCTTCATCACCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGcc                                                                 <  1:772660/140‑1 (MQ=255)                                                                   cccTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGt                                                     <  2:102981/139‑1 (MQ=255)                                                                        aTCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGtt                                                >  1:874510/1‑140 (MQ=255)                                                                         tCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGttt                                               >  2:205227/1‑140 (MQ=255)                                                                         tCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGttt                                               >  1:749238/1‑140 (MQ=255)                                                                          cgcgACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTc                                              >  1:905838/1‑140 (MQ=255)                                                                            cgACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCgg                                            >  1:639180/1‑140 (MQ=255)                                                                             gACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCggg                                           >  1:806143/1‑140 (MQ=255)                                                                               ccTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCa                                         >  1:439968/1‑140 (MQ=255)                                                                                  tcttgatcCCGCCATTTACCGTTCTTGATCTTGGGGATAATGATGCGCCCATAGGGGTTCGTCCCGAGGGTGCGGTAGTAGTAAAGCCCCCAGCCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAAc                                      <  2:639180/135‑1 (MQ=255)                                                                                  tcatcaTCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAAc                                      >  2:86196/1‑140 (MQ=255)                                                                                           agccagTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAAc                             <  2:439968/140‑1 (MQ=255)                                                                                            gccagTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAAcc                            <  2:874510/140‑1 (MQ=255)                                                                                             ccagTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCggg                                 <  2:527325/134‑1 (MQ=255)                                                                                             ccagTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCggg                                 >  1:527325/1‑134 (MQ=255)                                                                                             ccagTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCg                           >  1:622909/1‑140 (MQ=255)                                                                                             ccagTTAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCg                           >  1:52084/1‑140 (MQ=255)                                                                                                 ttAACGTTCTTCATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGtt                       <  1:5904/140‑1 (MQ=255)                                                                                                          ttcATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATgcc              <  2:283074/140‑1 (MQ=255)                                                                                                           tcATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATgccg             >  2:450402/1‑140 (MQ=255)                                                                                                           tcATCTTGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATgccg             <  2:749238/140‑1 (MQ=255)                                                                                                                 tGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCgcgc       >  2:473908/1‑140 (MQ=255)                                                                                                                 tGCGGATCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGACCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCgcgc       <  2:355429/140‑1 (MQ=255)                                                                                                                 ggCGGACCATGAGGAGCTGATGGCGGTTGGATCAGCCGGTCCCGTAGTAGTTAAAGCCCAGGACCAGTTGCCGGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCgcgc       <  2:644377/139‑1 (MQ=255)                                                                                                                      aTCATGATGAGCTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGcc  >  2:434251/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                  GATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCC  >  NC_000913/1538803‑1539059

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