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%	182.3 / NA	52	E414K (GAA→AAA)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/17);  total (35/17)

GGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCT  >  NC_000913/1541224‑1541486                                                                                                                                |                                                                                                                                        ggCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTAcc                                                                                                                             >  2:723546/1‑140 (MQ=255) ggCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTAcc                                                                                                                             >  1:149689/1‑140 (MQ=255) ggCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTACAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGCTTCCGTAcc                                                                                                                             >  2:1097393/1‑140 (MQ=255)  gCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACcg                                                                                                                            >  2:1315110/1‑140 (MQ=255)     aTTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACcggcg                                                                                                                         >  1:83071/1‑140 (MQ=255)        ggCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCa                                                                                                                      <  1:1435645/140‑1 (MQ=255)        ggCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCa                                                                                                                      <  2:1252063/140‑1 (MQ=255)            cACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATgg                                                                                                                  >  1:801572/1‑140 (MQ=255)              ccAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgat                                                                                                                >  1:1475163/1‑140 (MQ=255)                aGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgattc                                                                                                              >  2:780458/1‑140 (MQ=255)                       gATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgattccagattc                                                                                                       >  2:417019/1‑140 (MQ=255)                       gATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgattccagattc                                                                                                       >  2:1328725/1‑140 (MQ=255)                           gtgtTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCAtt                                                                                                   >  1:594451/1‑140 (MQ=255)                            tgtTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCAttt                                                                                                  >  1:21531/1‑140 (MQ=255)                            tgtTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCAttt                                                                                                  >  2:1443724/1‑140 (MQ=255)                              ttttACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGc                                                                                                <  2:788460/140‑1 (MQ=255)                                   cGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCtttt                                                                                           >  2:727844/1‑140 (MQ=255)                                    gCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCttttt                                                                                          <  1:1181524/140‑1 (MQ=255)                                      tGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTtct                                                                                        >  2:203080/1‑140 (MQ=255)                                      tGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTtct                                                                                        >  1:170089/1‑140 (MQ=255)                                        cGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTcc                                                                                      >  2:490808/1‑140 (MQ=255)                                              tGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCg                                                                                <  1:70518/140‑1 (MQ=255)                                              tGCGGATTTTCAATACCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCg                                                                                <  1:959246/140‑1 (MQ=255)                                                 ggATTTTCAATGCCGCCAAAGTAGGGGAAGGCCCCGCCTGCAACAGCGACACGTTGGCCGAGCAGGGTTGACGACAATTTGGTTTCCGTTCCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGaa                                                                             <  2:99626/140‑1 (MQ=255)                                                        cAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATc                                                                      >  2:858087/1‑140 (MQ=255)                                                               gccAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGt                                                               >  2:941016/1‑140 (MQ=255)                                                                ccAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGtt                                                              <  2:391581/140‑1 (MQ=255)                                                                    aGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGc                                                          <  1:1197764/140‑1 (MQ=255)                                                                    aGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGc                                                          <  2:594451/140‑1 (MQ=255)                                                                            aaGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCaa                                                  <  1:723546/140‑1 (MQ=255)                                                                            aaGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCaa                                                  <  2:149689/140‑1 (MQ=255)                                                                             aGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAAt                                                 >  1:162039/1‑140 (MQ=255)                                                                              ggCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAAtt                                                >  1:77305/1‑140 (MQ=255)                                                                              ggCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAAtt                                                >  2:1281564/1‑140 (MQ=255)                                                                              cgCCACGCCGGCAACAGCGGCATGTTGACCGAGCAGGGTTAACGACGATTTGGTTTCCGTACCGCCGGCAATGGATTCCAGTTTCCAACTGCCTTTT‑‑TCCCCAGCGGAGACCAATCGAACCGGTCGGCACTACCAATTCa                                              >  2:672558/2‑140 (MQ=255)                                                                               gCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTc                                               >  1:558842/1‑140 (MQ=255)                                                                                ccACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCa                                              >  1:263635/1‑140 (MQ=255)                                                                                    gCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGTGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTg                                          <  2:801572/140‑1 (MQ=255)                                                                                     ccTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGc                                         <  1:1315110/140‑1 (MQ=255)                                                                                         cAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTa                                     >  2:7637/1‑140 (MQ=255)                                                                                         cAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTa                                     >  1:640027/1‑140 (MQ=255)                                                                                           aCAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTAtt                                   <  2:646996/140‑1 (MQ=255)                                                                                              gCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAAc                                >  1:166107/1‑140 (MQ=255)                                                                                               cGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACt                               >  2:88424/1‑140 (MQ=255)                                                                                               cGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACt                               >  1:1003248/1‑140 (MQ=255)                                                                                                gTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTg                              >  1:97053/1‑140 (MQ=255)                                                                                                gTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTg                              >  2:376677/1‑140 (MQ=255)                                                                                                 tCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGc                             >  1:559824/1‑140 (MQ=255)                                                                                                 tCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGc                             <  2:170089/140‑1 (MQ=255)                                                                                                         aCCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGttt                     >  1:1001540/1‑140 (MQ=255)                                                                                                             aGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCa                 <  1:726679/140‑1 (MQ=255)                                                                                                                        cGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGAttgt      >  2:743836/1‑140 (MQ=255)                                                                                                                        cGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGAttgt      <  1:1443724/140‑1 (MQ=255)                                                                                                                            aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCt  >  1:252062/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                        GGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCT  >  NC_000913/1541224‑1541486

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: ‑