Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,592,788	C→T	L378L (CTG→TTG)	tsr →	methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,592,788	0	C	T	100.0%	151.2 / NA	48	L378L (CTG→TTG)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/18);  total (30/18)

CTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGA  >  NC_000913/4592654‑4592924                                                                                                                                       |                                                                                                                                         cTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATCATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCt                                                                                                                                     >  2:443836/1‑140 (MQ=255)  ttCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCtt                                                                                                                                    >  1:428801/1‑140 (MQ=255)       aaaCGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAAc                                                                                                                               >  2:702129/1‑140 (MQ=255)          cGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAAcgcg                                                                                                                            <  2:693819/140‑1 (MQ=255)             cgcAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACgcggcg                                                                                                                         >  2:341508/1‑140 (MQ=255)                          tAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTgcgc                                                                                                            >  1:499617/1‑140 (MQ=255)                                ggTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCgg                                                                                                      <  2:499617/140‑1 (MQ=255)                                         cGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCaa                                                                                              >  1:121163/1‑139 (MQ=255)                                         cGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCaa                                                                                              <  2:121163/139‑1 (MQ=255)                                         cGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCaa                                                                                              <  2:580085/139‑1 (MQ=255)                                         cGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCaa                                                                                              >  1:580085/1‑139 (MQ=255)                                         cGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAg                                                                                             >  2:462473/1‑140 (MQ=255)                                         cGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAg                                                                                             >  1:634511/1‑140 (MQ=255)                                           tggtgCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAggg                                                                                           >  2:771956/1‑140 (MQ=255)                                           tggtgCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAggg                                                                                           >  2:510365/1‑140 (MQ=255)                                              tgCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGgcgc                                                                                        >  1:287634/1‑140 (MQ=255)                                              tgCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGgcgc                                                                                        >  2:744220/1‑140 (MQ=255)                                                cAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCgg                                                                                      >  1:53349/1‑140 (MQ=255)                                                cAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCgg                                                                                      >  2:308371/1‑140 (MQ=255)                                                       tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCg                                                                               >  2:622532/1‑140 (MQ=255)                                                        gcgcGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCgg                                                                              >  2:513832/1‑140 (MQ=255)                                                        gcgcGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAAGGCGCGGTTTTGCgg                                                                              >  2:457297/1‑140 (MQ=255)                                                            gATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCggtggt                                                                          <  1:744220/140‑1 (MQ=255)                                                            gATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCggtggt                                                                          <  2:287634/140‑1 (MQ=255)                                                             atatCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCggtg                                                                            >  2:875207/1‑137 (MQ=255)                                                             atatCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCggtg                                                                            <  1:875207/137‑1 (MQ=255)                                                                   ccaccaGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAATGGCGCGGTTTTGCGGTGGTCGCGGga                                                                   >  2:840589/1‑140 (MQ=255)                                                                     accaGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGATGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGgaga                                                                 <  2:304178/140‑1 (MQ=255)                                                                     accaGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGgaga                                                                 <  1:341508/140‑1 (MQ=255)                                                                      ccaGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGaa                                                                <  2:53349/140‑1 (MQ=255)                                                                          ttCGCAGAAAATCGCCGATATTATAAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGc                                                            >  1:371465/1‑140 (MQ=255)                                                                           tCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAATGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCg                                                           <  1:840589/140‑1 (MQ=255)                                                                           tCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCg                                                           <  1:771956/140‑1 (MQ=255)                                                                           tCGCAGAAAATCGCCGATATTATAAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCg                                                           <  2:371465/140‑1 (MQ=255)                                                                            cGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGTGGTGGTCGCGGGAGAAGGGCGt                                                          >  1:309773/1‑140 (MQ=255)                                                                                           atatTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCcagcgc                                           <  1:510365/140‑1 (MQ=255)                                                                                           atatCATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCcagcgc                                           <  1:443836/140‑1 (MQ=255)                                                                                                     cGTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCatg                                 >  2:737478/1‑138 (MQ=255)                                                                                                      gTAATTGACGGCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAggc                                >  2:80343/1‑140 (MQ=255)                                                                                                               ggCATTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGa                       <  1:622532/140‑1 (MQ=255)                                                                                                                  aTTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAAt                    >  2:242283/1‑140 (MQ=255)                                                                                                                  aTTGCCTTCCAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAAt                    >  1:617253/1‑140 (MQ=255)                                                                                                                           cAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCt           >  1:350897/1‑140 (MQ=255)                                                                                                                           cAGACCAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCt           >  1:761720/1‑140 (MQ=255)                                                                                                                                cAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTaa                 >  1:447244/1‑129 (MQ=255)                                                                                                                                cAATATTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTaa                 <  2:447244/129‑1 (MQ=255)                                                                                                                                   tatTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAg   >  2:537753/1‑140 (MQ=255)                                                                                                                                    atTTTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGa  <  2:738019/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         CTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGA  >  NC_000913/4592654‑4592924

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑