Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,541,351	C→T	E414K (GAA→AAA)	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,541,351	0	C	T	100.0%	144.8 / NA	41	E414K (GAA→AAA)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/20);  total (21/20)

ACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCT  >  NC_000913/1541236‑1541486                                                                                                                    |                                                                                                                                        aCCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGga                                                                                                                 >  1:667877/1‑140 (MQ=255)           gATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgattccagattc                                                                                                       >  1:1071771/1‑140 (MQ=255)           gATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgattccagattc                                                                                                       >  2:138974/1‑140 (MQ=255)              tgtgTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCAt                                                                                                    >  1:54871/1‑140 (MQ=255)                tgtTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgattccaga                                                                                                          <  2:543092/132‑1 (MQ=255)                tgtTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgattccaga                                                                                                          >  1:543092/1‑132 (MQ=255)                       cGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCtttt                                                                                           >  1:602623/1‑140 (MQ=255)                        gCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCttttt                                                                                          <  2:602623/140‑1 (MQ=255)                             gAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTccc                                                                                     <  2:640255/140‑1 (MQ=255)                             gAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTccc                                                                                     <  1:772196/140‑1 (MQ=255)                                     ggATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGaa                                                                             <  2:465260/140‑1 (MQ=255)                                      gATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGaaa                                                                            >  1:334780/1‑140 (MQ=255)                                       aTTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAAc                                                                           <  2:83002/140‑1 (MQ=255)                                          ttCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCaa                                                                        <  2:1042271/140‑1 (MQ=255)                                               tGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGaa                                                                   <  1:1055915/140‑1 (MQ=255)                                               tGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGaa                                                                   <  2:1071771/140‑1 (MQ=255)                                                 ccgccAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGTGGAAACCAATCGAAcc                                                                 >  1:731975/1‑140 (MQ=255)                                                       aaGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCgg                                                           <  2:143374/140‑1 (MQ=255)                                                        aGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGc                                                          <  2:54871/140‑1 (MQ=255)                                                         gTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCa                                                         <  1:1226645/140‑1 (MQ=255)                                                                  ggCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGc                                                          <  1:488346/130‑1 (MQ=255)                                                                  ggCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGc                                                          <  1:528588/130‑1 (MQ=255)                                                                  ggCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGc                                                          >  2:488346/1‑130 (MQ=255)                                                                  ggCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGc                                                          >  2:528588/1‑130 (MQ=255)                                                                  ggCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAAtt                                                >  1:906128/1‑140 (MQ=255)                                                                   gCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTc                                               >  2:790821/1‑140 (MQ=255)                                                                   gCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTCTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTc                                               >  1:998229/1‑140 (MQ=255)                                                                       cGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCt                                           <  1:138974/140‑1 (MQ=255)                                                                        gCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTg                                          <  1:728832/140‑1 (MQ=255)                                                                        gCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTg                                          <  1:296672/140‑1 (MQ=255)                                                                         ccTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGc                                         <  2:103162/140‑1 (MQ=255)                                                                              aaCAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTAt                                    >  1:495105/1‑140 (MQ=255)                                                                                    gTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTg                              >  1:1143883/1‑140 (MQ=255)                                                                                        tGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTAc                          >  1:78061/1‑140 (MQ=255)                                                                                         gTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACg                         <  1:381964/140‑1 (MQ=255)                                                                                            gACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGtt                      >  2:371051/1‑140 (MQ=255)                                                                                                 aGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCa                 >  2:542681/1‑140 (MQ=255)                                                                                                 aGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCa                 >  2:375984/1‑140 (MQ=255)                                                                                                  gCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTAAGGTTTTCCAc                >  2:775673/1‑140 (MQ=255)                                                                                                                aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCt  >  1:331713/1‑140 (MQ=255)                                                                                                                aaTTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCt  <  1:790821/140‑1 (MQ=255)                                                                                                                    |                                                                                                                                        ACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCT  >  NC_000913/1541236‑1541486

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