Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,686,595 G→A 39.0% intergenic (+570/+248) sseA → / ← sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,686,5950GA39.0% 60.1 / 72.4 110intergenic (+570/+248)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base A (32/11);  total (64/46)
Fisher's exact test for biased strand distribution p-value = 6.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

TATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGTG  >  minE/1686530‑1686661
                                                                 |                                                                  
tattatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTgg                                                                 >  1:3143464/1‑69 (MQ=255)
tattatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTgg                                                                 >  1:1519798/1‑69 (MQ=255)
tattatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTgg                                                                 >  1:2103023/1‑69 (MQ=255)
tattatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTgg                                                                 >  1:2573279/1‑69 (MQ=255)
tattatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTgg                                                                 >  1:2775626/1‑69 (MQ=255)
tattatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTgg                                                                 >  1:772799/1‑69 (MQ=255)
tattatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTgg                                                                 >  1:831963/1‑69 (MQ=255)
tattatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTgg                                                                 >  1:1111751/1‑69 (MQ=255)
  ttatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGc                                                                >  1:455555/1‑68 (MQ=255)
  ttatCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTACGTGGc                                                                >  1:1041371/1‑68 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGc                                                       >  1:669705/1‑67 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCa                                                      >  1:443962/1‑68 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCa                                                      >  1:421254/1‑68 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCa                                                      >  1:2737942/1‑68 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCa                                                      >  1:3077918/1‑68 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCa                                                      >  1:2471740/1‑68 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCa                                                      >  1:2336251/1‑68 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCa                                                      >  1:924199/1‑68 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCa                                                      >  1:2015104/1‑68 (MQ=255)
            cGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCa                                                      >  1:1636483/1‑68 (MQ=255)
              gCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGgt                                                   <  1:3235327/69‑1 (MQ=255)
              gCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGgt                                                   <  1:3228337/69‑1 (MQ=255)
              gCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGgt                                                   <  1:858722/69‑1 (MQ=255)
              gCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGgt                                                   <  1:3432174/69‑1 (MQ=255)
              gCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGgt                                                   <  1:2793721/69‑1 (MQ=255)
              gCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGgt                                                   <  1:2785577/69‑1 (MQ=255)
              gCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGgt                                                   <  1:513838/69‑1 (MQ=255)
              gCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGgt                                                   <  1:1200982/69‑1 (MQ=255)
                                     gggCGCTGGCTTAAGGAAAGGATGTTCCATGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATcgcg                            <  1:2000935/69‑1 (MQ=255)
                                     gggCGCTGGCTTAAGGAAAGGATGTTCCATGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATcgcg                            <  1:1400434/69‑1 (MQ=255)
                                     gggCGCTGGCTTAAGGAAAGGATGTTCCATGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATcgcg                            <  1:149546/69‑1 (MQ=255)
                                     gggCGCTGGCTTAAGGAAAGGATGTTCCATGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATcgcg                            <  1:118556/69‑1 (MQ=255)
                                     gggCGCTGGCTTAAGGAAAGGATGTTCCATGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATcgcg                            <  1:756486/69‑1 (MQ=255)
                                     gggCGCTGGCTTAAGGAAAGGATGTTCCATGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATcgcg                            <  1:1141136/69‑1 (MQ=255)
                                     gggCGCTGGCTTAAGGAAAGGATGTTCCATGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATcgcg                            <  1:782447/69‑1 (MQ=255)
                                     gggCGCTGGCTTAAGGAAAGGATGTTCCATGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATcgcg                            <  1:2755511/69‑1 (MQ=255)
                                     gggCGCTGGCTTAAGGAAAGGATGTTCCATGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATcgcg                            <  1:2650765/69‑1 (MQ=255)
                                       gcgcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:3363755/69‑1 (MQ=255)
                                       gcgcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:3021439/69‑1 (MQ=255)
                                       gcgcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:241162/69‑1 (MQ=255)
                                       gcgcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:326888/69‑1 (MQ=255)
                                       gcgcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:3306383/69‑1 (MQ=255)
                                       gcgcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:3334910/69‑1 (MQ=255)
                                       gcgcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:639010/69‑1 (MQ=255)
                                       gcgcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:1093223/69‑1 (MQ=255)
                                       gcgcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTACCGCGGc                          <  1:437506/69‑1 (MQ=255)
                                         gcTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCaa                        <  1:1514847/69‑1 (MQ=255)
                                                    gAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAgt              <  1:1372925/68‑1 (MQ=255)
                                                    gAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGACAATATCGCCAgt              <  1:1317994/68‑1 (MQ=255)
                                                     aaaGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCGTGCTATCGCGGCAATATCGCCAgt              <  1:3190136/67‑1 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1740970/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:2007351/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1345691/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1979194/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:519062/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:2077256/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:558031/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1837287/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:2225926/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:2288466/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1828406/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1766945/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:2377379/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:2403534/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1127802/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:3151518/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1586384/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:3131543/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:854044/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1478564/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:282306/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1723445/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:886611/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1699707/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1410006/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:1475056/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:925021/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTcaca                                           >  1:948381/1‑38 (MQ=255)
                                                     aaaGGATGTTCCATGGCCGTAAATGCAGGTGTTTc                                              >  1:2721916/1‑35 (MQ=255)
                                                         gATGTTCCGTGGCAGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCt        <  1:2555595/69‑1 (MQ=255)
                                                          aTGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:100427/50‑1 (MQ=255)
                                                          aTGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:2794104/50‑1 (MQ=255)
                                                          aTGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:1852495/50‑1 (MQ=255)
                                                          aTGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:2463291/50‑1 (MQ=255)
                                                          aTGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:754980/50‑1 (MQ=255)
                                                          aTGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:330515/50‑1 (MQ=255)
                                                          aTGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:847347/50‑1 (MQ=255)
                                                          aTGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGc                          <  1:3158841/50‑1 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGt   >  1:2609570/1‑68 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGt   >  1:1513205/1‑68 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGt   >  1:1481704/1‑68 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGt   <  1:2093352/68‑1 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGt   >  1:402644/1‑68 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGTg  >  1:3152658/1‑69 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGTg  >  1:33495/1‑69 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGTg  >  1:580670/1‑69 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGTg  >  1:2545595/1‑69 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGTg  >  1:896616/1‑69 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGTg  >  1:357936/1‑69 (MQ=255)
                                                               ccGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGTg  >  1:928735/1‑69 (MQ=255)
                                                                 |                                                                  
TATTATCCACTGCGGCCCTTTCCGCCGTCTCGCAAACGGGCGCTGGCTTTAGGAAAGGATGTTCCGTGGCCGTAAATGCAGGTGTTTCACAGCGCTTGCTATCGCGGCAATATCGCCAGTGGTGCTGTCGTG  >  minE/1686530‑1686661

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