Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,539,530	A→G	S1021P (TCC→CCC)	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,539,530	0	A	G	100.0%	108.6 / NA	31	S1021P (TCC→CCC)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/14);  total (17/14)

CTGATAGAGCTGCTGGCGTCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGC  >  NC_000913/1539453‑1539607                                                                              |                                                                              cTGATAGAGCTGCTGGCGTCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGAc                                                                            >  2:291232/1‑81 (MQ=255)    aTAGAGCTGCTGGCGTCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGt                                                                          >  2:291241/1‑80 (MQ=255)    aTAGAGCTGCTGGCGTCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTg                                                                         <  1:291200/81‑1 (MQ=255)       gagCTGCTGGCGTCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGAt                                                                       <  2:291201/80‑1 (MQ=255)                 cGTCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGa                                                            >  1:291231/1‑81 (MQ=255)                  gTCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGAc                                                           >  2:291247/1‑81 (MQ=255)                   tCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGAcc                                                          <  2:291203/81‑1 (MQ=255)                   tCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGAcc                                                          <  2:291202/81‑1 (MQ=255)                       gACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGAc                                                      <  1:291204/81‑1 (MQ=255)                       gACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGAc                                                      >  2:291249/1‑81 (MQ=255)                       gACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGAc                                                      >  2:291235/1‑81 (MQ=255)                        aCAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGAc                                                      >  1:291238/1‑80 (MQ=255)                         cAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCa                                                    <  1:291206/81‑1 (MQ=255)                         cAGCGTGCGCCACGGAATTAACTCATGAACGATGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCa                                                    <  2:291205/81‑1 (MQ=255)                          aGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAg                                                   >  2:291250/1‑81 (MQ=255)                          aGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAg                                                   <  1:291207/81‑1 (MQ=255)                                    aCGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTgg                                         <  2:291208/81‑1 (MQ=255)                                    aCGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTg                                          >  1:291257/1‑80 (MQ=255)                                                 aTGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTa                              <  2:291209/79‑1 (MQ=255)                                                 aTGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTAcg                            >  1:291253/1‑81 (MQ=255)                                                 aTGAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTAcg                            >  2:291239/1‑81 (MQ=255)                                                   gAACGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTAcgcg                          >  1:291242/1‑81 (MQ=255)                                                    aaCGTTGGTGTATCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCgg                         <  2:291210/81‑1 (MQ=255)                                                               aTCCCGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGaa              >  1:291244/1‑81 (MQ=255)                                                                 cccGCGTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAata            >  1:291254/1‑81 (MQ=255)                                                                   cgcgTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATc          <  1:291211/81‑1 (MQ=255)                                                                   cgcgTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATc          <  1:747107/81‑1 (MQ=255)                                                                    gcgTTATAGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATcg         >  2:291245/1‑81 (MQ=255)                                                                         ataGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGaa     <  1:291212/80‑1 (MQ=255)                                                                           aGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGc  >  1:291240/1‑81 (MQ=255)                                                                           aGGGGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGc  >  2:291265/1‑81 (MQ=255)                                                                              |                                                                              CTGATAGAGCTGCTGGCGTCCCGACAGCGTGCGCCACGGAATTAACTCATGAACGTTGGTGTATCCCGCGTTATAGGAGACGTGATCGCTTTCCAGACCGGACCAGGTGGGGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGC  >  NC_000913/1539453‑1539607

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