Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,280,235	A→G	E124E (GAA→GAG)	narG →	nitrate reductase 1, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,280,235	0	A	G	100.0%	220.9 / NA	63	E124E (GAA→GAG)	narG	nitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/33);  total (30/33)

CGCGTACCCGTCCGGATCTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAAGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGACGCGGTGGATTTGTTCGTTCTTCCTGGCAGGA  >  NC_000913/1280099‑1280357                                                                                                                                         |                                                                                                                           cgcgTACCCGTCCGGATCTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGc                                                                                                                          >  2:436480/1‑139 (MQ=255)   cgTACCCGTCCGGATCTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGa                                                                                                                        >  1:797138/1‑139 (MQ=255)       cccGTCCGGATCTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGgc                                                                                                                    >  2:1193016/1‑139 (MQ=255)           tCCGGATCTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTg                                                                                                                >  2:1607806/1‑139 (MQ=255)                aTCTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACTGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAg                                                                                                           <  1:1193016/139‑1 (MQ=255)                aTCTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAg                                                                                                           <  1:1607806/139‑1 (MQ=255)                  cTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAgg                                                                                                                     >  1:787277/1‑127 (MQ=255)                  cTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAgg                                                                                                                     <  2:787277/127‑1 (MQ=255)                  cTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCg                                                                                                         >  2:1247923/1‑139 (MQ=255)                             tGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAg                                                                                              >  2:306641/1‑139 (MQ=255)                               aaCCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACGCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGc                                                                                            <  1:1247923/139‑1 (MQ=255)                                aCCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCa                                                                                           <  2:1253464/139‑1 (MQ=255)                                 ccTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCAt                                                                                          <  2:797138/139‑1 (MQ=255)                                      cgGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCt                                                                                     <  2:1207396/139‑1 (MQ=255)                                       ggCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCtt                                                                                    <  2:150099/139‑1 (MQ=255)                                                  cgGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGaa                                                                         <  2:549827/139‑1 (MQ=255)                                                               tCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGgcg                                                                                                                   >  2:743659/1‑84 (MQ=255)                                                               tCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGgcg                                                                                                                   <  1:743659/84‑1 (MQ=255)                                                                   ggTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGaa                                                        <  1:436480/139‑1 (MQ=255)                                                                                 gCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGAt                                                                                                       <  1:177376/78‑1 (MQ=255)                                                                                 gCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGAt                                                                                                       >  2:177376/1‑78 (MQ=255)                                                                                 gCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAg                                              <  2:1131056/135‑1 (MQ=255)                                                                                 gCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAg                                              >  1:1131056/1‑135 (MQ=255)                                                                                  ccAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGAt                                                                                                       >  2:215537/1‑77 (MQ=255)                                                                                  ccAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGAt                                                                                                       <  1:215537/77‑1 (MQ=255)                                                                                  ccAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCa                                            >  2:1621849/1‑136 (MQ=255)                                                                                  ccAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCa                                            <  1:1621849/136‑1 (MQ=255)                                                                                  ccAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGgc                                         >  2:1227235/1‑139 (MQ=255)                                                                                   cAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGa                                                                                                                        <  1:1363213/59‑1 (MQ=255)                                                                                   cAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGa                                                                                                                        >  2:1363213/1‑59 (MQ=255)                                                                                         ccTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTc                                                                                   <  2:27092/90‑1 (MQ=255)                                                                                         ccTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTc                                                                                   >  1:27092/1‑90 (MQ=255)                                                                                         ccTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGa                                  <  1:1227235/139‑1 (MQ=255)                                                                                           tGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGAt                                                                <  1:1354891/107‑1 (MQ=255)                                                                                           tGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGAt                                                                >  2:1354891/1‑107 (MQ=255)                                                                                                tACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGgcg                                                                                                                   >  1:1337261/1‑51 (MQ=255)                                                                                                tACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGgcg                                                                                                                   <  2:1337261/51‑1 (MQ=255)                                                                                                       tgatgCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGACGCGGTGGATTTgt                    >  2:832021/1‑139 (MQ=255)                                                                                                          tgCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAg                                              <  2:803428/110‑1 (MQ=255)                                                                                                          tgCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAg                                              >  1:803428/1‑110 (MQ=255)                                                                                                          tgCGCAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGACGCGGTGGATTTgttcg                 >  1:568511/1‑139 (MQ=255)                                                                                                            cgcAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGc                                                                                      >  1:1075181/1‑68 (MQ=255)                                                                                                            cgcAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGc                                                                                      <  2:1075181/68‑1 (MQ=255)                                                                                                             gcAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATc                                                                                                      >  1:753628/1‑51 (MQ=255)                                                                                                             gcAAACGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATc                                                                                                      <  2:753628/51‑1 (MQ=255)                                                                                                                 aCGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTc                                                                                   <  2:619989/66‑1 (MQ=255)                                                                                                                 aCGCCTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTc                                                                                   >  1:619989/1‑66 (MQ=255)                                                                                                                     cTGATGAAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGACGCGGTGGATTTGTTCGTTCTTCCTGGc      <  2:534178/139‑1 (MQ=255)                                                                                                                       gatgaAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGACGCGGTGGATTTGTTCGTTCTTCCTGGCag    <  1:306641/139‑1 (MQ=255)                                                                                                                        atgaAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAg                                           >  1:476370/1‑99 (MQ=255)                                                                                                                        atgaAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAg                                           <  2:476370/99‑1 (MQ=255)                                                                                                                        atgaAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGACGCGGTGGATTTGTTCGTTCTTCCTGGCagg   >  1:651958/1‑139 (MQ=255)                                                                                                                         tgaAAATGTGGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGACGCGGTGGATTTGTTCGTTCTTCCTGGCagga  <  1:832021/139‑1 (MQ=255)                                                                                                                                 tgGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTc                                                                                   <  1:1470651/50‑1 (MQ=255)                                                                                                                                 tgGCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTc                                                                                   >  2:1470651/1‑50 (MQ=255)                                                                                                                                  ggCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCtt                                                                                    <  2:532902/48‑1 (MQ=255)                                                                                                                                  ggCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCtt                                                                                    >  1:532902/1‑48 (MQ=255)                                                                                                                                   gCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCtt                                                                                    >  1:812636/1‑47 (MQ=255)                                                                                                                                   gCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCtt                                                                                    <  1:1606293/47‑1 (MQ=255)                                                                                                                                   gCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCtt                                                                                    >  2:1606293/1‑47 (MQ=255)                                                                                                                                   gCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCtt                                                                                    <  2:812636/47‑1 (MQ=255)                                                                                                                                   gCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGAcgcg                              <  2:714604/101‑1 (MQ=255)                                                                                                                                   gCGTGAGGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGAcgcg                              >  1:714604/1‑101 (MQ=255)                                                                                                                                         |                                                                                                                           CGCGTACCCGTCCGGATCTGCCAAACCATGAACCTCGCGGCTGCCCGCGCGGTGCCAGCTACTCCTGGTATCTTTACAGTGCCAACCGCCTGAAATACCCGATGATGCGCAAACGCCTGATGAAAATGTGGCGTGAAGCGAAGGCGCTGCATAGCGATCCGGTTGAGGCATGGGCTTCTATCATTGAAGACGCCGATAAAGCGAAAAGCTTTAAGCAGGCGCGTGGACGCGGTGGATTTGTTCGTTCTTCCTGGCAGGA  >  NC_000913/1280099‑1280357

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