Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,539,947	C→G	A882P (GCG→CCG)	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,947	0	C	G	100.0%	95.4 / NA	28	A882P (GCG→CCG)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (21/7);  total (21/7)

TTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGATGCATATCC  >  NC_000913/1539703‑1540191                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     ttGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCCATGGCGGGGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGGGGAATTGAGTTTGCCGAGGAAATCGATTTCCTCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCCATGGCCCGAGCGGGGCAAAGCGTTCCTACGTAGAAGGGTAGTCGCGCTCCCACGCCACAAAATTCGGCGGGGtttt                                                                                                                                                                                                                                                >  2:467875/1‑251 (MQ=255)                      agAATCACTTCTGAAGCGTCCATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGGTTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGGGGTAAAGCGTTCATAAGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGGCTGGAATCCGATCGCAttt                                                                                                                                                                                                                           >  1:175875/1‑250 (MQ=255)                                                                                             gtAATTGAGTTTGCCGAGGCAAACGATTTCATCCTGGGTCTTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCACCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGATGTGCGAGCTCTGCCGGAGCGTCATGCCGTAGGGGTTGt                                                                                                                                                     >  2:315164/1‑250 (MQ=255)                                                                                                         gCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCACTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTCGGGGGTGTAATACCACGTCgg                                                                                                                                        >  2:428659/1‑251 (MQ=255)                                                                                                               gAAATCGATTTCATCCTGGGAATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCtt                                                                                                                                  >  2:247215/1‑251 (MQ=255)                                                                                                                 aaTCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCCGTAGGGGTTGTAATACCACGTCGGTTTCTTTg                                                                                                                                >  2:359914/1‑251 (MQ=255)                                                                                                                         tCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATg                                                                                                                        >  1:445028/1‑251 (MQ=255)                                                                                                                         tCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGGCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGCGTCATGCAGGAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATg                                                                                                                        >  1:1852/1‑251 (MQ=255)                                                                                                                                         aCGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTAGCGCTCCACCGCCACAATATTCGGCGGGGGTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGACGAGCACTTCACACGGCTGTGCGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGcaaccttgc                                                                                                        >  1:24563/1‑244 (MQ=255)                                                                                                                                              cccctttccccGCCATGGCCAATTTTGTCCATCAAGGGCCGGAGAGAGGAAAAGCGTTAATGCGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGaaaa                                                                                                   <  1:359914/242‑1 (MQ=255)                                                                                                                                                ccttttccccccATTGCAAGGTTTCTCCATCAGGGGCCCGAGCGGGTTAAAGCGTTCATCCGTAGCGGGGTAGTCGCCCCCCCCCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATg                                                                                                 <  2:122457/241‑1 (MQ=255)                                                                                                                                                  cTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGCATCAGAt                                                                                                                                                                                                                                  >  2:495599/1‑120 (MQ=255)                                                                                                                                                         gCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCAACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGCTCGCATTAGCCTTTGCGCCAGTCGAACACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGCTGACCCACGCCCACTTGCGAAACTGCTTTGGCgg                                                                                        >  1:495735/1‑250 (MQ=255)                                                                                                                                                              tGCCAAGATTGACCACCCATGGCACGAGCGAGGTACAGCGATCATACGTAGCCGGGGAGTCCCGCTCCCCCGCCACAATATTCGGAGGGGTCCTGAGTGGACTCAGATCGCATTCGCCTATGCAACAGCCGAACAATTCACACGGCTGTGAGAACTCTGACGGAGAGTCATGCAGTAGGGGTTGTAATAACACGTCGGTTCCTTTGCCAATATGACCCACGCCCACTAGCGAAAATGCTTTGGcactacct                                                                                   >  2:510571/1‑244 (MQ=255)                                                                                                                                                                    gTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGt                                                                               >  1:369175/1‑249 (MQ=255)                                                                                                                                                                    gTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGATATTCGCGCTCCACCGCAACAATATTCGGCGGAGTTTTGACAGGAATCAGATCGCATTAGCATTTGCGCCAGTCGAGCACTTAACACAGCTGTTAGAGCACTGCCGGAGAGTCATGCAATATTGGTTATAATACCAAGTCGGTTTCTTTGCCACGATTACCCCAGCAAACTTGCGAAACTGCTTTAGCGATACCATTGt                                                                               >  1:495111/1‑249 (MQ=255)                                                                                                                                                                       tGTCCATCAATGGCCCGAGCGAGGTAAAGAGTTCATACGTAGCAGGGAAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTAGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGAtt                                                                          >  1:365771/1‑251 (MQ=255)                                                                                                                                                                               agtGGCCCGAGCGAGGCAAAGCGTGCATGCGGGGCGGGGAAGTCCGGCTCTCCCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCGTTGCGCCAGTCGAGCACTTCATACGGCTGTGAGAGCTCTGCCGGAGATGCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTc                                                                  <  1:290419/249‑1 (MQ=255)                                                                                                                                                                                             ggTAAAGGTTTCATACTGAGCATGGGAGTCGCCCTCCACCGCCACAATATGCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTTCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGGGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTGTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTccc                                                    <  1:227260/251‑1 (MQ=255)                                                                                                                                                                                                 aaGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGGGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGCTACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGc                                                >  2:245051/1‑251 (MQ=255)                                                                                                                                                                                                  aGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGTGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAAAGCTCTGCCGGTGAGTCATGCAGTAGGGGTTGTAATATCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGAGAAAATGCTTTAGCGATAATTTTGTCGGTTTCCCAGTCGCTGCGTGATTCCCACGcc                                               >  2:286348/1‑251 (MQ=255)                                                                                                                                                                                                     gTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGCGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCCCACTTTCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGAGATTCCCACGCCggc                                            >  2:358675/1‑250 (MQ=255)                                                                                                                                                                                                                  aGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTAGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCCCACGGCTGTGAGAGCTCTGCCGGCGAGTCATGCAGTAGTGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGATAATGCTTTGGCGATACCTTTTTAGGTTTCCCCGTCGCTGCGTGATTCCCACGCCGGCTCGACcgccgc                                 >  1:326139/1‑249 (MQ=255)                                                                                                                                                                                                                                   ccGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAAc              >  2:55676/1‑251 (MQ=255)                                                                                                                                                                                                                                     gCCACAATATTCGGCGGGTTTTTGCCTGGATTCAGATCGCATTCGCCTTTGCGCCAGGCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACCCTTGCGAAAATGCTTTGGCGATCCCTTTGTAGATTTCCCAGTCGCTGCGGGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACgg            <  1:326908/251‑1 (MQ=255)                                                                                                                                                                                                                                     cccACAATATTCGGGGGGGTTTTGCCTGGAAGCAGTGCGCATTCGCCTTTGCGCCAGTCGAGCACTTCCCACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACCCTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGTTGCGTGATTCCCACGCCGGATCGTCCGCCGCAGAAAGCGGATGAATAAACgg            <  2:43703/250‑1 (MQ=255)                                                                                                                                                                                                                                       cacaATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGAt          <  1:402121/251‑1 (MQ=255)                                                                                                                                                                                                                                             atTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGACGCAGAAAGCGGATGAATAAACGGATGCatat    >  2:227953/1‑251 (MQ=255)                                                                                                                                                                                                                                                cggcggGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCAGGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTGGCGAAAATGCTTTGGCGATACCTTTGGAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGATGCATATcc  <  1:497957/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     TTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGATGCATATCC  >  NC_000913/1539703‑1540191

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