Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,542,804	C→T	P419P (CCG→CCA)	narU ←	nitrate/nitrite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,542,804	0	C	T	100.0%	167.7 / NA	53	P419P (CCG→CCA)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/27);  total (26/27)

TTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGCGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCTTTATGAGCTTGCTCATCACTACCGCCTT  >  NC_000913/1542678‑1542907                                                                                                                               |                                                                                                        ttGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAg                                                                                             >  1:1478639/1‑139 (MQ=255)   ggCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAgcc                                                                                           >  1:486289/1‑139 (MQ=255)   ggCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAgcc                                                                                           >  2:1400675/1‑139 (MQ=255)   ggCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAgcc                                                                                           >  2:967046/1‑139 (MQ=255)    gCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAgccg                                                                                          <  1:967046/139‑1 (MQ=255)      tGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAgccgcc                                                                                        <  1:374477/139‑1 (MQ=255)        aaCTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCa                                                                                      >  2:1405424/1‑139 (MQ=255)            tCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTgc                                                                                  <  2:1478639/139‑1 (MQ=255)             ccgccgACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAgc                                                                                            >  2:648106/1‑128 (MQ=255)             ccgccgACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAgc                                                                                            <  1:648106/128‑1 (MQ=255)                ccgACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCa                                                                              <  1:1405424/139‑1 (MQ=255)                      aTAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCt                                                                        <  2:1196458/139‑1 (MQ=255)                               ccagGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATaaaa                                                                                                <  1:721217/106‑1 (MQ=255)                               ccagGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATaaaa                                                                                                >  2:721217/1‑106 (MQ=255)                               ccagGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAg                                                               >  1:1317509/1‑139 (MQ=255)                                   gTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCa                                                           >  2:49323/1‑139 (MQ=255)                                    tCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCag                                                          <  2:486289/139‑1 (MQ=255)                                                      tGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGt                                        >  2:1019668/1‑139 (MQ=255)                                                           aagaTTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAgccgc                                                    >  2:904467/1‑122 (MQ=255)                                                           aagaTTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAgccgc                                                    <  1:904467/122‑1 (MQ=255)                                                            agaTTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGc                                  <  1:1400675/139‑1 (MQ=255)                                                              aTTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAgccgc                                                    >  2:852054/1‑119 (MQ=255)                                                              aTTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAgccgc                                                    <  1:852054/119‑1 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTgg                                                                                                       <  2:1254272/66‑1 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCa                                                                              >  2:169300/1‑91 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCa                                                                              <  1:169300/91‑1 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGc                                                                         >  2:114146/1‑96 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGc                                                                         <  1:114146/96‑1 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCt                                                                        <  2:234809/97‑1 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCt                                                                        >  1:234809/1‑97 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGcttct                              >  2:1679654/1‑139 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGcttct                              >  2:103862/1‑139 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGctt                                >  2:529666/1‑137 (MQ=255)                                                                tAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGctt                                <  1:529666/137‑1 (MQ=255)                                                                  aaaaCACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCttt                            >  2:1123067/1‑139 (MQ=255)                                                                         tttCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGcc                                                                               <  2:1283499/81‑1 (MQ=255)                                                                         tttCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGcc                                                                               >  1:1283499/1‑81 (MQ=255)                                                                          ttCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATaaaa                                                                                                <  2:1416199/63‑1 (MQ=255)                                                                          ttCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATaaaa                                                                                                >  1:1416199/1‑63 (MQ=255)                                                                           tCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGATGGCTTCTTTATGAGCTTg                   <  1:1019668/139‑1 (MQ=255)                                                                             aTCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCTTTATGAGCTTGCt                 <  2:209502/139‑1 (MQ=255)                                                                                 cgcCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGctt                                <  2:1790749/120‑1 (MQ=255)                                                                                 cgcCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGctt                                >  1:1790749/1‑120 (MQ=255)                                                                                   ccGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCTTTATGAGCTTGCTCATCAc           >  1:696548/1‑139 (MQ=255)                                                                                      aCCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCTTTATGAGCTTGCTCATCACTAc        <  1:1123067/139‑1 (MQ=255)                                                                                       ccGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCTTTATGAGCTTGCTCATCACTAcc       <  1:103862/139‑1 (MQ=255)                                                                                            gagCCGGTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCTTTATGAGCTTGCTCATCACTACCGCCtt  <  2:1317509/139‑1 (MQ=255)                                                                                                  gTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCt                                                                        <  1:1645940/63‑1 (MQ=255)                                                                                                  gTCATATTGAGCGACATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCt                                                                        >  2:1645940/1‑63 (MQ=255)                                                                                                               aCATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGa                                      >  2:1783635/1‑84 (MQ=255)                                                                                                               aCATGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGa                                      <  1:1783635/84‑1 (MQ=255)                                                                                                                 aTGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCTTTATGAGCTTGCt                 <  1:931961/103‑1 (MQ=255)                                                                                                                 aTGCCAAACGCCTGTGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCTTTATGAGCTTGCt                 >  2:931961/1‑103 (MQ=255)                                                                                                                               |                                                                                                        TTGGCTGAACTTCCGCCGACCATAAACCAGCCAGGTCAGCAGCACACAAACGATGTAGAAGATTAAAAACACTTTCATCGCGCCGACCGGAGAGCCGGTCATATTGAGCGACATGCCAAACGCCTGCGGAATAAAAAAGCCGCCCACTGCGCCAATGGCTGAGATAAAGCCCAGAGCCGCCGCCGTTTCGGTGACGGCTTCTTTATGAGCTTGCTCATCACTACCGCCTT  >  NC_000913/1542678‑1542907

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