Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,542,486	G→A	S35S (AGC→AGT)	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,486	0	G	A	100.0%	141.0 / NA	41	S35S (AGC→AGT)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/18);  total (23/18)

GTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGG  >  NC_000913/1542349‑1542620                                                                                                                                          |                                                                                                                                       gTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACt                                                                                                                                      >  1:512337/1‑140 (MQ=255)  tGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTg                                                                                                                                     >  2:438456/1‑140 (MQ=255)            gTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGt                                                                                                                           <  1:438456/140‑1 (MQ=255)                 tCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCgg                                                                                                                      <  1:510022/140‑1 (MQ=255)                  cTGTTGGATTTCCCAGGTCACCAGACCATTTTTAATGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGt                                                                                                                     <  2:808520/140‑1 (MQ=255)                             cccAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATc                                                                                                          >  1:668462/1‑140 (MQ=255)                                  gTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCAcc                                                                                                       >  2:726677/1‑138 (MQ=255)                                  gTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCAcc                                                                                                       <  1:726677/138‑1 (MQ=255)                                       cAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGt                                                                                                >  2:426556/1‑140 (MQ=255)                                       cAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGt                                                                                                >  1:977136/1‑140 (MQ=255)                                           ccATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGAcc                                                                                            <  2:657270/140‑1 (MQ=255)                                                        aGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGtt                                                                               >  2:228554/1‑140 (MQ=255)                                                                cAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCtt                                                                       <  2:512337/140‑1 (MQ=255)                                                                   cTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTtttg                                                                    >  1:1140479/1‑140 (MQ=255)                                                                   cTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCCtttg                                                                    >  1:1109597/1‑140 (MQ=255)                                                                   cTACAGGAGCCTGTACAGTTAACACCGTGGGTGG‑ACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTtttgt                                                                   >  2:433245/1‑140 (MQ=255)                                                                      cAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTtttgttt                                                                 <  1:686281/140‑1 (MQ=255)                                                                              tGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGa                                                               >  2:199872/1‑134 (MQ=255)                                                                              tGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGa                                                               <  1:199872/134‑1 (MQ=255)                                                                              tGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGc                                                         >  2:952217/1‑140 (MQ=255)                                                                                  cAGTTAACACCGTGGGTGG‑ACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAg                                                    <  1:433245/140‑1 (MQ=255)                                                                                     ttAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCg                                                  >  1:802825/1‑140 (MQ=255)                                                                                     ttAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCg                                                  >  1:630806/1‑140 (MQ=255)                                                                                     ttAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCg                                                  >  2:365551/1‑140 (MQ=255)                                                                                        acacCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCCTTTGTTTGAAGTAGCGAAAGCGATc                                               <  2:1109597/140‑1 (MQ=255)                                                                                          acCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCa                                             >  2:1174517/1‑140 (MQ=255)                                                                                             gtgggtggAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCaaaa                                          >  2:914715/1‑140 (MQ=255)                                                                                             gtgggtggAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCaaaa                                          >  1:27925/1‑140 (MQ=255)                                                                                             gtgggtggAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCaaaa                                          >  2:961653/1‑140 (MQ=255)                                                                                               gggtggAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGt                                        >  1:224730/1‑140 (MQ=255)                                                                                                          aCGATTTTGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTcatgac                             <  1:225859/140‑1 (MQ=255)                                                                                                                 tGTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTc                      >  2:1136454/1‑140 (MQ=255)                                                                                                                  gTCGAACTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTcc                     <  1:914715/140‑1 (MQ=255)                                                                                                                       aCTGCCAACGCTGGTGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTc                <  2:178740/140‑1 (MQ=255)                                                                                                                        cTGCCAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATAACATGACTCCTGCTcc               <  2:302777/140‑1 (MQ=255)                                                                                                                            cAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGga           <  2:180357/140‑1 (MQ=255)                                                                                                                            cAACGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGga           <  1:228554/140‑1 (MQ=255)                                                                                                                              aCGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGgaat         <  2:224730/140‑1 (MQ=255)                                                                                                                              aCGCTGGCGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGgaat         <  1:961653/140‑1 (MQ=255)                                                                                                                                     cGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGTGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTgg  >  2:992532/1‑140 (MQ=255)                                                                                                                                     cGATAACTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTgg  >  2:260177/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       GTGCGCGGGTAGTCGGTCTGTTGGATTTCCCAGGTCACCAGACCATTTTTAACGTAGATTTTCCAGCTACAGGAGCCTGTACAGTTAACACCGTGGGTGGAACGCACGATTTTGTCGAACTGCCAACGCTGGCGATAGCTGTCCTCCCAGTCGCGGTTGCTATGCATCACCTGTCCGTGACCATCGGCAAAGGTTTCGCCCTTTTGTTTGAAGTAGCGAAAGCGATCCAAAAGTTTACTCATGACATGACTCCTGCTCCAGGAATGAATTGG  >  NC_000913/1542349‑1542620

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