Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,541,351	C→T	E414K (GAA→AAA)	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,541,351	0	C	T	100.0%	126.0 / NA	36	E414K (GAA→AAA)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/17);  total (19/17)

GTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTT  >  NC_000913/1541214‑1541487                                                                                                                                          |                                                                                                                                         gTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTgg                                                                                                                                        >  2:647560/1‑140 (MQ=255)                  ggCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCa                                                                                                                       <  2:612621/140‑1 (MQ=255)                        ccAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgat                                                                                                                 >  2:684287/1‑140 (MQ=255)                         cAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgatt                                                                                                                <  1:345341/140‑1 (MQ=255)                                         tttACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGcc                                                                                                <  2:262600/140‑1 (MQ=255)                                             cGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGGTTCCATTTGCCtttt                                                                                            >  2:437974/1‑140 (MQ=255)                                                tGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGcc                                                                                                <  2:287730/133‑1 (MQ=255)                                                tGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGcc                                                                                                >  1:287730/1‑133 (MQ=255)                                                tGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTtct                                                                                         >  2:918679/1‑140 (MQ=255)                                                  cGAAAATGCGGATTTTCAATGCCGCCGAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTcc                                                                                       >  2:885734/1‑140 (MQ=255)                                                   gAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTccc                                                                                      <  2:858644/140‑1 (MQ=255)                                                            gATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGaa                                                                              <  1:638733/139‑1 (MQ=255)                                                            gATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGaa                                                                              <  1:893974/139‑1 (MQ=255)                                                            gATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGaa                                                                              >  2:638733/1‑139 (MQ=255)                                                            gATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGaa                                                                              >  2:893974/1‑139 (MQ=255)                                                                         gccAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGt                                                                >  1:671/1‑140 (MQ=255)                                                                          ccAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGtt                                                               <  1:304046/140‑1 (MQ=255)                                                                                    ggAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTAcc                                                     >  1:259247/1‑140 (MQ=255)                                                                                        ggCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAAtt                                                 >  1:975734/1‑140 (MQ=255)                                                                                        ggCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAAtt                                                 >  2:941312/1‑140 (MQ=255)                                                                                         gCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTc                                                >  1:253711/1‑140 (MQ=255)                                                                                         gCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTc                                                >  1:58403/1‑140 (MQ=255)                                                                                            aCGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCAcc                                             >  1:493404/1‑140 (MQ=255)                                                                                              gCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTg                                           <  1:379686/140‑1 (MQ=255)                                                                                                cTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCg                                         <  1:1056596/140‑1 (MQ=255)                                                                                                cTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCg                                         <  2:590654/140‑1 (MQ=255)                                                                                                cTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCg                                         <  1:501619/140‑1 (MQ=255)                                                                                                      cAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTa                                   <  1:941312/140‑1 (MQ=255)                                                                                                       aGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTaa                                  <  2:493404/140‑1 (MQ=255)                                                                                                                   aCCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGttt                      >  2:73418/1‑140 (MQ=255)                                                                                                                      gAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTcc                   >  1:883906/1‑140 (MQ=255)                                                                                                                        gCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCAc                 >  1:852102/1‑140 (MQ=255)                                                                                                                                 aCGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGAttg        >  1:735980/1‑140 (MQ=255)                                                                                                                                      aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCt   <  2:852102/140‑1 (MQ=255)                                                                                                                                      aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCt   <  2:883906/140‑1 (MQ=255)                                                                                                                                      aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCt   >  2:492262/1‑140 (MQ=255)                                                                                                                                       aTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCtt  <  1:647560/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         GTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTT  >  NC_000913/1541214‑1541487

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