Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,653,266	A→G	N258S (AAC→AGC)	sseA →	3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,653,266	0	A	G	100.0%	167.9 / NA	48	N258S (AAC→AGC)	sseA	3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/27);  total (21/27)

TAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGC  >  NC_000913/2653131‑2653401                                                                                                                                        |                                                                                                                                          tAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTg                                                                                                                                       <  1:782850/140‑1 (MQ=255)  aaaaaCGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGa                                                                                                                                      <  1:826802/140‑1 (MQ=255)  aaaaaCGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGa                                                                                                                                      <  1:193060/140‑1 (MQ=255)             aTGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTacgac                                                                                                                           <  1:1199397/140‑1 (MQ=255)             aTGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCCGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTacgac                                                                                                                           >  2:508024/1‑140 (MQ=255)                  cTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGgcgc                                                                                                                      >  2:1319762/1‑140 (MQ=255)                   tGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCa                                                                                                                     <  2:421247/140‑1 (MQ=255)                     gATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATg                                                                                                                   >  2:285442/1‑140 (MQ=255)                      aTGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATgg                                                                                                                  <  1:457814/140‑1 (MQ=255)                       tGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGa                                                                                                                 >  2:566876/1‑140 (MQ=255)                       tGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGa                                                                                                                 >  2:1006561/1‑140 (MQ=255)                        gCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAg                                                                                                                >  2:852189/1‑140 (MQ=255)                            tatTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTacgac                                                                                                                           <  1:9502/125‑1 (MQ=255)                            tatTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTacgac                                                                                                                           >  2:9502/1‑125 (MQ=255)                            tatTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGaa                                                                                                            <  1:85338/140‑1 (MQ=255)                                         cggcgTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGgcgcg                                                                                                   <  2:628672/136‑1 (MQ=255)                                         cggcgTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGgcgcg                                                                                                   >  1:628672/1‑136 (MQ=255)                                          ggcgTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGc                                                                                              <  1:688138/140‑1 (MQ=255)                                           gcgTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCa                                                                                             <  1:285442/140‑1 (MQ=255)                                                 gCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTa                                                                                       <  1:1140040/140‑1 (MQ=255)                                                          aaCCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAg                                                                              <  1:508024/140‑1 (MQ=255)                                                           aCCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTa                                                                                       >  2:985669/1‑130 (MQ=255)                                                           aCCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTa                                                                                       <  1:985669/130‑1 (MQ=255)                                                           aCCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGc                                                                             <  2:823068/140‑1 (MQ=255)                                                             cAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCa                                                                           <  1:852189/140‑1 (MQ=255)                                                               aTTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGt                                                                         >  2:313405/1‑140 (MQ=255)                                                               aTTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGt                                                                         >  2:694777/1‑140 (MQ=255)                                                                ttATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCgg                                                                                   >  2:1158826/1‑129 (MQ=255)                                                                ttATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCgg                                                                                   <  1:1158826/129‑1 (MQ=255)                                                                  aTCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGaa                                                                      <  1:694777/140‑1 (MQ=255)                                                                  aTCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGaa                                                                      >  2:1073594/1‑140 (MQ=255)                                                                    cGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCCGGATGTACCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAAt                                                                    <  1:289315/140‑1 (MQ=255)                                                                        aGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGt                                                                <  2:79197/140‑1 (MQ=255)                                                                          cTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTAt                                                              >  2:749705/1‑140 (MQ=255)                                                                             cGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTAtttt                                                           >  1:234660/1‑140 (MQ=255)                                                                              ggCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTa                                                          <  1:1006561/140‑1 (MQ=255)                                                                                ctctGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACa                                                        <  2:234660/140‑1 (MQ=255)                                                                                  ctGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAgg                                                      <  2:231066/140‑1 (MQ=255)                                                                                               aGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAAccg                                       >  1:283023/1‑140 (MQ=255)                                                                                               aGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAAccg                                       >  1:93398/1‑140 (MQ=255)                                                                                               aGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAAccg                                       >  2:1134178/1‑140 (MQ=255)                                                                                                         ttGTTAGCACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAACCGCCGAATTTgg                             <  2:93398/140‑1 (MQ=255)                                                                                                                cACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAAAACCGCCAAATTcggcgg                          <  1:1376485/138‑1 (MQ=255)                                                                                                                cACTCGCGACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAAAACCGCCAAATTcggcgg                          >  2:1376485/1‑138 (MQ=255)                                                                                                                    cgcgACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAACCGCCGAATTTGGCGGTTTTTTAt                  <  1:566876/140‑1 (MQ=255)                                                                                                                      cgACGCTGGATGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAACCGCCGAATTTGGCGGTTTTTTATTg                <  1:313405/140‑1 (MQ=255)                                                                                                                               aTGTGCCAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTgg       >  1:288892/1‑140 (MQ=255)                                                                                                                                    ccAAGCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgc  >  1:1320035/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                          TAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAAT‑‑AAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGC  >  NC_000913/2653131‑2653401

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