Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,592,948	G→C	G431A (GGC→GCC)	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,948	0	G	C	100.0%	173.3 / NA	49	G431A (GGC→GCC)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (26/23);  total (26/23)

CGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCG  >  NC_000913/4592816‑4593084                                                                                                                                     |                                                                                                                                         cgcgTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACg                                                                                                                                   <  1:1157940/140‑1 (MQ=255) cgcgTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACg                                                                                                                                   <  1:1034025/140‑1 (MQ=255)          gTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGaa                                                                                                                          <  2:433351/140‑1 (MQ=255)             aGCAAGGGCGCGGTTTTGCGG‑GGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAgcg                                                                                                                      <  2:477052/140‑1 (MQ=255)                   ggCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCCGGTCGAAAGCGCCggg                                                                                                                 <  2:320960/140‑1 (MQ=255)                       cgGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGaaa                                                                                                             >  2:690954/1‑140 (MQ=255)                                    gtCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGAt                                                                                                >  2:762762/1‑140 (MQ=255)                                    gtCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGAt                                                                                                >  2:611123/1‑140 (MQ=255)                                    gtCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGAt                                                                                                <  1:1182649/140‑1 (MQ=255)                                    gtCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGAt                                                                                                <  1:690954/140‑1 (MQ=255)                                    gtCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGAt                                                                                                >  1:739753/1‑140 (MQ=255)                                      cgcgGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTg                                                                                              >  2:1231212/1‑140 (MQ=255)                                                     tAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGAccc                                                                               <  2:805725/140‑1 (MQ=255)                                                           ggCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGa                                                                         <  2:739753/140‑1 (MQ=255)                                                           ggCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGa                                                                         <  1:762762/140‑1 (MQ=255)                                                            gCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGAc                                                                        >  2:544008/1‑140 (MQ=255)                                                             cccAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACg                                                                       >  2:217880/1‑140 (MQ=255)                                                                agcgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGAc                                                                    >  2:355939/1‑140 (MQ=255)                                                                agcgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGAc                                                                    >  1:370231/1‑140 (MQ=255)                                                                  cgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACAt                                                                  >  1:714013/1‑140 (MQ=255)                                                                  cgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACAt                                                                  >  1:216566/1‑140 (MQ=255)                                                                  cgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACAt                                                                  >  1:385090/1‑140 (MQ=255)                                                                  cgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACAt                                                                  >  2:512790/1‑140 (MQ=255)                                                                  cgcagcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACAt                                                                  >  2:652750/1‑140 (MQ=255)                                                                      gcgcCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATg                                                              >  2:692890/1‑140 (MQ=255)                                                                         cccAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGc                                                           <  1:286480/140‑1 (MQ=255)                                                                                  cTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATtgct                                                  <  2:134377/140‑1 (MQ=255)                                                                                    cGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATtgcttc                                                <  2:385090/140‑1 (MQ=255)                                                                                        aaaTTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATtgcttctgct                                            >  1:1187982/1‑140 (MQ=255)                                                                                            tAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGACTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGcttctgcttctg                                        <  2:919425/140‑1 (MQ=255)                                                                                              aaaGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCtgat                                      <  2:1207995/140‑1 (MQ=255)                                                                                                aGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCtgatga                                    >  2:810489/1‑140 (MQ=255)                                                                                                      aTTGAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCagag                              >  2:547426/1‑140 (MQ=255)                                                                                                         gAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCg                           <  1:355939/140‑1 (MQ=255)                                                                                                         gAAGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCg                           <  2:640351/140‑1 (MQ=255)                                                                                                          aaGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGt                          >  2:968598/1‑140 (MQ=255)                                                                                                           aGACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTg                         <  1:217880/140‑1 (MQ=255)                                                                                                            gACTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTgg                        <  2:370231/140‑1 (MQ=255)                                                                                                             aCTCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGt                       <  2:590517/140‑1 (MQ=255)                                                                                                               tCGGTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGATATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTAt                     >  2:849556/1‑140 (MQ=255)                                                                                                                cGGTGGGGAGAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATc                    >  2:534601/1‑140 (MQ=255)                                                                                                                 ggTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCg                   >  2:99845/1‑140 (MQ=255)                                                                                                                 ggTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCg                   >  2:97864/1‑140 (MQ=255)                                                                                                                 ggTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCg                   >  1:28300/1‑140 (MQ=255)                                                                                                                  gTGGGGAAAGTGGATGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGa                  <  1:836777/140‑1 (MQ=255)                                                                                                                                tGTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAg    >  2:545517/1‑140 (MQ=255)                                                                                                                                 gTTGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGc   >  1:632560/1‑140 (MQ=255)                                                                                                                                 cttGCCTCTACGCTGGTCGAAAGCGCGGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGc   >  2:202365/2‑140 (MQ=255)                                                                                                                                  ttGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCg  <  1:692890/140‑1 (MQ=255)                                                                                                                                  ttGCCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCg  <  1:512790/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                         CGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCG  >  NC_000913/4592816‑4593084

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