Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,543,610	T→C	N151D (AAT→GAT)	narU ←	nitrate/nitrite transporter NarU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,543,610	0	T	C	100.0%	121.0 / NA	34	N151D (AAT→GAT)	narU	nitrate/nitrite transporter NarU
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/18);  total (16/18)

AGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCA  >  NC_000913/1543477‑1543740                                                                                                                                      |                                                                                                                                   agaAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCAt                                                                                                                              >  1:481606/1‑140 (MQ=255)           aTACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCa                                                                                                                    <  2:484196/140‑1 (MQ=255)                     aaaaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCAc                                                                                                          <  1:537936/140‑1 (MQ=255)                                gCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCg                                                                                               <  2:311802/140‑1 (MQ=255)                                  ggTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCa                                                                                             <  2:375566/140‑1 (MQ=255)                                     gCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAg                                                                                          <  2:335194/140‑1 (MQ=255)                                           aGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGc                                                                                    >  1:506455/1‑140 (MQ=255)                                                        ttacacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAat                                                                       <  2:315143/140‑1 (MQ=255)                                                          acacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAatat                                                                     >  2:323496/1‑140 (MQ=255)                                                          acacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAatat                                                                     >  1:62853/1‑140 (MQ=255)                                                          acacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAatat                                                                     >  2:23193/1‑140 (MQ=255)                                                           cacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATc                                                                    >  1:503970/1‑140 (MQ=255)                                                           cacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATc                                                                    >  2:567836/1‑140 (MQ=255)                                                                 aGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCaaaa                                                              >  1:573292/1‑140 (MQ=255)                                                                     tCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAg                                                          <  1:323496/140‑1 (MQ=255)                                                                     tCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAg                                                          >  2:552821/1‑140 (MQ=255)                                                                       cTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTa                                                        >  1:463159/1‑140 (MQ=255)                                                                             cGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGa                                                  <  2:463159/140‑1 (MQ=255)                                                                                  ttAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTg                                             <  2:337361/140‑1 (MQ=255)                                                                                        ccAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGc                                       <  2:481606/140‑1 (MQ=255)                                                                                        ccAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGAATTTGCACGGc                                       <  2:573292/140‑1 (MQ=255)                                                                                            gagCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAAtt                                   >  1:511991/1‑140 (MQ=255)                                                                                               cgcTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCg                                >  2:284785/1‑140 (MQ=255)                                                                                                 cTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAg                              >  1:169785/1‑140 (MQ=255)                                                                                                   cccTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGcc                            <  1:284785/140‑1 (MQ=255)                                                                                                     cTTGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAg                          <  2:503970/140‑1 (MQ=255)                                                                                                       tGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGAc                        >  1:415366/1‑140 (MQ=255)                                                                                                       tGTTTGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGAc                        >  1:197268/1‑140 (MQ=255)                                                                                                         tttGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGc                      >  1:413510/1‑140 (MQ=255)                                                                                                          ttGGCTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCa                     <  1:23193/140‑1 (MQ=255)                                                                                                              cTTTTGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCGGACGCAAGGa                 <  2:169785/140‑1 (MQ=255)                                                                                                                 ttGGAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGaata              <  1:177593/140‑1 (MQ=255)                                                                                                                   ggAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGTATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGaataat            <  2:415366/140‑1 (MQ=255)                                                                                                                   ggAAAGAAGAAACTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGaataat            <  2:197268/140‑1 (MQ=255)                                                                                                                             aaCTGATATCGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCa  <  2:413510/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   AGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCA  >  NC_000913/1543477‑1543740

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