Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,541,058	G→A	T511T (ACC→ACT)	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,058	0	G	A	100.0%	80.9 / NA	24	T511T (ACC→ACT)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/13);  total (11/13)

TGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACG  >  NC_000913/1540780‑1541344                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                               ttgcgcGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCt                                                                                                                                                                                                                                                                           <  2:244537/298‑1 (MQ=255) tGGCGCGGTGGTCGGTTCCGGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCt                                                                                                                                                                                                                                                                           <  2:734422/300‑1 (MQ=255) tGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCt                                                                                                                                                                                                                                                                           <  1:88496/300‑1 (MQ=255)                                 ggCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCttt                                                                                                                                                                                                                                           >  1:81354/1‑300 (MQ=255)                                                                 gTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGa                                                                                                                                                                                                           >  2:350060/1‑300 (MQ=255)                                                                                  tAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCAGTAATTCGCCAGt                                                                                                                                                                                          >  2:460476/1‑300 (MQ=255)                                                                                                       cTTTGCCAGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGTCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGacc                                                                                                                                                                     <  2:314946/300‑1 (MQ=255)                                                                                                                   caACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCCTGTCCATGTGATACCATTGGTTAACACCTGCGCCTAGGAGAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGTGCGATGGTGTCGATATACTGGCGCGGCACGCCAGTAATTTGCCCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtg                                                                                                                                                         <  2:218893/300‑1 (MQ=255)                                                                                                                   caACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtg                                                                                                                                                         <  1:183218/300‑1 (MQ=255)                                                                                                                        gCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGcc                                                                                                                                                    <  1:558147/300‑1 (MQ=255)                                                                                                                                                      gTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGAATAACGTCAGGTTTTTAACGGGCaa                                                                                                                      >  2:572377/1‑300 (MQ=255)                                                                                                                                                         gTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCCCCGGACAGGTGTTGCCATCGACTACCGTCCGGTTTTTACCGGGCAATTg                                                                                                                   >  2:94754/1‑300 (MQ=255)                                                                                                                                                               gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCAc                                                                                                             >  1:642056/1‑300 (MQ=255)                                                                                                                                                                   aTGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCAATCGACTAACGTCAGGTTTTTAAAGGGCAATTGGCGCACCAGc                                                                                                         >  2:727147/1‑300 (MQ=255)                                                                                                                                                                    tgtgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCAACAGCa                                                                                                        >  1:314582/1‑300 (MQ=255)                                                                                                                                                                       gATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCg                                                                                                     >  2:459762/1‑300 (MQ=255)                                                                                                                                                                          aCCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGAt                                                                                                  <  1:473314/300‑1 (MQ=255)                                                                                                                                                                                      cacCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGc                                                                                      >  1:99472/1‑300 (MQ=255)                                                                                                                                                                                                                  cccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCANCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATgccgcc                                                          >  2:655284/1‑300 (MQ=255)                                                                                                                                                                                                                           tttATGGGTGGTATTGGTAGATTCACGGTCGATGGTTTAGATGTACTGGTACGACACGCCAGTAATTTGTTCACCCCAGGCTGGGGTGTACGGTTTGATTGCAGCGTAATCTTTCGCACTGTTTTCATCGTCCGGCCCGCGGTCGAGGCCGTAATTCGCCAGTACCAGATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAGTGCGGATTTTCAATGTCGCCAGAGTgggg                                                 <  2:91899/300‑1 (MQ=255)                                                                                                                                                                                                                            ttATGGGCGGTATAGGCAAATTCACGGTCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGACAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAgggg                                                <  1:367292/300‑1 (MQ=255)                                                                                                                                                                                                                                               aTTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAg                             <  1:559997/300‑1 (MQ=255)                                                                                                                                                                                                                                                      ggCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATg                      <  1:42308/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          ggCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGGGTACGGTTTGATTTCAGCGTAATCTTTCGCACGGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACg  <  2:267440/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                               TGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACG  >  NC_000913/1540780‑1541344

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