Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,538,868	T→C	V1241V (GTA→GTG)	narZ ←	nitrate reductase Z subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,538,868	0	T	C	100.0%	139.6 / NA	41	V1241V (GTA→GTG)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/20);  total (21/20)

TGCGTACTCCATGCCTTCGCGCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATA  >  NC_000913/1538737‑1539003                                                                                                                                    |                                                                                                                                        tGCGTACTCCATGCCTTCGCGCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATc                                                                                                                                 >  2:469819/1‑140 (MQ=255) tGCGTACTCCATGCCTTCGCGCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGGGAACGGAATTTGATTTTTTCGGCTCCTGCACCTGATc                                                                                                                                 >  1:290503/1‑140 (MQ=37)            tGCCTTCGCGCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTtcat                                                                                                                      >  2:537144/1‑140 (MQ=255)            tGCCTTCGCGCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTtcat                                                                                                                      >  2:252478/1‑140 (MQ=255)            tGCCTTCGCGCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTtcat                                                                                                                      >  2:623900/1‑140 (MQ=255)              ccTTCGCGCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTtcatca                                                                                                                    <  1:469819/140‑1 (MQ=255)                  cgcgCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATcca                                                                                                                >  1:244374/1‑140 (MQ=255)                  cgcgCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATcca                                                                                                                <  1:596448/140‑1 (MQ=255)                  cgcgCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATcca                                                                                                                >  2:925589/1‑140 (MQ=255)                  cgcgCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGCATTTTCATGTTTACTCCTCCTGCACCTGATCCCGACCTTAATCATcca                                                                                                                >  2:459420/1‑140 (MQ=21)                    cgcCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATccagc                                                                                                              >  2:554901/1‑140 (MQ=255)                        cGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGt                                                                                                          <  2:749661/140‑1 (MQ=255)                            ccAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAAc                                                                                                      <  2:244374/140‑1 (MQ=255)                                          aGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGc                                                                                        <  1:535084/140‑1 (MQ=255)                                             tCACCGAACAGGTATGGTAGCCTATACATTTGTCGAGGTTAAGCACCAAGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGa                                                                                     <  2:410995/140‑1 (MQ=255)                                             tCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGa                                                                                     >  2:894287/1‑140 (MQ=255)                                              cACCGAACAGGTATGGCAGCCGTTACATTTGTCTAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGAt                                                                                    <  2:449551/140‑1 (MQ=255)                                               aCCGAACAGGTATGGCAGCCGAAACTTTTCTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATc                                                                                   <  1:339183/140‑1 (MQ=255)                                               aCCGAACAGGTATGGAAGCCGATACATTTGTCGAGGTGAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATc                                                                                   <  1:413356/140‑1 (MQ=255)                                                        gTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAACCAGTTAACGTCCTTCATCTTGCGGATCATGATGAAc                                                                          >  1:371799/1‑140 (MQ=255)                                                           tGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCCGTTAACGTTCTGCATCTTGCGTATCATGATGAACTCa                                                                       >  1:281393/1‑140 (MQ=255)                                                             caggCCGATACATTTGTCGAGGTTAAGCACCATGCCGCCTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATc                                                                     <  1:493097/137‑1 (MQ=255)                                                                gCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATcgcg                                                                  >  2:594969/1‑140 (MQ=255)                                                                 ccGATACATTTGTCGAGTTTAACCACCATGCTGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCgg                                                                 <  2:387772/140‑1 (MQ=255)                                                                  cGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGt                                                                <  2:934936/140‑1 (MQ=255)                                                                  cGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGt                                                                >  1:934936/1‑140 (MQ=255)                                                                   gATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGtt                                                               <  2:685362/140‑1 (MQ=255)                                                                           tGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAGCGTTCTTCGTCTTGAGGATCATGATAAACTAATCGCGGTGCGATcctg                                                       >  2:493097/1‑138 (MQ=38)                                                                           tGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGa                                                       <  1:594969/140‑1 (MQ=255)                                                                            gTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCGTGCGGATGATGGTGTACTCATTGCGGTTCGATCCGAc                                                      >  1:295955/1‑140 (MQ=38)                                                                              cGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACgg                                                    <  1:252478/140‑1 (MQ=255)                                                                                      aGCACCATGCCGACTTGCGAACTGATTTTCACTTTATCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCgtag                                            <  2:290503/140‑1 (MQ=255)                                                                                           cATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGtt                                       >  2:989932/1‑140 (MQ=255)                                                                                           cATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGtt                                       >  1:628838/1‑140 (MQ=255)                                                                                            aTGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGcc                                                >  2:101635/1‑130 (MQ=255)                                                                                            aTGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGcc                                                <  1:101635/130‑1 (MQ=255)                                                                                                       tgAACGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGc                           >  1:649089/1‑140 (MQ=255)                                                                                                          aCGGATTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAg                        <  2:119244/140‑1 (MQ=255)                                                                                                              aTTTTCATTTTTTCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTgc                    >  1:675865/1‑140 (MQ=255)                                                                                                                          tCGCCTCCTGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCa        <  1:554901/140‑1 (MQ=255)                                                                                                                                cctGCACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCata  <  1:537144/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        TGCGTACTCCATGCCTTCGCGCCCGGTCCAGACGTTTTTACAGGTCACCGAACAGGTATGGCAGCCGATACATTTGTCGAGGTTAAGCACCATGCCGACTTGTGAACGGATTTTCATTTTTTCGCCTCCTGTACCTGATCGCGACCTTCATCATCCAGCCAGTTAACGTTCTTCATCTTGCGGATCATGATGAACTCATCGCGGTTCGATCCGACGGTGCCGTAGTAGTTAAAGCCCCAGGCCAGCTGCGCGTAACCGCCAATCATA  >  NC_000913/1538737‑1539003

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