Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,539,343	T→C	E1083G (GAA→GGA)	narZ ←	nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,539,343	0	T	C	100.0%	140.0 / NA	38	E1083G (GAA→GGA)	narZ	nitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/18);  total (20/18)

ATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCGACAGCGTGCG  >  NC_000913/1539198‑1539485                                                                                                                                                  |                                                                                                                                               aTCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCggg                                                                                                                                           >  1:539640/1‑151 (MQ=255)      cAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGcc                                                                                                                                      >  1:437321/1‑151 (MQ=255)       aaTGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCg                                                                                                                                     <  1:243627/151‑1 (MQ=255)        aTGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGt                                                                                                                                    <  1:435601/151‑1 (MQ=255)               gTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCgg                                                                                                                             <  1:434657/151‑1 (MQ=255)                    cGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTAt                                                                                                                          <  1:617701/149‑1 (MQ=255)                     gggCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTg                                                                                                                       <  2:181908/151‑1 (MQ=255)                                   cTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTg                                                                                                          >  1:614287/1‑150 (MQ=255)                                     gATCCAGACCATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGAc                                                                                                        >  1:543637/1‑150 (MQ=255)                                      aTCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGc                                                                                                      <  2:29395/151‑1 (MQ=255)                                        ccAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCt                                                                                                     >  1:524485/1‑150 (MQ=255)                                                  ggTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCg                                                                                           >  1:410247/1‑150 (MQ=255)                                                    tCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGAt                                                                                         >  2:120661/1‑150 (MQ=255)                                                         cgcgAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcgg                                                                                   <  1:288631/151‑1 (MQ=255)                                                          gcgAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAggcggg                                                                                  >  1:516990/1‑151 (MQ=255)                                                               agagCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGAt                                                                              >  2:495281/1‑150 (MQ=255)                                                                   cGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCTGGCGATAAGcc                                                                         <  2:77133/151‑1 (MQ=255)                                                                         cATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAgg                                                                   >  1:575821/1‑151 (MQ=255)                                                                                 ggTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCAcc                                                           >  1:585672/1‑151 (MQ=255)                                                                                        aCTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCa                                                    <  2:281190/151‑1 (MQ=255)                                                                                           gTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCAc                                                   <  2:584374/149‑1 (MQ=255)                                                                                           gTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGc                                                 <  2:215687/151‑1 (MQ=255)                                                                                                   aGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACgg                                         <  2:581507/151‑1 (MQ=255)                                                                                                    gTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGga                                        <  2:222936/151‑1 (MQ=255)                                                                                                           ccccATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCt                                 <  1:289012/151‑1 (MQ=255)                                                                                                             ccATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTg                                >  2:282055/1‑150 (MQ=255)                                                                                                                    tGGTGCGGCGTCAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAgct                        >  1:102723/1‑151 (MQ=255)                                                                                                                               cAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTccc             >  1:382966/1‑151 (MQ=255)                                                                                                                               cAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTccc             >  2:15225/1‑151 (MQ=255)                                                                                                                               cAGGAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTcc              <  1:248393/150‑1 (MQ=255)                                                                                                                                 ggAAGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCg            >  2:205668/1‑150 (MQ=255)                                                                                                                                    aGTTAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCGAc          >  1:646533/1‑149 (MQ=255)                                                                                                                                       tAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCGACAGCg      >  1:101158/1‑150 (MQ=255)                                                                                                                                       tAAGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCGACAGCGt     >  2:168190/1‑151 (MQ=255)                                                                                                                                        aaGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCGACAGCGTg    <  1:461017/151‑1 (MQ=255)                                                                                                                                         aGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCGACAGCGTgc   >  2:144893/1‑151 (MQ=255)                                                                                                                                         aGTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCGACAGCGTgc   <  2:587455/151‑1 (MQ=255)                                                                                                                                          gTGCTTTTCCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCGACAGCGTgcg  <  2:172134/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               ATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTGGCGCATCTCACTGACGCTACGGGTGTCGATAGGCGGGCGATAAGCCACCAGGCTTTCACCAAAAGCACGCATCCACGGATGATCCTGATAGAGCTGCTGGCGTCCCGACAGCGTGCG  >  NC_000913/1539198‑1539485

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