Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,591,444	C→T	intergenic (‑165/‑213)	btsT ← / → tsr	pyruvate:H(+) symporter/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,444	0	C	T	100.0%	179.8 / NA	53	intergenic (‑165/‑213)	btsT/tsr	pyruvate:H(+) symporter/methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/32);  total (21/32)

AAGCCGTTCTAGGTTTGCGTTGCGTTTGAGGAGGTAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAA  >  NC_000913/4591316‑4591568                                                                                                                                 |                                                                                                                             aaGCCGTTCTAGGTTTGCGTTGCGTTTGAGGAGGTAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGAAAATAACTAAGTGGTTATTTTAGGGTTGTAAGCGGTTAGtt                                                                                                                          >  2:161027/1‑133 (MQ=255)    ccGTTCTAGGTTTGCGTTGCGTTTGAGGAGGTAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAGTAACTGAGTGGTAATTTTAGGGATGGAAGCGGTTAGTTTTg                                                                                                                       >  2:136611/1‑133 (MQ=255)                  cgttgcgttTGAGGAGGTAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTGTTTTAGGGATGTCAGCGGTTAGTTTTGCGGCTGAGCGGCaa                                                                                                         >  2:164385/1‑133 (MQ=255)                        gttTGAGGAGGTAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGcgt                                                                                                   >  2:496193/1‑133 (MQ=255)                                aggTAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAg                                                                                           <  2:513222/133‑1 (MQ=255)                                 ggTAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTTCGGTTGGGCGTCAAAGGCGTCGTTTAAGt                                                                                          <  2:81610/133‑1 (MQ=255)                                 ggTAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGt                                                                                          <  2:126781/133‑1 (MQ=255)                                   tAAATTGACCTCTATACTGTAGTGATGTCGCCCGAGTTACGTGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGa                                                                                        <  2:136364/133‑1 (MQ=255)                                   tAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGa                                                                                        <  1:489250/133‑1 (MQ=255)                                   tAAATTGACCGCTATCCTGTAGTGATGTGGCCGGAGTTACGGGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGa                                                                                        <  2:154771/133‑1 (MQ=255)                                     aaTTGACCGCTATCCTGTAGTGATCTTGTCAGAGTTACGCGTGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACg                                                                                      <  2:154570/133‑1 (MQ=255)                                     aaTTGACCCCTATCCAGTTGTGGTGTTGCCTGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACg                                                                                      <  1:164385/133‑1 (MQ=255)                                     aaTTGACCACTATCCCGTAGTGGTGTTGCCTGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACg                                                                                      <  1:183200/133‑1 (MQ=255)                                          aCCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATaaa                                                                                 >  1:493942/1‑133 (MQ=255)                                           ccGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAat                                                                                <  2:272857/133‑1 (MQ=255)                                            cGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGGGTCGTTGAAGTGACGATAAata                                                                               >  1:157068/1‑133 (MQ=255)                                               tATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATAtgt                                                                            >  2:436666/1‑133 (MQ=255)                                                     gTAGTGATGTTGCCACAGTTACTCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCa                                                                      <  1:159234/133‑1 (MQ=255)                                                        gTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCatat                                                                   <  2:404385/133‑1 (MQ=255)                                                        gTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCatat                                                                   <  2:53114/133‑1 (MQ=255)                                                         tGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATc                                                                  >  1:440312/1‑133 (MQ=255)                                                          gATGTTGCCAGAGTTTCGCGCGAGTAAAGTCAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATca                                                                 <  1:213601/133‑1 (MQ=255)                                                          gATGTTGCCAGAGTTCCGAGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATca                                                                 <  1:161381/133‑1 (MQ=255)                                                          gATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATca                                                                 <  1:436666/133‑1 (MQ=255)                                                          gATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATca                                                                 <  2:493942/133‑1 (MQ=255)                                                           aTGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTGAAGTGACGATAGATGTGTGATACATATcac                                                                >  1:149724/1‑133 (MQ=255)                                                           aTGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTAAAGTGACTATATATGTGTGATGCATATcac                                                                >  1:211622/1‑133 (MQ=255)                                                            tGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATcaca                                                               <  1:277511/133‑1 (MQ=255)                                                            tGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATcaca                                                               <  2:321380/133‑1 (MQ=255)                                                            tGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATcaca                                                               <  1:346325/133‑1 (MQ=255)                                                            tGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTAATATcaca                                                               <  2:111613/133‑1 (MQ=255)                                                            tGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGCTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGGGACGATAAATATGTGATTCATATcaca                                                               <  2:400930/133‑1 (MQ=255)                                                              ttGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGACGTCGTTTCAGTGACGGTAAATATGGGATTGATATCACata                                                             >  2:135584/1‑133 (MQ=255)                                                                gCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGGTCAATATGTGATTCAAATCACATAtt                                                           >  2:235417/1‑133 (MQ=255)                                                                gCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATAtt                                                           >  1:272999/1‑133 (MQ=255)                                                                    gagTTACGTGCGCGTTAAGTAAGGTAAATAGCTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTtata                                                       <  1:235417/133‑1 (MQ=255)                                                                     agTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTatat                                                      <  2:409269/133‑1 (MQ=255)                                                                               gcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGaataat                                            >  1:310590/1‑133 (MQ=255)                                                                               gcgTAAAGAAAGGTACATCACTAAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGaataat                                            <  1:161027/133‑1 (MQ=255)                                                                                cgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAAtt                                           >  1:84960/1‑133 (MQ=255)                                                                                cgTAAAGTAAGGTAAATAACTGAGTGGGTATTTTAGGGCTGTAAGCGGTTAGTTTTGCGGCTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCGCATGTTAATATTGGGAATAAtt                                           >  1:152191/1‑133 (MQ=255)                                                                                   aaaGTAAGTTAAATATCTCAGTCGTTATTTCAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTAt                                        <  1:136611/133‑1 (MQ=255)                                                                                   aaaGTAAGGTAAATAACTGAGTGGTTATTTTATGGATGTAATCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTAt                                        <  2:216307/133‑1 (MQ=255)                                                                                    aaGTAAGGTAAATAAGTGATTGGTTATTTTAGCGATGTAAGGGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATg                                       <  2:211622/133‑1 (MQ=255)                                                                                      gTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTa                                     <  2:497482/133‑1 (MQ=255)                                                                                               aTAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGc                            >  1:6380/1‑133 (MQ=255)                                                                                               aTAACTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGc                            <  2:6380/133‑1 (MQ=255)                                                                                                   cTGAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCtttt                        <  2:84960/133‑1 (MQ=255)                                                                                                     gAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTa                      >  2:23174/1‑133 (MQ=255)                                                                                                     gAGTGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTa                      >  2:66591/1‑133 (MQ=255)                                                                                                        tGGTTATTTTAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAt                   >  1:11765/1‑133 (MQ=255)                                                                                                           ttATTTTAGGGATGTGAGCGGTTAGTTTTGCGTTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                <  2:157068/133‑1 (MQ=255)                                                                                                               tttAGGGATGTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAgt            <  1:23174/133‑1 (MQ=255)                                                                                                                        gTAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAGTAATGGCGTTTTAAATACAGATTGTGAATaaa   >  2:208576/1‑133 (MQ=255)                                                                                                                         tAAGCGGTTAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATaaaa  <  2:440312/133‑1 (MQ=255)                                                                                                                                 |                                                                                                                             AAGCCGTTCTAGGTTTGCGTTGCGTTTGAGGAGGTAAATTGACCGCTATCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAA  >  NC_000913/4591316‑4591568

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