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

GTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCT  >  NC_000913/1541214‑1541486                                                                                                                                          |                                                                                                                                        gTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTgg                                                                                                                                       >  2:27008/1‑140 (MQ=255)     ttAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTc                                                                                                                                   >  2:1083612/1‑140 (MQ=255)           ggCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTAcc                                                                                                                             >  2:685104/1‑140 (MQ=255)             cAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACcgg                                                                                                                           <  2:228877/140‑1 (MQ=255)                ttGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACcggcgg                                                                                                                        <  1:958042/140‑1 (MQ=255)                   gcgcACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCaa                                                                                                                     <  2:698852/140‑1 (MQ=255)                                 gATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATATGGTTTCCGTACCGGCGGCAATGgattccagat                                                                                                         <  2:673454/138‑1 (MQ=255)                                 gATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATATGGTTTCCGTACCGGCGGCAATGgattccagat                                                                                                         >  1:673454/1‑138 (MQ=255)                                    tgtgTTTTACGCTACGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCAt                                                                                                    <  2:1063233/140‑1 (MQ=255)                                         tttACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGcc                                                                                               <  1:685104/140‑1 (MQ=255)                                                tGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTtct                                                                                        >  2:419486/1‑140 (MQ=255)                                                  cGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTcc                                                                                      >  1:381936/1‑140 (MQ=255)                                                   gAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTccc                                                                                     <  1:1083612/140‑1 (MQ=255)                                                                          ccAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGtt                                                              <  1:607816/140‑1 (MQ=255)                                                                          ccAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATCTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGtt                                                              >  2:909065/1‑140 (MQ=255)                                                                          ccAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATCTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGtt                                                              <  1:909065/140‑1 (MQ=255)                                                                              aGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGc                                                          <  2:769134/140‑1 (MQ=255)                                                                                         gCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTc                                               >  2:732559/1‑140 (MQ=255)                                                                                              gCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTg                                          <  1:569065/140‑1 (MQ=255)                                                                                               ccTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGc                                         <  1:419486/140‑1 (MQ=255)                                                                                                   cAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTa                                     >  1:797255/1‑140 (MQ=255)                                                                                                   cAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTa                                     >  1:528472/1‑140 (MQ=255)                                                                                                         cGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACt                               >  1:1013519/1‑140 (MQ=255)                                                                                                         cGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACt                               >  1:1069739/1‑140 (MQ=255)                                                                                                         cGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACt                               >  1:322924/1‑140 (MQ=255)                                                                                                             aTGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTa                           >  2:260439/1‑140 (MQ=255)                                                                                                             aTGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTa                           <  2:797255/140‑1 (MQ=255)                                                                                                                ttGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACgg                        <  1:601141/140‑1 (MQ=255)                                                                                                                     cGAGCATGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTc                   >  2:654735/1‑140 (MQ=255)                                                                                                                      gAGCATGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTcc                  <  1:654735/140‑1 (MQ=255)                                                                                                                       aGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCa                 >  2:107131/1‑140 (MQ=255)                                                                                                                        gCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCAc                >  1:816096/1‑140 (MQ=255)                                                                                                                           gggTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGc             <  2:528472/140‑1 (MQ=255)                                                                                                                            ggTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCg            <  2:322924/140‑1 (MQ=255)                                                                                                                                  cGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGAttgt      <  1:732559/140‑1 (MQ=255)                                                                                                                                      aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCGCCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCt  >  1:154378/1‑140 (MQ=255)                                                                                                                                      aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTc                   >  2:804516/1‑123 (MQ=255)                                                                                                                                      aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTc                   <  1:804516/123‑1 (MQ=255)                                                                                                                                      aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCt  >  2:723472/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        GTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCT  >  NC_000913/1541214‑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: ‑