Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,282,055	A→G	Y731C (TAC→TGC)	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,282,055	0	A	G	100.0%	145.9 / NA	41	Y731C (TAC→TGC)	narG	nitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/21);  total (20/21)

AATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTC  >  NC_000913/1281920‑1282189                                                                                                                                        |                                                                                                                                       aaTCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCctg                                                                                                                                    >  2:24752/1‑140 (MQ=255)   tCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCctgct                                                                                                                                  >  2:152670/1‑140 (MQ=255)          aaGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACg                                                                                                                           <  1:24752/140‑1 (MQ=255)              gggTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGa                                                                                                                             >  2:662141/1‑134 (MQ=255)              gggTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGa                                                                                                                             <  1:662141/134‑1 (MQ=255)               ggTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCa                                                                                                                      >  1:340495/1‑140 (MQ=255)                 ttCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACg                                                                                                                    <  1:204410/140‑1 (MQ=255)                      tCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATc                                                                                                               <  2:209858/140‑1 (MQ=255)                        cGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCa                                                                                                             <  2:192131/140‑1 (MQ=255)                           tttGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATc                                                                                                               >  2:685456/1‑135 (MQ=255)                           tttGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATc                                                                                                               <  1:685456/135‑1 (MQ=255)                              gCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTaa                                                                                                       >  2:687273/1‑140 (MQ=255)                              gCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTaa                                                                                                       <  2:569284/140‑1 (MQ=255)                                        aacCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTgggg                                                                                             <  1:40704/140‑1 (MQ=255)                                               aaaCGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAgg                                                                                      <  2:425286/140‑1 (MQ=255)                                                      aaaaaCCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGt                                                                               >  2:576047/1‑140 (MQ=255)                                                       aaaaCCACCCGCACAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTg                                                                              <  1:746224/140‑1 (MQ=255)                                                             accCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCg                                                                        >  2:51637/1‑140 (MQ=255)                                                             accCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGGTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCg                                                                        >  2:808262/1‑140 (MQ=255)                                                              cccGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCgg                                                                       <  1:687273/140‑1 (MQ=255)                                                              cccGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCgg                                                                       <  2:340495/140‑1 (MQ=255)                                                                  cgcAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGgaaga                                                                   >  2:129921/1‑140 (MQ=255)                                                                            tCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTgg                                                         <  1:152670/140‑1 (MQ=255)                                                                             cATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGc                                                        >  1:400811/1‑140 (MQ=255)                                                                                 tGGCGTTCTAACCTGCTCGGTTCTTCCAGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGa                                                    >  1:212649/1‑140 (MQ=255)                                                                                  ggCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGAc                                                   >  2:464532/1‑140 (MQ=255)                                                                                          aaCCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACa                                                  <  2:395748/133‑1 (MQ=255)                                                                                          aaCCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACa                                                  >  1:395748/1‑133 (MQ=255)                                                                                               gCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAg                                      <  1:51637/140‑1 (MQ=255)                                                                                               gCTCGGTTCTTCCGGTAAAGGTCATGAGGTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAg                                      <  1:808262/140‑1 (MQ=255)                                                                                                cTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAgg                                     >  1:310044/1‑140 (MQ=255)                                                                                                       cttcCGGTAAAGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTg                              <  2:400811/140‑1 (MQ=255)                                                                                                               aaaGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACTGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtggt                      >  2:207204/1‑140 (MQ=255)                                                                                                               aaaGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtggt                      >  2:375177/1‑140 (MQ=255)                                                                                                                 aGGTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCtggtgg                    <  1:129921/140‑1 (MQ=255)                                                                                                                   gTCATGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGtt                  >  1:663654/1‑140 (MQ=255)                                                                                                                     cATGAGTTTATGCTCAAGTGCCTGCTGGGGACTGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTAc                <  1:207204/140‑1 (MQ=255)                                                                                                                      aTGAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACg               >  2:77246/1‑140 (MQ=255)                                                                                                                        gAGTTTATGCTCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACGGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCt             <  2:68567/140‑1 (MQ=255)                                                                                                                                  tCAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGt   <  1:464532/140‑1 (MQ=255)                                                                                                                                   cAAGTGCCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTc  >  2:411299/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       AATCCCTGAAAGAGGGTTCCATCCGTTTTGCGGCAGAACAACCAGAAAACGGTAAAAACCACCCGCGCAACCTGTTCATCTGGCGTTCTAACCTGCTCGGTTCTTCCGGTAAAGGTCATGAGTTTATGCTCAAGTACCTGCTGGGGACGGAGCACGGTATCCAGGGTAAAGATCTGGGGCAACAGGGCGGCGTGAAGCCGGAAGAAGTGGACTGGCAGGACAATGGTCTGGAAGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTC  >  NC_000913/1281920‑1282189

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