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%	86.9 / NA	26	A882P (GCG→CCG)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (19/7);  total (19/7)

GTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTG‑AAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCC‑ATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGATGCATATCCGA  >  NC_000913/1539702‑1540193                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        gTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTG‑AAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTTTGCCGGGCC‑ATCACGCCTGGGGTAATTGAGTTTTCCGAGGGAATCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCCCC‑ATTGCCCAGTTTGTCCATCAATGGCCCGGGCGAGGGAAAGCGTTCCTACGTAGCCAGGG‑AGTCGCGCCCCACCGCCCCCATATTCGGCGGGGttt                                                                                                                                                                                                                                                    >  2:169785/1‑251 (MQ=255)                 cagtcccaaaatccTTCTGGAAGCGTAAATGGCGGTGTAATTAAGGGCACCCCCCTGT‑CCGGGCCCATCTCGTTTGGTTAAATGGATTTTGCCGAGAAATTCTATTTCACCGGGGGTTT‑CCCAC‑AAATCCCCTTTCCCGCC‑ATTGCCAAGTGTTTCCATCAATGCCCGGTGCGAGGTAAGGTTTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGAt                                                                                                                                                                                                                                     <  1:325337/238‑1 (MQ=255)                    gccagAATCACTTCTG‑AAGCCTCAATGGCGGTGTCAATCAGCGGACGCCCCTATGCCGGGCC‑ATCACGCTTGGAGTACTTGAGTTTGCCGAGGAAATCGGTTTCCTCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGGCGCGCTCCACCGCCACACAATTCGGCCGGGTTTTGCCTGGCATTAGAttgc                                                                                                                                                                                                                                  >  2:246146/1‑247 (MQ=255)                              cTTCTG‑AAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCC‑ATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGGATCAGATCGCAATCGCCTTTg                                                                                                                                                                                                                       >  1:78800/1‑251 (MQ=255)                                                                   gCCCCTGTGCCGGGCC‑ATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGctct                                                                                                                                                                                   >  1:28201/1‑251 (MQ=255)                                                                   gCCCCTGTGCCGGGCC‑ATCAAGCTTGGAGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGAAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGctct                                                                                                                                                                                   >  2:503248/1‑251 (MQ=255)                                                                                        cGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGt                                                                                                                                                                >  1:64325/1‑249 (MQ=255)                                                                                             ggTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGATGTGAGAGCTCTGCCGGAG‑AGTCATGCGGCAGGGGtt                                                                                                                                                         >  1:250369/1‑251 (MQ=255)                                                                                                  aaTTGAGTTTGCTGAGGAAATCGATTTCATCCTGGGTAT‑TCCACGAACTCCC‑TTTACCGCA‑CTTGCCACGTTTGTCCATCAATGGACCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATCCGGCGGGGCTTTGCCTGGCATCAGATCGCCTTCGGCATTGAGCCAGTCGAGCACTTCACACGGCTGTGAGCGCTTTGCCGGGGCAG‑CACGCCGTAGGGGTTGTAAt                                                                                                                                                    >  2:408787/1‑251 (MQ=255)                                                                                                           ggcgaagaaaatgattttCATCC‑GGGTATGCCCACAAAATCC‑‑TTTACCCCCCATTGCCAAGTGTGTCCATCAATGGCCGGAGGGAGGAAAAGCGTTCAGACTGAGC‑GGGGGAGCCGCGCCCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCg                                                                                                                                           <  2:28201/236‑1 (MQ=255)                                                                                                               gAGGAAATCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGttt                                                                                                                                       >  1:192318/1‑251 (MQ=255)                                                                                                                   aaaTCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATAGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTCATACCACGTCGGTTTCttt                                                                                                                                   >  2:82285/1‑251 (MQ=255)                                                                                                                                          tCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACt                                                                                                             >  1:378948/1‑251 (MQ=255)                                                                                                                                                                              gTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGAttt                                                                           >  1:106554/1‑251 (MQ=255)                                                                                                                                                                              gTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACACACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGAttt                                                                           >  1:283590/1‑251 (MQ=255)                                                                                                                                                                                    cAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATAACTTTGTAGATTTCCCAGt                                                                     >  2:184541/1‑251 (MQ=255)                                                                                                                                                                                            cgagcgagGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑GGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCACGATGACCCACGCACCCTTGCGAAAATGCTTTGGCGATAACTTTGTAGATTTCCCAGTCGCTGCGt                                                             >  2:466423/1‑251 (MQ=255)                                                                                                                                                                                                ggagGAAAAGCGTTCAGCCGTAGC‑AGGGTAGTCGCGCCCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGAtt                                                         <  1:82285/250‑1 (MQ=255)                                                                                                                                                                                                         gcgttaacaggTAGC‑GGGGGAGCCGCGCCCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCg                                                <  2:192318/241‑1 (MQ=255)                                                                                                                                                                                                                 aCGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGAcc                                        >  1:155323/1‑251 (MQ=255)                                                                                                                                                                                                                                 cgcgCTCCACCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCggg                         >  2:235810/1‑250 (MQ=255)                                                                                                                                                                                                                                    ccTCCCCCGCCACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGa                      <  1:385985/250‑1 (MQ=255)                                                                                                                                                                                                                                            gCCACAATATCCGGCGGCGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGCG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCCGAAAGCGGATGAATAAACgg              >  2:228295/1‑251 (MQ=255)                                                                                                                                                                                                                                             ccACAATATTCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGa             <  2:378948/251‑1 (MQ=255)                                                                                                                                                                                                                                              cacaATATCCGGCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTCCCCAGTCGCTGCGTGATTCCCACGCCGGAGCGACCGCCGCAGAAAGCGGATGAATAAACGGAt            <  2:106554/251‑1 (MQ=255)                                                                                                                                                                                                                                                        ggCGGGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCAAACGGCTGTGCGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGCTTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGCCCGCCGCAGAAAGCGGATGAATAAACGGCTGCATATCCgg  >  2:448520/1‑250 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        GTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTG‑AAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCC‑ATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTAT‑TCCACGAAATCCC‑TTTACCGCC‑ATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGC‑AGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAG‑AGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGATGCATATCCGA  >  NC_000913/1539702‑1540193

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