Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,283,325	T→C	R1154R (CGT→CGC)	narG →	nitrate reductase A subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,283,325	0	T	C	100.0%	157.0 / NA	45	R1154R (CGT→CGC)	narG	nitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/22);  total (23/22)

CTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCGAAACCGAC  >  NC_000913/1283188‑1283450                                                                                                                                          |                                                                                                                              ctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCgcgt                                                                                                                             >  1:726823/1‑140 (MQ=255)    aCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTcc                                                                                                                          >  2:775884/1‑140 (MQ=255)     cTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCg                                                                                                                         >  2:213108/1‑140 (MQ=255)            tCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGa                                                                                                                  <  2:159732/140‑1 (MQ=255)            tCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGa                                                                                                                  <  1:781869/140‑1 (MQ=255)              gcgGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATg                                                                                                                >  1:228864/1‑140 (MQ=255)              gcgGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATg                                                                                                                >  2:543132/1‑140 (MQ=255)               cggtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGa                                                                                                               <  2:788042/140‑1 (MQ=255)               cggtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGa                                                                                                               <  2:780675/140‑1 (MQ=255)               cggtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGa                                                                                                               <  1:225825/140‑1 (MQ=255)               cggtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGa                                                                                                               <  1:301830/140‑1 (MQ=255)               cggtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGa                                                                                                               >  2:75027/1‑140 (MQ=255)               cggtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGa                                                                                                               <  2:726823/140‑1 (MQ=255)                            gtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCa                                                                                                  >  2:213053/1‑140 (MQ=255)                              cTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCAcg                                                                                                <  2:526904/140‑1 (MQ=255)                              cTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCAcg                                                                                                >  2:35549/1‑140 (MQ=255)                              cTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCAcg                                                                                                <  1:213053/140‑1 (MQ=255)                                ggTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCAcgcg                                                                                              <  2:557910/140‑1 (MQ=255)                                       tGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAAc                                                                                       >  1:775858/1‑140 (MQ=255)                                       tGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAAc                                                                                       >  2:338748/1‑140 (MQ=255)                                        gAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACg                                                                                      <  2:713427/140‑1 (MQ=255)                                           gCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTAt                                                                                   >  1:470240/1‑140 (MQ=255)                                               aTGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                                                               <  1:318786/140‑1 (MQ=255)                                                      aGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctgc                                                                        <  1:338748/140‑1 (MQ=255)                                                              gTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCg                                                                <  1:543132/140‑1 (MQ=255)                                                                    ccGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAAtt                                                          <  1:35549/140‑1 (MQ=255)                                                                     cGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAAtca                                                         >  2:45630/1‑138 (MQ=255)                                                                          aCGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCaa                                                    <  2:27459/140‑1 (MQ=255)                                                                          aCGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCaa                                                    <  2:470240/140‑1 (MQ=255)                                                                                 gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATCACCCAACAGCgtg                                             <  1:45630/140‑1 (MQ=255)                                                                                  aTTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGATCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTTGGAAATTACCCAACAGCgagg                                            >  1:53095/1‑140 (MQ=255)                                                                                      aaGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCg                                                                                                                         <  1:100508/59‑1 (MQ=255)                                                                                      aaGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCg                                                                                                                         >  2:100508/1‑59 (MQ=255)                                                                                           ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCa                                   >  1:91000/1‑140 (MQ=255)                                                                                           ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGGTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCa                                   >  1:704834/1‑140 (MQ=255)                                                                                            tCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCAt                                  >  1:195751/1‑140 (MQ=255)                                                                                                   cAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTTCCCAACAGCGTGGTGGTATCCATAACTAgg                           >  2:493026/1‑140 (MQ=255)                                                                                                   cAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCgg                           >  2:327147/1‑140 (MQ=255)                                                                                                          gCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCg                    >  1:317847/1‑140 (MQ=255)                                                                                                             ctGACTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTAt                 <  2:704834/140‑1 (MQ=255)                                                                                                                aCTGCCCGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCAc              >  1:262211/1‑140 (MQ=255)                                                                                                                     ccGTGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCGa         >  2:581691/1‑140 (MQ=255)                                                                                                                        tGCGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCGAAAc      >  2:73743/1‑140 (MQ=255)                                                                                                                          cGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCGAAACCg    <  2:262211/140‑1 (MQ=255)                                                                                                                          cGGTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCGAAACCg    <  1:75027/140‑1 (MQ=255)                                                                                                                            gTTGTCAGCCAGCGCGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCGAAACCGAc  >  2:440244/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                              CTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCGTGGTGGTATCCATAACTCGGTCACCCGTATCACGCCGAAACCGAC  >  NC_000913/1283188‑1283450

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: ‑