Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,369,780 A→G Q31R (CAG→CGG)  pspE → thiosulfate:cyanide sulfurtransferase (rhodanese)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,369,7800AG100.0% 222.6 / NA 62Q31R (CAG→CGG) pspEthiosulfate:cyanide sulfurtransferase (rhodanese)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/29);  total (33/29)

ATAATGGATGGCAATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAG  >  NC_000913/1369650‑1369914
                                                                                                                                  |                                                                                                                                      
ataATGGATGGCAATTTTCATCCATAGAAGGACGCTTACATGATTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGc                                                                                                                                >  1:1002746/1‑139 (MQ=255)
      gatggCAATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCa                                                                                                                               >  2:1235321/1‑134 (MQ=255)
      gatggCAATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCa                                                                                                                               <  1:1235321/134‑1 (MQ=255)
      gatggCAATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGc                                                                                                                          >  2:950320/1‑139 (MQ=255)
       atggCAATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCa                                                                                                                         >  1:1298517/1‑139 (MQ=255)
         ggCAATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACg                                                                                                                       >  2:1386854/1‑139 (MQ=255)
         ggCAATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACg                                                                                                                       >  1:789762/1‑139 (MQ=255)
                ttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAggg                                                                                                                <  2:1002746/139‑1 (MQ=255)
                   aTCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGc                                                                                                             <  2:1298517/139‑1 (MQ=255)
                        tAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCa                                                                                                        >  1:69280/1‑139 (MQ=255)
                         aGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCaa                                                                                                       >  2:1312158/1‑139 (MQ=255)
                           aaGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAata                                                                                                     <  2:69280/139‑1 (MQ=255)
                             ggACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATAtt                                                                                                   <  1:291293/139‑1 (MQ=255)
                                cGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTccc                                                                                                >  1:559205/1‑139 (MQ=255)
                                 gCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTcccc                                                                                               >  2:1243054/1‑139 (MQ=255)
                                  cTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCa                                                                                                                               <  1:347039/106‑1 (MQ=255)
                                  cTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCa                                                                                                                               >  2:347039/1‑106 (MQ=255)
                                  cTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCt                                                                                              >  2:168169/1‑139 (MQ=255)
                                  cTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCt                                                                                              >  2:809670/1‑139 (MQ=255)
                                  cTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCt                                                                                              >  1:1198060/1‑139 (MQ=255)
                                   ttACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTg                                                                                             >  2:1311674/1‑139 (MQ=255)
                                     aCATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaa                                                                                           <  1:1386854/139‑1 (MQ=255)
                                     aCATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaa                                                                                           <  1:1312158/139‑1 (MQ=255)
                                      cATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaag                                                                                          >  1:1434347/1‑139 (MQ=255)
                                           ttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAgggg                                                                                                               >  1:717408/1‑113 (MQ=255)
                                           ttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAgggg                                                                                                               <  2:717408/113‑1 (MQ=255)
                                                       tACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGAGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt                                                                         <  2:789762/139‑1 (MQ=255)
                                                              tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAcc                                                                  >  2:730255/1‑139 (MQ=255)
                                                               ctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccg                                                                 <  1:1311674/139‑1 (MQ=255)
                                                                       ggTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCg                                                         >  2:448133/1‑139 (MQ=255)
                                                                          gtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                                                      <  2:1198060/139‑1 (MQ=255)
                                                                                  tGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGa                                              >  2:1374984/1‑139 (MQ=255)
                                                                                  tGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGa                                              >  2:1190943/1‑139 (MQ=255)
                                                                                       gTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCg                                         <  2:997732/139‑1 (MQ=255)
                                                                                       gTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCg                                         <  2:1434347/139‑1 (MQ=255)
                                                                                        ttttCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccg                                                                 >  2:635594/1‑114 (MQ=255)
                                                                                        ttttCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccg                                                                 <  1:635594/114‑1 (MQ=255)
                                                                                                    aaCACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccg                                                                 <  1:1450096/102‑1 (MQ=255)
                                                                                                    aaCACTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccg                                                                 >  2:1450096/1‑102 (MQ=255)
                                                                                                     acacTGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGc                           <  1:1243054/139‑1 (MQ=255)
                                                                                                         tGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGa                                      >  2:807432/1‑124 (MQ=255)
                                                                                                         tGGATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGa                                      <  1:807432/124‑1 (MQ=255)
                                                                                                           gATCGATGTTCGTGTTCCAGAGCGGTATCAGCAagag                                                                                                                           <  1:1289594/37‑1 (MQ=255)
                                                                                                           gATCGATGTTCGTGTTCCAGAGCGGTATCAGCAagag                                                                                                                           >  2:1289594/1‑37 (MQ=255)
                                                                                                           gATCGATGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGcc                     >  1:1090566/1‑139 (MQ=255)
                                                                                                                aTGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtg                                 <  1:388428/122‑1 (MQ=255)
                                                                                                                aTGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtg                                 >  2:388428/1‑122 (MQ=255)
                                                                                                                aTGTTCGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACg                <  1:1374984/139‑1 (MQ=255)
                                                                                                                     cGTGTTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGACACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGt           >  2:930958/1‑139 (MQ=255)
                                                                                                                       tgtTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGcc                     <  1:1142640/127‑1 (MQ=255)
                                                                                                                       tgtTCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGcc                     >  2:1142640/1‑127 (MQ=255)
                                                                                                                         ttCCAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAgg       <  1:930958/139‑1 (MQ=255)
                                                                                                                           ccAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaagg                                                                                         >  1:1128986/1‑55 (MQ=255)
                                                                                                                           ccAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTgaagg                                                                                         <  2:1128986/55‑1 (MQ=255)
                                                                                                                           ccAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGt           <  1:405664/133‑1 (MQ=255)
                                                                                                                           ccAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGt           >  2:405664/1‑133 (MQ=255)
                                                                                                                           ccAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGgc     <  1:1190943/139‑1 (MQ=255)
                                                                                                                            cAGAGCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGgca    <  1:168169/139‑1 (MQ=255)
                                                                                                                             agagCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGgcaa   >  1:518031/1‑139 (MQ=255)
                                                                                                                              gagCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGgcaag  <  1:730255/139‑1 (MQ=255)
                                                                                                                               agCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTc                                                           >  1:1173750/1‑81 (MQ=255)
                                                                                                                               agCGGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTc                                                           <  2:1173750/81‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                      
ATAATGGATGGCAATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAG  >  NC_000913/1369650‑1369914

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