Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,542,454	A→G	V46A (GTG→GCG)	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,542,454	0	A	G	100.0%	155.8 / NA	45	V46A (GTG→GCG)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/20);  total (25/20)

GGGCAGGTCAGGGCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGT  >  NC_000913/1542333‑1542582                                                                                                                          |                                                                                                                                 gggCAGGTCAGGGCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGcc                                                                                                                >  1:84748/1‑140 (MQ=255)  ggCAGGTCAGGGCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCa                                                                                                               >  2:242010/1‑140 (MQ=255)  ggCAGGTCAGGGCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCa                                                                                                               >  2:434919/1‑140 (MQ=255)     aGGTCAGGGCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACg                                                                                                            <  1:197434/140‑1 (MQ=255)       gTCAGGGCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCt                                                                                                          >  1:141514/1‑140 (MQ=255)       gTCAGGGCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCt                                                                                                          <  2:275581/140‑1 (MQ=255)            ggCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGa                                                                                                     >  1:324433/1‑140 (MQ=255)            ggCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGa                                                                                                     >  1:533170/1‑140 (MQ=255)                aGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCCGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGc                                                                                                 <  2:635656/140‑1 (MQ=255)                 gTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCt                                                                                                >  1:366001/1‑140 (MQ=255)                                 tCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCt                                                                                                >  1:523608/1‑124 (MQ=255)                                 tCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCt                                                                                                <  2:523608/124‑1 (MQ=255)                                  cTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGt                                                                               <  2:324433/140‑1 (MQ=255)                                  cTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGt                                                                               <  2:533170/140‑1 (MQ=255)                                     ttGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTCGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGc                                                                            <  2:271862/140‑1 (MQ=255)                                       ggATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGt                                                                               <  2:759620/135‑1 (MQ=255)                                       ggATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGt                                                                               >  1:759620/1‑135 (MQ=255)                                       ggATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTa                                                                          >  2:563801/1‑140 (MQ=255)                                       ggATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTa                                                                          <  1:434919/140‑1 (MQ=255)                                             cccAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATc                                                                    >  2:624421/1‑140 (MQ=255)                                               cAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAGCGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCAc                                                                  >  2:576297/1‑140 (MQ=255)                                                aGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCAcc                                                                 <  1:284502/140‑1 (MQ=255)                                                aGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCAcc                                                                 <  1:372784/140‑1 (MQ=255)                                                    cACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTc                                                             >  2:14217/1‑140 (MQ=255)                                                            cATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCa                                                     >  1:733256/1‑140 (MQ=255)                                                                  tAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCa                                               <  2:126653/140‑1 (MQ=255)                                                                      gTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAgg                                           >  1:772230/1‑140 (MQ=255)                                                                       tAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGt                                          >  2:283110/1‑140 (MQ=255)                                                                         gATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGttt                                        >  1:188301/1‑140 (MQ=255)                                                                         gATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGttt                                        >  2:383367/1‑140 (MQ=255)                                                                                cAGCTACAGGAGTCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCtt                                 <  1:738753/140‑1 (MQ=255)                                                                                cAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCtt                                 <  1:576297/140‑1 (MQ=255)                                                                                cAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCtt                                 <  2:141514/140‑1 (MQ=255)                                                                                    tACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTtttgt                             <  1:189912/140‑1 (MQ=255)                                                                                        ggAGCCTGTACAGTTAACACCGTGGGTGGAACCCGCGATTCTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGa                         <  2:733256/140‑1 (MQ=255)                                                                                            ccTGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTa                     >  1:829114/1‑140 (MQ=255)                                                                                              tGTACAGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGc                   >  2:781702/1‑140 (MQ=255)                                                                                                   aGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAg              <  2:84748/140‑1 (MQ=255)                                                                                                   aGTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAg              <  2:188301/140‑1 (MQ=255)                                                                                                    gTTAACACCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGc             >  1:553870/1‑140 (MQ=255)                                                                                                          acCGTGGGTGGAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCATTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCa       >  1:175825/1‑140 (MQ=255)                                                                                                             gtgggtggAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCaaaa    >  2:476672/1‑140 (MQ=255)                                                                                                             gtgggtggAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCaaaa    >  1:679593/1‑140 (MQ=255)                                                                                                               gggtggAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGt  <  2:389833/140‑1 (MQ=255)                                                                                                               gggtggAACGCGCGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGt  >  2:296503/1‑140 (MQ=255)                                                                                                                          |                                                                                                                                 GGGCAGGTCAGGGCGAGTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGT  >  NC_000913/1542333‑1542582

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