Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,593,116	A→G	N487S (AAC→AGC)	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,593,116	0	A	G	100.0%	105.1 / NA	32	N487S (AAC→AGC)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/16);  total (16/16)

GCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAACAGCAGCGTGAAACATCGGCTGTGGTAAAAACCGTGAC  >  NC_000913/4592983‑4593251                                                                                                                                      |                                                                                                                                        gCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGcgccgc                                                                                                                                   <  1:1197227/140‑1 (MQ=255)   ggAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCgcgc                                                                                                                                 <  2:1124457/140‑1 (MQ=255)           tCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGcg                                                                                                                                       <  1:142615/126‑1 (MQ=255)           tCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGcg                                                                                                                                       >  2:142615/1‑126 (MQ=255)                 ccGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTct                                                                                                                   >  1:388815/1‑140 (MQ=255)                               aCGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTgccgccgc                                                                                                     <  2:1177530/140‑1 (MQ=255)                                       tATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTgg                                                                                             <  2:991872/140‑1 (MQ=255)                                            gCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAagag                                                                                        <  1:1531865/140‑1 (MQ=255)                                               aaaTTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAg                                                                                     >  1:1187584/1‑140 (MQ=255)                                               aaaTTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAg                                                                                     >  2:995927/1‑140 (MQ=255)                                                              ctgATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGAcc                                                                      <  1:1471734/140‑1 (MQ=255)                                                               tgatgaGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCg                                                                     <  2:1187584/140‑1 (MQ=255)                                                                gatgaGCAGAGCCGTGAGATCGATCAGGTTGTCTTAGCGCTTGCTGAGTTGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGa                                                                    <  2:681409/140‑1 (MQ=255)                                                                  tgaGCAGAGCCGTGGTACCGATAAGGTTGGCTTAGCGGTCGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAg                                                                  <  2:551858/140‑1 (MQ=255)                                                                      cAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAgcagt                                                              >  2:1469373/1‑140 (MQ=255)                                                                        gagCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAgcagtgg                                                            >  2:678659/1‑140 (MQ=255)                                                                         agCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAgcagtggc                                                           >  2:761261/1‑140 (MQ=255)                                                                          gCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAgcagtggca                                                          <  2:982958/140‑1 (MQ=255)                                                                           ccGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGGTGGGAGAGGAGGCCAGGAGCCGGACCGAAggagtggcag                                                         >  2:738105/1‑140 (MQ=255)                                                                                   tCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCgg                                                 >  2:964543/1‑140 (MQ=255)                                                                                      aTCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTAGAAGAGCAAGCCGGTCACCTGACCGCAGCAATGGCAGTGTTCCGGAtt                                              >  2:472973/1‑140 (MQ=255)                                                                                         aGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGAt                                               <  1:1178152/136‑1 (MQ=255)                                                                                         aGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGAt                                               >  2:1178152/1‑136 (MQ=255)                                                                                               gCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAAcag                                     >  2:95371/1‑140 (MQ=255)                                                                                                cTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAGTCAGTGGCAGTGTTCCGGATTCAGCAAcagc                                    >  2:751235/1‑140 (MQ=255)                                                                                                cTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGCCCAGTCGCCTGATCAAGGCAGTGGCTGTGTGCCGGATTCAGCAAcagc                                    >  2:683203/1‑140 (MQ=255)                                                                                                   aGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAACagcagc                                 >  2:931281/1‑140 (MQ=255)                                                                                                    gCGGTTGCTGAGATGGATCGGGGAACTCAACCGAGGGCCACGCTGGTGTAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAACAGCAGCg                                <  2:653801/140‑1 (MQ=255)                                                                                                               gaTGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAACAGCAGCGTGAAACATCgg                     <  2:388815/140‑1 (MQ=255)                                                                                                               gaTGGACCGGGTAACTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAACAGCAGCGTGAAACATCgg                     <  1:931281/140‑1 (MQ=255)                                                                                                                 tGGACCGGGTAACTCTACAGAGCGCCGCGCTGGTGGACGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAACAGCAGCGTGAAACATCGGCt                   <  1:472973/140‑1 (MQ=255)                                                                                                                          tAACTCAACAGAGCGTCGCGCTAGTCGAATAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGAACGAAGCAGTGGCAGTGTTCCGGATTCAGCAACAGCAGCGTGAAACATCGGCTGTGGTaaaa          <  1:722298/140‑1 (MQ=255)                                                                                                                             cTCAACAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAACAGCAGCGTGAAACATCGGCTGTGGTAAAAAcc       >  2:1115073/1‑140 (MQ=255)                                                                                                                                  cAGAGCGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAACAGCAGCGTGAAACATCGGCTGTGGTAAAAACCGTGAc  >  2:1221417/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                        GCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAGCAGTGGCAGTGTTCCGGATTCAGCAACAGCAGCGTGAAACATCGGCTGTGGTAAAAACCGTGAC  >  NC_000913/4592983‑4593251

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