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%	62.0 / NA	22	L142L (CTG→CTA)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/12);  total (10/12)

TCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGC  >  NC_000913/4591871‑4592308                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   tCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCCGGCAAAATCAACGAGTTCTTTGATCAGCCg                                                                                                                                                                                             >  2:75145/1‑251 (MQ=255)   cGCTACATGATGGATCAGAATAATATTGGTAGCGGTGCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGAc                                                                                                                                                                                           <  1:9848/251‑1 (MQ=255)                   aGAATAATATTGGGAGCGGTCCAACCGTTGATGGGCAGAGTAAGAGTGCCAGTATTTCGCGGAAACAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTGGCCGCGTGACCCGGGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGATCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGAACCAGGGATAGCAGGAc                                                                                                                                                                           <  2:174832/251‑1 (MQ=255)                          tatTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGCGTGCCAGTATTTCGCTGAAACAGGCGGACAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCg                                                                                                                                                                    >  1:236951/1‑251 (MQ=255)                                      ttCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAACCAGGCGGAAAAAAACTGGGCAGATTACGACGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCAAACTATTAGGTGCAGACAAAATCAACGAGTTTTTTGATCAGCCGACGCAGGGATATCAGGACGGTTTCGCGAAGCAGtctg                                                                                                                                                        >  1:130615/1‑248 (MQ=255)                                                          ggcagagTGCCAGTATTTCGCTGAAACAGGCGGAAAAAACCTGGGCAGATTACGAAGCGTTGCGGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAAc                                                                                                                                    <  1:38923/248‑1 (MQ=255)                                                                   ccAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGAAGAGATCACACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCATGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTc                                                                                                                           >  1:103796/1‑251 (MQ=255)                                                                    cAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTct                                                                                                                          >  1:52934/1‑250 (MQ=255)                                                                    cAGTATCTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTGGCCGCGTGACCCGCGTCAGATCACCGCCGCAGCGGTAGAGATCAAACGTAATTACGATATTTATCACAAGGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGTCGAACCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCACAACGATCGGCTct                                                                                                                          <  2:150703/251‑2 (MQ=255)                                                                            ggCTGAAACAGGCGGAAAAAAACTGGGCAATTTACGAAGCGTTGCCCCGTTACCCTCGGCAGAGCCCCGCCGCAGCGGCAGAGATCAAACGTAAATACGATATTTATCACAATGCGCTGGCGGATCTGATCCAACTATTAAGTGCAGGCAAAATCAACGAGTTCTTTTATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATc                                                                                                                  <  1:123198/250‑1 (MQ=255)                                                                                gAAACAGGCGGAAAAAAACTGGGCAGATTACGACGCGTTGCCGCGTGACCCGCGTCAGCGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGCTATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGatctctccg                                                                                                              >  2:86401/1‑244 (MQ=255)                                                                                  aaCAGGCGGAAAAAAACTGGGCAGATTCCGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTc                                                                                                            <  1:126957/251‑1 (MQ=255)                                                                                              aaaccggggCAGATTACGAAGCGTTGCCGCGTGACCCGCGTAGGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAAt                                                                                                <  2:121145/245‑1 (MQ=255)                                                                                                     ggcagATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGcctcct                                                                                         >  1:91444/5‑251 (MQ=255)                                                                                                                                                       cGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTGTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCGGGTGGGCGTGATGATCGGCGTACTGGCg                                       <  2:257111/251‑1 (MQ=255)                                                                                                                                                        ggCAGAGATCAAACGTAATTACGCTATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAACCGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCgg                                      >  2:90011/1‑251 (MQ=255)                                                                                                                                                                      gTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATc                                                                                                                                >  2:223223/1‑147 (MQ=255)                                                                                                                                                                      gTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATc                                                                                                                                <  1:223223/147‑1 (MQ=255)                                                                                                                                                                          ttACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCTggg                    >  1:233477/1‑250 (MQ=255)                                                                                                                                                                                     aTCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCATCGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTaa          <  1:75145/250‑1 (MQ=255)                                                                                                                                                                                          aaTGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTc    <  1:60450/251‑1 (MQ=255)                                                                                                                                                                                             gcgcTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGc  <  2:233477/250‑1 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   TCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGC  >  NC_000913/4591871‑4592308

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