Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,541,275	T→G	N439T (AAT→ACT)	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,541,275	0	T	G	100.0%	122.7 / NA	36	N439T (AAT→ACT)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (18/18);  total (18/18)

AGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAA  >  NC_000913/1541143‑1541411                                                                                                                                     |                                                                                                                                         aGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAAt                                                                                                                                   <  1:284013/140‑1 (MQ=255)        aaTTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGTAGTTTCAATGCCGCCa                                                                                                                            <  1:81118/140‑1 (MQ=255)                 gTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAgggg                                                                                                                   <  1:156727/140‑1 (MQ=255)                 gTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAgggg                                                                                                                   <  2:246788/140‑1 (MQ=255)                  tACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGa                                                                                                                  >  1:217985/1‑140 (MQ=255)                       aaTCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGcc                                                                                                             <  1:16837/140‑1 (MQ=255)                            tAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGcc                                                                                                        >  2:5829/1‑140 (MQ=255)                            tAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGcc                                                                                                        >  2:194414/1‑140 (MQ=255)                                    gACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAg                                                                                                >  1:288813/1‑140 (MQ=255)                                    gACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAg                                                                                                >  2:194595/1‑140 (MQ=255)                                     accaccGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCAGCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGc                                                                                               >  2:566577/1‑140 (MQ=255)                                     accaccGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGc                                                                                               <  1:330965/140‑1 (MQ=255)                                               aGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGa                                                                                     <  1:194595/140‑1 (MQ=255)                                                 gtgtTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGa                                                                                     <  1:159209/138‑1 (MQ=255)                                                 gtgtTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGa                                                                                     >  2:159209/1‑138 (MQ=255)                                                 gtgtTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGAcc                                                                                   >  1:420533/1‑140 (MQ=255)                                                 gtgtTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGAcc                                                                                   <  2:420533/140‑1 (MQ=255)                                                      gCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCa                                                                              >  1:656857/1‑140 (MQ=255)                                                      gCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCa                                                                              <  1:566577/140‑1 (MQ=255)                                                           cGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGtt                                                                         >  1:794390/1‑140 (MQ=255)                                                           cGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGtt                                                                         >  1:614006/1‑140 (MQ=255)                                                           cGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGtt                                                                         >  1:513636/1‑140 (MQ=255)                                                           cGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGtt                                                                         >  2:55900/1‑140 (MQ=255)                                                            gACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGa                                                                                     <  1:290545/127‑1 (MQ=255)                                                            gACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGa                                                                                     >  2:290545/1‑127 (MQ=255)                                                             aCTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTaa                                                                       <  2:288813/140‑1 (MQ=255)                                                               tAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACg                                                                     <  2:481980/140‑1 (MQ=255)                                                                      aGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAAttt                                                              <  2:472958/140‑2 (MQ=255)                                                                        gTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTgg                                                            >  1:134856/1‑140 (MQ=255)                                                                        gTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTgg                                                            >  1:605388/1‑140 (MQ=255)                                                                                                     ccGGATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgattccaga                               <  2:656857/140‑1 (MQ=255)                                                                                                        gATTGTGTTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGgattccagattc                            >  2:525531/1‑140 (MQ=255)                                                                                                             tgtTTTACGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCAttt                       >  1:436458/1‑140 (MQ=255)                                                                                                                    cGCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCtttt                >  2:277230/1‑140 (MQ=255)                                                                                                                     gCTGCGAAAATGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCttttt               <  1:277230/140‑1 (MQ=255)                                                                                                                             aaTGCGGAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAg       <  1:731074/140‑1 (MQ=255)                                                                                                                                  ggAGTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTTGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCTATTTGCCTTTTTCTCCCCAGCGGaa  <  1:5829/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                         AGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAA  >  NC_000913/1541143‑1541411

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