Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,539,241	A→G	V1117A (GTC→GCC)	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,539,241	0	A	G	100.0%	156.0 / NA	44	V1117A (GTC→GCC)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/23);  total (21/23)

GATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTG  >  NC_000913/1539104‑1539368                                                                                                                                          |                                                                                                                                gATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCa                                                                                                                               >  2:558905/1‑140 (MQ=255)   tACATCATGGTCTTGCCCGGCGTTACACGTTGGCTGACCCCCGCGCGGGCGGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAAt                                                                                                                             <  2:653617/140‑1 (MQ=255)        catGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTc                                                                                                                        <  1:815341/140‑1 (MQ=255)        catGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTc                                                                                                                        >  1:738123/1‑140 (MQ=255)          tGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCa                                                                                                                      <  2:21876/140‑1 (MQ=255)          tGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCa                                                                                                                      >  1:657791/1‑140 (MQ=255)           ggTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTc                                                                                                                        >  2:235993/1‑137 (MQ=255)           ggTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTc                                                                                                                        <  1:235993/137‑1 (MQ=255)            gTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCAcc                                                                                                                    >  1:384430/1‑140 (MQ=255)               aTGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACcgcg                                                                                                                 >  1:790430/1‑140 (MQ=255)               aTGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACcgcg                                                                                                                 >  1:617743/1‑140 (MQ=255)                tGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTc                                                                                                                        <  1:704849/132‑1 (MQ=255)                tGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTc                                                                                                                        >  2:704849/1‑132 (MQ=255)                tGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCga                                                                                                                <  1:90161/140‑1 (MQ=255)                tGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCga                                                                                                                <  2:767990/140‑1 (MQ=255)                       cggTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCg                                                                                                         <  1:542394/140‑1 (MQ=255)                       cggTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCg                                                                                                         <  2:617743/140‑1 (MQ=255)                        ggTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGt                                                                                                        <  1:715967/140‑1 (MQ=255)                                 tGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAg                                                                                               >  1:745563/1‑140 (MQ=255)                                 tGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAg                                                                                               >  2:531982/1‑140 (MQ=255)                                        ccaccGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGTAGGttt                                                                                        <  1:558905/140‑1 (MQ=255)                                        ccaccGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGttt                                                                                        <  2:169143/140‑1 (MQ=255)                                        ccaccGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGttt                                                                                        <  1:205660/140‑1 (MQ=255)                                        ccaccGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGttt                                                                                        <  2:790430/140‑1 (MQ=255)                                             gcgcgGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACt                                                                                   <  2:384430/140‑1 (MQ=255)                                             gcgcgGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACt                                                                                   <  1:531982/140‑1 (MQ=255)                                              cgcgGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCAc                                                                                    <  2:803438/138‑1 (MQ=255)                                              cgcgGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCAc                                                                                    >  1:803438/1‑138 (MQ=255)                                               gcgGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGt                                                                                 <  2:809474/140‑1 (MQ=255)                                                       tCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAg                                                                         >  1:277361/1‑140 (MQ=255)                                                           cgcgcCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGaa                                                                     <  1:569470/140‑1 (MQ=255)                                                             cgcCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATg                                                                   >  1:163812/1‑140 (MQ=255)                                                                           aaTACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTg                                                     <  2:657791/140‑1 (MQ=255)                                                                            aTACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTgc                                                    >  2:352622/1‑140 (MQ=255)                                                                                  ccACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTg                                                     >  2:852758/1‑133 (MQ=255)                                                                                  ccACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTg                                                     <  1:852758/133‑1 (MQ=255)                                                                                     cccAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAgg                                           <  2:163812/140‑1 (MQ=255)                                                                                          tCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGtt                                      >  1:170183/1‑140 (MQ=255)                                                                                            gTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTaa                                    <  2:738123/140‑1 (MQ=255)                                                                                                    cAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCtttt                            >  2:440009/1‑140 (MQ=255)                                                                                                          tCAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCggg                      >  2:299749/1‑140 (MQ=255)                                                                                                          gcAGTTCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCggg                      >  2:687303/2‑140 (MQ=255)                                                                                                              ttCTCGGGCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGc                  >  2:803396/1‑140 (MQ=255)                                                                                                                   gggCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGttt             >  2:625368/1‑140 (MQ=255)                                                                                                                   gggCATCTGTTTCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGttt             >  2:206538/1‑140 (MQ=255)                                                                                                                              tCGCTGATCCAGGCAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTg  <  1:206538/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                GATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATGCCCCATTTCTGGTGCGGCGTCAGGAAGTTAAGTGCTTTTTCCGGGAAGCCGTTTGGCGGTATCTG  >  NC_000913/1539104‑1539368

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