Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,285,372	A→G	E78E (GAA→GAG)	narJ →	molybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,285,372	0	A	G	100.0%	139.4 / NA	38	E78E (GAA→GAG)	narJ	molybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/18);  total (20/18)

CAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCG  >  NC_000913/1285245‑1285504                                                                                                                                |                                                                                                                                     cAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACgg                                                                                                                           <  1:1415943/139‑1 (MQ=255) cAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACgg                                                                                                                           <  1:1468087/139‑1 (MQ=255)          aTGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCcg                                                                                                                  >  2:411830/1‑139 (MQ=255)          aTGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCcg                                                                                                                  >  2:178475/1‑139 (MQ=255)               cATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGAcc                                                                                                             <  1:1135669/139‑1 (MQ=255)                     cTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGcc                                                                                                       >  1:591778/1‑139 (MQ=255)                                         tttAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgctg                                                                                   >  2:248675/1‑139 (MQ=255)                                         tttAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgctg                                                                                   >  1:638385/1‑139 (MQ=255)                                           tAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgc                                                                                 <  2:638385/139‑1 (MQ=255)                                                         ccGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcagc                                                                   >  2:1354184/1‑139 (MQ=255)                                                               gATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt                                                             <  1:767481/139‑1 (MQ=255)                                                                 tGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTg                                                           >  1:341134/1‑139 (MQ=255)                                                                 tGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGCTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTg                                                           >  2:873247/1‑139 (MQ=255)                                                                       ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTa                                                     >  2:721213/1‑139 (MQ=255)                                                                              tACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa                                                                                          >  2:339511/1‑95 (MQ=255)                                                                              tACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa                                                                                          <  1:339511/95‑1 (MQ=255)                                                                                       cTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGc                                     <  1:873247/139‑1 (MQ=255)                                                                                        tGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGcc                                    >  1:931484/1‑139 (MQ=255)                                                                                               ccGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGAcc                                                                                                             <  1:768286/59‑1 (MQ=255)                                                                                               ccGTGGCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGAcc                                                                                                             >  2:768286/1‑59 (MQ=255)                                                                                                    gCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGcc                        <  2:1025348/139‑1 (MQ=255)                                                                                                    gCCGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGcc                        >  2:339710/1‑139 (MQ=255)                                                                                                     ccGCGCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCg                       <  2:817036/139‑1 (MQ=255)                                                                                                      cgcgCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGc                      <  2:1409733/139‑1 (MQ=255)                                                                                                      cgcgCCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGc                      <  1:1354184/139‑1 (MQ=255)                                                                                                       gcgcCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCt                     <  1:178475/139‑1 (MQ=255)                                                                                                       gcgcCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCt                     <  2:341134/139‑1 (MQ=255)                                                                                                         gcCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcgcg                                           >  1:569830/1‑115 (MQ=255)                                                                                                         gcCACGTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcgcg                                           <  2:569830/115‑1 (MQ=255)                                                                                                              gTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGc                      <  2:1345446/131‑1 (MQ=255)                                                                                                              gTCACTGTTGCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGc                      >  1:1345446/1‑131 (MQ=255)                                                                                                                       gCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTg     >  2:790078/1‑139 (MQ=255)                                                                                                                       gCTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTg     >  1:224213/1‑139 (MQ=255)                                                                                                                        cTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt                                                             <  1:1305806/82‑1 (MQ=255)                                                                                                                        cTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt                                                             >  2:1305806/1‑82 (MQ=255)                                                                                                                        cTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTg                          >  2:1509755/1‑117 (MQ=255)                                                                                                                        cTGTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTg                          <  1:1509755/117‑1 (MQ=255)                                                                                                                          gTTTGAGCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGgcg  >  2:1187045/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                     CAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCG  >  NC_000913/1285245‑1285504

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