Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,276,251	Δ1 bp	coding (932/3744 nt)	narG →	nitrate reductase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,276,249	0	G	.	100.0%	221.8 / NA	50	coding (930/3744 nt)	narG	nitrate reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (25/25);  total (25/25)

CACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCC  >  NZ_CP009273/1276162‑1276332                                                                                        |                                                                                    cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGTTGGCAACGAACCAGGGCACCAATGCGGCAATGGCGCTGGCGATggc                                                                                   <  2:1565529/90‑2 (MQ=255) cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  2:2814789/90‑2 (MQ=255) cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   >  1:31443/1‑89 (MQ=255) cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  2:2483159/90‑2 (MQ=255) cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  2:283892/90‑2 (MQ=255) cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  1:1535504/90‑2 (MQ=255) cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  1:1022922/90‑2 (MQ=255)       aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTa                                                                            >  2:2031114/1‑90 (MQ=255)        gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTaa                                                                           <  2:2126048/90‑1 (MQ=255)        gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTaa                                                                           <  2:197686/90‑1 (MQ=255)        gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTaa                                                                           >  2:1301864/1‑90 (MQ=255)            aCGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGACACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAAtgct                                                                       <  2:776072/90‑1 (MQ=255)             cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAAtgctg                                                                      >  1:1604056/1‑90 (MQ=255)             cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAAtgctg                                                                      >  2:2894095/1‑90 (MQ=255)             cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAAtgctg                                                                      >  2:2226042/1‑90 (MQ=255)                 gAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTg                                                                  >  1:1349342/1‑90 (MQ=255)                     tCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAAtt                                                              <  2:961374/90‑1 (MQ=255)                     tCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAAtt                                                              <  1:1494556/90‑1 (MQ=255)                               gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  2:2553358/60‑2 (MQ=255)                               gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   >  1:2553358/1‑59 (MQ=255)                               gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                                    <  1:1020983/90‑1 (MQ=255)                               gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                                    <  2:1604056/90‑1 (MQ=255)                                    aTCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAAccc                                               <  2:2182158/90‑1 (MQ=255)                                    aTCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAAccc                                               <  2:1496395/90‑1 (MQ=255)                                      cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGt                                                                             >  2:615474/1‑58 (MQ=255)                                      cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGt                                                                             <  1:615474/58‑1 (MQ=255)                                      cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCaa                                             >  2:2724357/1‑90 (MQ=255)                                       tgtgGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGTTGCGTGAATTCCACCTCGACAACCCAAg                                            <  1:498525/90‑1 (MQ=255)                                         tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                          >  2:490340/1‑90 (MQ=255)                                         tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                          >  2:1339515/1‑90 (MQ=255)                                         tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                          <  2:352821/90‑1 (MQ=255)                                         tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                          <  1:177726/90‑1 (MQ=255)                                         tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                          >  1:180099/1‑90 (MQ=255)                                              ggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                          <  1:2984714/85‑1 (MQ=255)                                              ggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                          >  2:2984714/1‑85 (MQ=255)                                              ggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt                                     >  2:1745546/1‑90 (MQ=255)                                               gcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAtt                                    >  1:1869773/1‑90 (MQ=255)                                                  ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa                                 >  2:366929/1‑90 (MQ=255)                                                   cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAc                                <  1:46339/90‑1 (MQ=255)                                                   cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAc                                <  1:2556485/90‑1 (MQ=255)                                                          gggCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                                    >  2:743370/1‑63 (MQ=255)                                                          gggCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                                    <  1:743370/63‑1 (MQ=255)                                                          gggCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAt                         >  2:895656/1‑90 (MQ=255)                                                          gggCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAt                         >  1:987707/1‑90 (MQ=255)                                                            gCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTAtgt                       >  2:650098/1‑90 (MQ=255)                                                              aCCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc                     <  1:2724357/90‑1 (MQ=255)                                                              aCCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc                     <  1:2796828/90‑1 (MQ=255)                                                               ccGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCg                    >  1:92661/1‑90 (MQ=255)                                                               ccGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCg                    >  1:1069191/1‑90 (MQ=255)                                                                cGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGt                   <  2:1813052/90‑1 (MQ=255)                                                                cGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCGCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGt                   <  2:1063459/90‑1 (MQ=255)                                                                  aTGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCg                 >  2:2821030/1‑90 (MQ=255)                                                                      ggCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTac             <  2:711466/90‑1 (MQ=255)                                                                      ggCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTac             <  1:2894095/90‑1 (MQ=255)                                                                        cAATGGCGCTGGCGAT‑GGCCCCGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTacac           <  1:2031114/90‑1 (MQ=255)                                                                         aaTGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACAcc          >  2:2783561/1‑90 (MQ=255)                                                                         aaTGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACAcc          >  1:1828984/1‑90 (MQ=255)                                                                            ggCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGAc       >  1:76778/1‑90 (MQ=255)                                                                                 tGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCCCCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGcc  <  1:2226042/90‑1 (MQ=255)                                                                                        |                                                                                    CACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCC  >  NZ_CP009273/1276162‑1276332

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