Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,592,969	C→T	A438V (GCC→GTC)	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,969	0	C	T	100.0%	70.8 / NA	22	A438V (GCC→GTC)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/16);  total (6/16)

TGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATG  >  NC_000913/4592844‑4593096                                                                                                                              |                                                                                                                                tGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATCTTGCCTCTACGCTGGTAGAAAACGTCGGGGaa                                                                                                                          >  2:258071/1‑133 (MQ=255)                               ggCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTAAGCGCCGTGAcc                                                                                            >  1:73602/1‑133 (MQ=255)                                  ccAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCcgc                                                                                         <  1:359719/133‑1 (MQ=255)                                    agcgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTGAGCGCCGTGACCCGCGt                                                                                       <  1:145327/133‑1 (MQ=255)                                    agcgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGt                                                                                       <  1:109057/133‑1 (MQ=255)                                     gcgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTg                                                                                      <  2:378709/133‑1 (MQ=255)                                      cgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTAAGCGCCGTGACCCGCGTGa                                                                                     <  2:73602/133‑1 (MQ=255)                                          gcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGa                                                                                 >  1:416999/1‑133 (MQ=255)                                                aGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTAt                                                                           <  2:17284/133‑1 (MQ=255)                                                             aaTTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATtgct                                                              <  2:397849/133‑1 (MQ=255)                                                                 aaaaGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGGCGTGACCCACGTAACGGGCATTATGAGCGAGATtgcttctg                                                          >  2:297205/1‑133 (MQ=255)                                                                    aGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATtgcttctgctt                                                       <  1:18377/133‑1 (MQ=255)                                                                       cTGATTGAAGCCTCGGTGGGGAAATTGGCTGTTGGCTCTACACTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGcttctgcttctg                                                    <  2:248870/133‑1 (MQ=255)                                                                       cTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGcttctgcttctg                                                    <  2:416999/133‑1 (MQ=255)                                                                            tGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAg                                               <  2:462024/133‑1 (MQ=255)                                                                              aaGACTCGATGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCa                                             <  1:241359/133‑1 (MQ=255)                                                                                     ggTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGt                                      >  2:462952/1‑133 (MQ=255)                                                                                                    tGTTGGCTCTACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGtt                       <  1:610192/133‑1 (MQ=255)                                                                                                           tctACGCTGGTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAg                <  1:16056/133‑1 (MQ=255)                                                                                                                 cTGGTCGATAGCGTCTGGTAAAGAATGGCGTAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTg          <  1:258071/133‑1 (MQ=255)                                                                                                                    gTCGAAAGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTg       >  2:575014/1‑133 (MQ=255)                                                                                                                         aaGCGTCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATg  <  2:362832/133‑1 (MQ=255)                                                                                                                              |                                                                                                                                TGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATG  >  NC_000913/4592844‑4593096

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: ‑