Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,592,082	G→A	L142L (CTG→CTA)	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,082	0	G	A	100.0%	102.5 / NA	32	L142L (CTG→CTA)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/14);  total (18/14)

CGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAAT  >  NC_000913/4591861‑4592325                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     cgcgcgGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCtt                                                                                                                                                                                                                         >  1:166110/1‑251 (MQ=255)   ggcgggttccccctccatgaggaaTCAAAATATTATTGGTGGCGATTCAACCGTTGCGGATCTGATGCG‑GATTCCCAGTATTCCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCTTTGCCTCGTGTCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACGATTAGGTGCAGGCAAAATCAACGAGTTCTTTg                                                                                                                                                                                                                       <  2:249375/228‑1 (MQ=255)    gcgGGTTTCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGGGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGa                                                                                                                                                                                                                      <  1:493244/251‑1 (MQ=255)                                     atTGGTAGCGGTTCACCCGTTGCGGAGCTG‑TGGAAGGGGGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTTGCGTCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCga                                                                                                                                                                                     <  2:135781/251‑1 (MQ=255)                                           aGCGGTTCAACCGTTGCTGAGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGc                                                                                                                                                                                >  1:136955/1‑250 (MQ=255)                                                      cGTTGCTGGGCTGATGCA‑GAGTGCCAGTATTTCGCTGAAACGGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTTATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCtt                                                                                                                                                                    <  1:247056/251‑1 (MQ=255)                                                              ggctgaTGCA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAg                                                                                                                                                            <  1:411632/250‑1 (MQ=255)                                                                       a‑gagTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGAt                                                                                                                                                   <  1:131974/251‑1 (MQ=255)                                                                            tGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGc                                                                                                                                               >  1:36683/1‑251 (MQ=255)                                                                                     ttCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGat                                                                                                                                       >  1:176451/1‑250 (MQ=255)                                                                                        gCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCg                                                                                                                                   <  1:134409/251‑1 (MQ=255)                                                                                                 ggCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTcta                                                                                                                                           <  1:309714/234‑3 (MQ=255)                                                                                                 ggCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGGGGCTTACCTGGGGCAAAACGATCGGCTctc                                                                                                                                           >  2:309714/1‑232 (MQ=255)                                                                                                    ggAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGAt                                                                                                                        >  1:450426/1‑250 (MQ=255)                                                                                                         cgtgcggggaagATTACGAAGCGTGGCCGCTTGACCCGCGCCAGATCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAAt                                                                                                                  <  2:400388/240‑1 (MQ=255)                                                                                                               gggcagtttACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGcctc                                                                                                             <  1:449140/243‑1 (MQ=255)                                                                                                               gggCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAAt                                                                                                                  >  1:400388/1‑245 (MQ=255)                                                                                                                                   ccGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCCACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGCAATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGAt                                                                                        >  1:34551/1‑251 (MQ=255)                                                                                                                                                     agCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTgatgat                                                                      <  1:142926/251‑1 (MQ=255)                                                                                                                                                      cccccgccgCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtc                                                                     <  2:450426/248‑1 (MQ=255)                                                                                                                                                         ccgccgCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGACTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCTTGATGAtcgtc                                                                  >  1:342616/1‑251 (MQ=255)                                                                                                                                                             cgcAGCGGCAGAGAACAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTAGGATATCGCCGTAAGCGATAACAATGCCTCCTACAGCCAGGGGGTGGGGGTTCTGGTGGGCGGGCTGATGGTCGTAc                                                              >  2:10024/1‑251 (MQ=255)                                                                                                                                                                      agagATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTAAAGCCAGGCGATTGGGATTCTGGGGGGCGAGATGATCGTCGTACTGGCGGTc                                                      >  2:495881/1‑250 (MQ=255)                                                                                                                                                                         gaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTAc                                                                                                         >  1:177647/1‑196 (MQ=255)                                                                                                                                                                         gaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTAc                                                                                                         <  2:177647/196‑1 (MQ=255)                                                                                                                                                                                      tACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCAGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGAGGATTCTGGTGGGCGAGGTGATCGTCGTACTGGCGGGAATCTTTGCCG‑TCTgggt                                    >  1:241026/1‑249 (MQ=255)                                                                                                                                                                                           tatTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCCTCTTCCCCG‑GCTGGGTCgggt                               >  1:127996/1‑249 (MQ=255)                                                                                                                                                                                            atTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGCTAACAATGCCTCCTACAGCCAGGCGATGTGAATTCTGGTGGGCGTGATGCTCGTCGTACTGGCgggacttttctccggcgcgtgggtggg                               >  2:188815/1‑226 (MQ=255)                                                                                                                                                                                               tATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATATTCTCCG‑GCTGGGTCGGTAttt                            >  1:415850/1‑249 (MQ=255)                                                                                                                                                                                                  cacaATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCGGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGCTTCTGGTGGGCGTGATGATCGTCGTACTTGCGGTCATCTTCGCCGCTCGGGTTGGGT‑TTAAAGc                        >  1:40546/1‑251 (MQ=255)                                                                                                                                                                                                        gcgcTGGCGGAGCTGATCCAACTAATAGGTGCAGGCAAAATCAACGAGGTCTTTGATCAGCCGGCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCGTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCt                  <  2:75260/251‑1 (MQ=255)                                                                                                                                                                                                                        tCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGCATCTTCGCCG‑CCTGGGTTGGTATTAAAGCCTCTCTGGTAGCGcccataagt  >  1:439129/1‑244 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     CGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGA‑GAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCG‑TCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAAT  >  NC_000913/4591861‑4592325

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