Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,591,974	G→A	A106A (GCG→GCA)	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,591,974	0	G	A	100.0%	70.6 / NA	24	A106A (GCG→GCA)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/12);  total (12/12)

GGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAAT  >  NC_000913/4591732‑4592214                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 ggCGGTCTGTGCTTTAATGCCGTAAAGAATGACAAAGAAAATGTCACTGTTTTACAAACCAGTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGATGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCGGAGCTGATGCAGAGTGCCAGTATGTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGa                                                                                                                                                                                                                                          <  2:72650/251‑1 (MQ=255)                aaTGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTg                                                                                                                                                                                                                            >  1:121145/1‑250 (MQ=255)                                                gttgttccaacccATTCGCCAGCAGCAATCCCCGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATGGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATGACGAAGCGTTGCCTCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCagag                                                                                                                                                                                            <  2:70799/239‑1 (MQ=255)                                                                cGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGTCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACg                                                                                                                                                                            >  1:174832/1‑250 (MQ=255)                                                                 gCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGa                                                                                                                                                                           >  2:60450/1‑250 (MQ=255)                                                                            tCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGGCCCGCGTCAGAGAACCGCCGCAGCGGCAGAGATAAAACGCAATACCGATATTTATcaa                                                                                                                                                                 >  2:231897/1‑248 (MQ=255)                                                                                    gAATGGCAGCGGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGTCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTg                                                                                                                                                        <  1:79590/250‑1 (MQ=255)                                                                                              tGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGCATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGCAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAAtt                                                                                                                                                                               >  1:146777/1‑217 (MQ=255)                                                                                              tGGGTCGCGTTGTGGCAGACGCGTAACCCCCTCAACCGCGCGGGTATCCGCTACATGATGGATAAGAATAATATTGGTAGCGGTTCAACCGTGGCTGAGCTGATGCAGAGTGCCAGTATGTCGCTGAAACAGGCGGAAAAAAACGGGGCAGATTACGAGGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAAtt                                                                                                                                                                               <  2:146777/217‑1 (MQ=255)                                                                                                         gttgCAGACGTGTAACACCCTAAACCGCGCGGGTATCCGCTACCTGATGGCTCAGAATAATATTGGCAGCGGTTCCACCGTTGCTGAGCTGCTGCAGCGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGAccc                                                                                                                                                                                                                        >  2:41312/1‑165 (MQ=255)                                                                                                         gttgCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCGGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGAccc                                                                                                                                                                                                                        <  1:41312/165‑1 (MQ=255)                                                                                                                                            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)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 GGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAAT  >  NC_000913/4591732‑4592214

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