Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,652,948	A→C	E152A (GAA→GCA)	sseA →	3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,652,948	0	A	C	100.0%	159.4 / NA	48	E152A (GAA→GCA)	sseA	3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/27);  total (21/27)

GAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACT  >  NC_000913/2652801‑2653092                                                                                                                                                    |                                                                                                                                                 gAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGc                                                                                                                                               >  1:133704/1‑151 (MQ=255)     aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCgt                                                                                                                                            <  1:601116/150‑1 (MQ=255)          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGaa                                                                                                                                      <  2:503377/151‑1 (MQ=255)               ctgcGCACCTTTGGTGTCGAGAAAGTGCCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTa                                                                                                                                  >  1:255248/1‑150 (MQ=255)                 gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTaa                                                                                                                                 <  1:838705/149‑1 (MQ=255)                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGAt                                                                                                                            >  1:817212/1‑151 (MQ=255)                     aCCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATg                                                                                                                           >  2:469570/1‑151 (MQ=255)                            gtgtAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTattatt                                                                                                                    <  1:16624/151‑1 (MQ=255)                                           cGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGaaa                                                                                                      >  2:677823/1‑150 (MQ=255)                                               tCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATAc                                                                                                  <  1:673795/150‑1 (MQ=255)                                               tCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACg                                                                                                 >  2:629391/1‑151 (MQ=255)                                                          gACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAattat                                                                                      <  2:401782/151‑1 (MQ=255)                                                              tGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTg                                                                                    <  1:352667/149‑1 (MQ=255)                                                              tGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGAt                                                                                  >  2:665157/1‑151 (MQ=255)                                                               gCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGAt                                                                                  <  2:923924/150‑1 (MQ=255)                                                                   gCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATgcccg                                                                             <  1:497922/151‑1 (MQ=255)                                                                   gCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATgcccg                                                                             <  1:858653/151‑1 (MQ=255)                                                                     tGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATgcccgcc                                                                           >  2:546622/1‑151 (MQ=255)                                                                      ggCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATgcccgcc                                                                           <  2:44097/150‑1 (MQ=255)                                                                        cAGCGCGATGATCTGTTGTTTGGTGATGGTGCATTAGTGCTGCCGGTGGCGGGATGTATCGCCCCGTTTAATCTTGCAGGCGTATTGATGGTAACCGAGGTATTATTGGCAAGCCATGAAAATACGGCACAAATTATTGATgcccgc                                                                            <  1:44397/147‑1 (MQ=255)                                                                           cgcgATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTg                                                                     <  2:600433/151‑1 (MQ=255)                                                                             cgATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGc                                                                    >  2:825198/1‑150 (MQ=255)                                                                              gatgatCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCAc                                                                  >  2:673172/1‑151 (MQ=255)                                                                                      gctgTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGtt                                                               <  1:280330/146‑1 (MQ=255)                                                                                      gctgTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTa                                                            <  1:815635/149‑1 (MQ=255)                                                                                         gTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTACCGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACg                                                         <  1:249914/149‑1 (MQ=255)                                                                                         gTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGc                                                        <  2:429890/150‑1 (MQ=255)                                                                                            agaagaAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGaa                                                    <  2:655431/151‑1 (MQ=255)                                                                                             gaagaaGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGa                                                     >  2:507968/1‑149 (MQ=255)                                                                                              aagaagGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGt                                                  <  1:478852/151‑1 (MQ=255)                                                                                                gaagGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTGTTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGtt                                                 >  1:758055/1‑150 (MQ=255)                                                                                                gaagGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTtg                                                >  2:277995/1‑151 (MQ=255)                                                                                                    gTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTtgatg                                             >  2:880254/1‑150 (MQ=255)                                                                                                      gCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAAc                                          >  2:17032/1‑151 (MQ=255)                                                                                                       cAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAAcc                                         <  1:893839/151‑1 (MQ=255)                                                                                                          tAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTc                                       >  2:677247/1‑150 (MQ=255)                                                                                                            gagCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGc                                     <  1:893228/150‑1 (MQ=255)                                                                                                                tGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAgg                                >  2:589792/1‑151 (MQ=255)                                                                                                                 gCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAgg                                <  2:379643/150‑1 (MQ=255)                                                                                                                   cGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGttt                             >  2:221814/1‑151 (MQ=255)                                                                                                                    ggAAGGAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGttt                             <  1:524416/150‑1 (MQ=255)                                                                                                                        ggAGAGTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTc                        <  1:726207/151‑1 (MQ=255)                                                                                                                            agTTTAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCgg                    >  1:889483/1‑151 (MQ=255)                                                                                                                            agTTTAACGCCGCGCTTCATCCTGCAGGTGTGGTGAAAGTAACCGCTGTATTATTGGCAAGCCATGAAAATACGGCGCATATTCTTGAAGCACGCACGGCTGCACGTTTTAACGCAGAACTTTTTGAACGTCGCCCAGGTTTCCGCCGCgg                    >  1:266430/1‑151 (MQ=255)                                                                                                                                tAACGCCGCGTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACa                 <  2:454746/150‑1 (MQ=255)                                                                                                                                      cgcgTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTccc          <  1:485122/151‑1 (MQ=255)                                                                                                                                         gTTTAATCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCg         >  1:576725/1‑149 (MQ=255)                                                                                                                                              aTCCTGCAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACt  <  1:38068/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACT  >  NC_000913/2652801‑2653092

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