Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,542,644	A→G	intergenic (‑54/+28)	narZ ← / ← narU	nitrate reductase Z subunit alpha/nitrate/nitrite transporter NarU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,542,644	0	A	G	100.0%	185.2 / NA	53	intergenic (‑54/+28)	narZ/narU	nitrate reductase Z subunit alpha/nitrate/nitrite transporter NarU
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/27);  total (26/27)

ATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGACACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCAT  >  NC_000913/1542514‑1542779                                                                                                                                   |                                                                                                                                        aTCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGt                                                                                                                                <  1:383367/140‑1 (MQ=255)     ccTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACAt                                                                                                                            >  1:160906/1‑140 (MQ=255)     ccTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACAt                                                                                                                            >  2:699915/1‑140 (MQ=255)     ccTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACAt                                                                                                                            <  1:563801/140‑1 (MQ=255)      cTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGgca                                                                                                                                       >  1:219697/1‑128 (MQ=255)      cTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGgca                                                                                                                                       <  2:219697/128‑1 (MQ=255)      cTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGgca                                                                                                                                       >  1:609386/1‑128 (MQ=255)      cTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATc                                                                                                                           <  2:553870/140‑1 (MQ=255)           cGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCa                                                                                                                      >  1:565665/1‑140 (MQ=255)            gTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCAt                                                                                                                     <  1:781702/140‑1 (MQ=255)             tGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATc                                                                                                                    >  1:188772/1‑140 (MQ=255)                   tCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGAtctc                                                                                                              >  2:139595/1‑140 (MQ=255)                   tCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGAtctc                                                                                                              >  2:122823/1‑140 (MQ=255)                        aaaGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTAtt                                                                                                         <  2:306188/140‑1 (MQ=255)                            gTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTg                                                                                                     >  2:544852/1‑140 (MQ=255)                              ttCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGc                                                                                                   <  1:544852/140‑1 (MQ=255)                                          tttgAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTccgcc                                                                                       <  2:188772/140‑1 (MQ=255)                                          tttgAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTccgcc                                                                                       <  1:699915/140‑1 (MQ=255)                                               aGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCa                                                                                  >  1:722612/1‑140 (MQ=255)                                               aGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCa                                                                                  >  1:344584/1‑140 (MQ=255)                                               aGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCa                                                                                  >  1:63991/1‑140 (MQ=255)                                                 tAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATa                                                                                >  1:699282/1‑140 (MQ=255)                                                  aGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATaa                                                                               <  1:283110/140‑1 (MQ=255)                                                      aaaGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAAcca                                                                           <  2:294956/140‑1 (MQ=255)                                                        aGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAAccagc                                                                         <  1:438684/140‑1 (MQ=255)                                                             tCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGt                                                                    >  1:312912/1‑140 (MQ=255)                                                                aaaaGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcag                                                                 <  2:344584/140‑1 (MQ=255)                                                                aaaaGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcag                                                                 <  2:722612/140‑1 (MQ=255)                                                                 aaaGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcagc                                                                <  1:476672/140‑1 (MQ=255)                                                                 aaaGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcagc                                                                <  2:160906/140‑1 (MQ=255)                                                                     tttACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAgcagca                                                            <  2:312912/140‑1 (MQ=255)                                                                       tACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcagcag                                                              >  1:683075/1‑136 (MQ=255)                                                                       tACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGcaca                                                          >  2:416125/1‑140 (MQ=255)                                                                       tACTCATGACATGACTCCTGCTCCAGGAATG‑ATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcag                                                                 >  2:659386/1‑132 (MQ=255)                                                                       tACTCATGACATGACTCCTGCTCCAGGAATG‑ATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcag                                                                 <  1:659386/132‑1 (MQ=255)                                                                       tACTCATGACATGACTCCTGCTCCAAGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcagcag                                                              <  2:683075/136‑1 (MQ=255)                                                                         cTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGcacaca                                                        <  2:565665/140‑1 (MQ=255)                                                                          tCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACaa                                                       >  2:644604/1‑140 (MQ=255)                                                                           catgacatgaCTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACaaa                                                      <  1:624421/140‑1 (MQ=255)                                                                                    acTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcagca                                                               >  2:778589/1‑122 (MQ=255)                                                                                    acTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTcagca                                                               <  1:778589/122‑1 (MQ=255)                                                                                         tGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGAtt                                        >  2:230934/1‑140 (MQ=255)                                                                                               aGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAac                                  >  2:29100/1‑140 (MQ=255)                                                                                                ggAATGAATTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTA‑TTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAacac                                >  2:478212/1‑140 (MQ=255)                                                                                                 gaatgaatTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAacac                                <  1:644604/140‑1 (MQ=255)                                                                                                      aatTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCa                           >  2:290443/1‑140 (MQ=255)                                                                                                      aatTGGGATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCa                           >  1:255565/1‑140 (MQ=255)                                                                                                           ggATCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATcgcg                      >  1:665085/1‑140 (MQ=255)                                                                                                             aTCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGcc                    <  1:416125/140‑1 (MQ=255)                                                                                                             aTCTGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGcc                    <  1:29100/140‑1 (MQ=255)                                                                                                                tGATGCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGAc                 >  1:501583/1‑140 (MQ=255)                                                                                                                    gCCTTGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGga             <  2:676844/140‑1 (MQ=255)                                                                                                                       ttGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTA‑TTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGc         <  1:478212/140‑1 (MQ=255)                                                                                                                        tGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGc         <  2:665085/140‑1 (MQ=255)                                                                                                                        tGTCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGc         <  2:699282/140‑1 (MQ=255)                                                                                                                          tCGGATGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCg       >  2:462492/1‑140 (MQ=255)                                                                                                                               tGCGGCACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCat  >  1:494734/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                        ATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGGGATCTGATGCCTTGTCGGATGCGACACGGTTGTACATCAGGCATCGATCTCTTATTTTTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCAT  >  NC_000913/1542514‑1542779

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