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%	92.9 / NA	28	T511T (ACC→ACT)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/15);  total (13/15)

AACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGG  >  NC_000913/1540822‑1541302                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     aaCCAGCCGGTTTGGGGGCGCAGTTTTCCCTGGCCGACATAGTGTGCCCAGCCGCCCCCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCAGCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGGTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTTTCGCGGCACGCCAGTAATTTGCTCAcc                                                                                                                                                                                                                                        <  2:471786/251‑1 (MQ=255)        cGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCt                                                                                                                                                                                                                                 <  1:14424/251‑1 (MQ=255)                   cgcAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCTATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACgg                                                                                                                                                                                                                      <  1:24041/251‑1 (MQ=255)                      aGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGttt                                                                                                                                                                                                                   <  1:522369/251‑1 (MQ=255)                                   cGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACGCCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTa                                                                                                                                                                                                       >  2:501966/1‑250 (MQ=255)                                    gACATAGTGTGCCCAGCCGCCACCTCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTGCGATAGACTGGCGCGGCACGCCAGTAAGTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAAt                                                                                                                                                                                                     <  1:519768/251‑1 (MQ=255)                                        aaGTGTCCCCAGCCGCCACCGCTTTGCCCGACCCAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAAGCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCttt                                                                                                                                                                                                 <  2:533178/250‑1 (MQ=255)                                         aGTGTGCCCAGCCCCCACCGCTGTGCCCGACACAGCCGCAGAAGATCAGCATATTGATCACCCCGCGGTAGTTCATGTCCATGTGAGACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCTTCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTGTCGATAGACTTGCGCGGCACGCCAGTGATTTGCTCTCCCCAGGCTGTTGTTTATGGTTTGATCTCattgaaatcttt                                                                                                                                                                                                 <  2:88941/250‑12 (MQ=255)                                          gtgtGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCg                                                                                                                                                                                               <  1:185933/251‑1 (MQ=255)                                                    ccgccACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTc                                                                                                                                                                                     <  1:113763/251‑1 (MQ=255)                                                       ccACCGCTCTGCCCGCCAAAACCGAAAAAGATCAGCATATTGATCATCACGCGGCAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCAAGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCACAGGCTGGGGCGTACGGTTTGACTTCAGCGTAATCTTTCGCACTGTTTTCATc                                                                                                                                                                                  >  2:66582/1‑251 (MQ=255)                                                          ccGCTTTGCCCGTCACACCCGCAGAAGATCAGCATATTGATCCTCCCGCGGTAGTTCATGTCCTTGTGATACCTGTGGTTAACACCTGCGTCGAGGATAATCTTCGAGCGCCCTTGCGTTTTTTGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGAGATATGGGCGCGGTACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCTGCGTAATTTTTCGCACTGTGTTCATCTTc                                                                                                                                                                               <  2:520945/251‑1 (MQ=255)                                                          ccGCTTTGCCCGACACACCCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc                                                                                                                                                                               >  1:173126/1‑251 (MQ=255)                                                                        acaACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTTGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGa                                                                                                                                                                 <  1:44885/251‑1 (MQ=255)                                                                                                 atcataCCGCGGTAGTTCAAGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGTATAATCATCGAGCGCCCATGCGTTTTATGGTGGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATGTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATc                                                                                                                                        <  2:5958/251‑1 (MQ=255)                                                                                                                        cATGTGATACCAGTGGTTAACACCTGCGCAGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAgg                                                                                                                 >  2:521701/1‑251 (MQ=255)                                                                                                                                 ccAGTGGTTAACACCTGCGACGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCCAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGCCAGGTGTTGCCAt                                                                                                        >  2:347475/1‑251 (MQ=255)                                                                                                                                      ggTTAACACCTACGCCGAGGATACTCATCGAGCGCCCATGCGCTTTCTGGGCGGTATCGGCAAATTCACGGGCGCTGGTTTAGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATACACGCTGACCACCGGACCGGTGTTGCCCTCGACt                                                                                                   >  2:80399/1‑251 (MQ=255)                                                                                                                                       gTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACt                                                                                                   >  1:77432/1‑250 (MQ=255)                                                                                                                                                   gcCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCCGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGtttt                                                                                      >  1:362611/1‑251 (MQ=255)                                                                                                                                                                        cccATGCGTTTTATGGGAGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCAcc                                                                  >  1:440809/1‑250 (MQ=255)                                                                                                                                                                                       ggCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACCGGGTTTGCCATCGACTAACGTCAGGTTTTTAACGGGAAATTGGCGCCCCAGCACCGGATTGTGtt                                                  >  2:522155/1‑251 (MQ=255)                                                                                                                                                                                               cGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCt                                           <  2:173126/250‑1 (MQ=255)                                                                                                                                                                                                  cAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCg                                        >  1:8074/1‑250 (MQ=255)                                                                                                                                                                                                         aCGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCg                                >  1:292768/1‑251 (MQ=255)                                                                                                                                                                                                                             aCTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCaa              >  1:14620/1‑249 (MQ=255)                                                                                                                                                                                                                                   gcgGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTGTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAgggg      <  1:99149/251‑1 (MQ=255)                                                                                                                                                                                                                                       cACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAgg  <  1:214706/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     AACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGG  >  NC_000913/1540822‑1541302

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑