Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,540,626	T→C	L655L (CTA→CTG)	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,540,626	0	T	C	100.0%	58.6 / NA	17	L655L (CTA→CTG)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/10);  total (7/10)

ATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTG  >  NC_000913/1540491‑1540755                                                                                                                                        |                                                                                                                                  aTCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCagcc                                                                                                                               >  2:248425/1‑140 (MQ=255)   cTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCagccag                                                                                                                             <  1:218457/140‑1 (MQ=255)         gttCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCagc                                                                                                                                >  1:228409/1‑131 (MQ=255)         gttCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCagc                                                                                                                                <  2:228409/131‑1 (MQ=255)                               gCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAaat                                                                                                 >  2:279550/1‑140 (MQ=255)                                cccGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatc                                                                                                >  1:260976/1‑140 (MQ=255)                                  cGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaa                                                                                              >  2:344019/1‑140 (MQ=255)                                          gcgcCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGa                                                                                      <  1:358447/140‑1 (MQ=255)                                          gcgcCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGa                                                                                      <  1:344019/140‑1 (MQ=255)                                           cgcCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGAc                                                                                     <  1:388400/140‑1 (MQ=255)                                           cgcCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGAc                                                                                     <  1:279550/140‑1 (MQ=255)                                             ccTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCg                                                                                   <  1:248425/140‑1 (MQ=255)                                                           tCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGc                                                                     >  2:23412/1‑140 (MQ=255)                                                                                            gCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCaga                                    <  1:149895/140‑1 (MQ=255)                                                                                             cTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCagag                                   <  1:23412/140‑1 (MQ=255)                                                                                                       gAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCAt                         >  2:307683/1‑140 (MQ=255)                                                                                                        aGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATa                        <  1:307683/140‑1 (MQ=255)                                                                                                                              gcgcAGAAGGCAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTg  >  2:373137/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                  ATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTG  >  NC_000913/1540491‑1540755

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