Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 899,672 T→C 8.3% L15P (CTG→CCG)  pspE → thiosulfate:cyanide sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE899,6720TC 8.3% 194.9 / 9.3 72L15P (CTG→CCG) pspEthiosulfate:cyanide sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (27/39);  new base C (4/2);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTT  >  minE/899606‑899739
                                                                  |                                                                   
ttcattcaTAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc                                                                 >  1:1984197/1‑71 (MQ=255)
ttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc                                                                 >  1:1418498/1‑71 (MQ=255)
ttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc                                                                 >  1:3079218/1‑71 (MQ=255)
ttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc                                                                 >  1:3275012/1‑71 (MQ=255)
  cATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGt                                                               >  1:1828104/1‑71 (MQ=255)
      cATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTc                                                           <  1:3005134/71‑1 (MQ=255)
      cATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTc                                                           <  1:520483/71‑1 (MQ=255)
        tAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTc                                                           <  1:3014207/69‑1 (MQ=255)
          gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgccg                                                       <  1:49542/71‑1 (MQ=255)
          gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgccg                                                       <  1:398955/71‑1 (MQ=255)
           aaGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgccg                                                       <  1:3491836/70‑1 (MQ=255)
               aCGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGtt                                                              <  1:661903/59‑1 (MQ=255)
                cGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGaa                                                  >  1:3876223/1‑70 (MQ=255)
                cGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGaa                                                  >  1:3659352/1‑70 (MQ=255)
                 gCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGtttt                                                            >  1:3254336/1‑59 (MQ=255)
                    tACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTg                                             >  1:2800224/1‑71 (MQ=255)
                    tACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTg                                             >  1:3760399/1‑71 (MQ=255)
                     aCATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTg                                             >  1:3264434/1‑70 (MQ=255)
                     aCATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTg                                             >  1:61421/1‑70 (MQ=255)
                          tttAAAAAATGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa                                       <  1:3738658/71‑1 (MQ=255)
                          tttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTgg                                            <  1:2966973/66‑1 (MQ=255)
                          tttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa                                       >  1:2384868/1‑71 (MQ=255)
                          tttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa                                       <  1:2547826/71‑1 (MQ=255)
                            tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGAt                                          >  1:2123320/1‑66 (MQ=255)
                            tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGAt                                          >  1:282300/1‑66 (MQ=255)
                            tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                     >  1:1930860/1‑71 (MQ=255)
                            tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                     >  1:209792/1‑71 (MQ=255)
                            tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                     >  1:3051776/1‑71 (MQ=255)
                              aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGt                                    <  1:1036882/70‑1 (MQ=255)
                              aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGt                                    <  1:1003257/70‑1 (MQ=255)
                              aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGt                                    <  1:168166/70‑1 (MQ=255)
                               aaaaGGCTTACTTGCTCTGTCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTTGATCGATGTTc                                  <  1:1717593/71‑1 (MQ=255)
                                      ttACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTc                            >  1:1688884/1‑70 (MQ=255)
                                           tGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                       >  1:3362055/1‑70 (MQ=255)
                                            gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTc                            >  1:3097274/1‑64 (MQ=255)
                                            gCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                       >  1:3352755/1‑69 (MQ=255)
                                                tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc                                         <  1:3894283/47‑1 (MQ=255)
                                                tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc                                         <  1:1799197/47‑1 (MQ=255)
                                                tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc                                         <  1:1803889/47‑1 (MQ=255)
                                                tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc                                         <  1:1935866/47‑1 (MQ=255)
                                                tGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc                                         <  1:1280034/47‑1 (MQ=255)
                                                tGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGt                   <  1:3243145/69‑1 (MQ=255)
                                                  gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                 <  1:486017/69‑1 (MQ=255)
                                                  gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                 <  1:1312394/69‑1 (MQ=255)
                                                  gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                 <  1:2936017/69‑1 (MQ=255)
                                                  gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                 <  1:2918211/69‑1 (MQ=255)
                                                  gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                 >  1:2877704/1‑69 (MQ=255)
                                                  gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                 <  1:1780944/69‑1 (MQ=255)
                                                  gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                 >  1:3680247/1‑69 (MQ=255)
                                                  gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCa               <  1:1379736/71‑1 (MQ=255)
                                                     cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa                                       <  1:918905/44‑1 (MQ=255)
                                                     cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa                                       <  1:1539828/44‑1 (MQ=255)
                                                     cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa                                       <  1:3299230/44‑1 (MQ=255)
                                                     cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa                                       <  1:2597291/44‑1 (MQ=255)
                                                     cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGa                                       <  1:2388946/44‑1 (MQ=255)
                                                        gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCaa           >  1:496611/1‑69 (MQ=255)
                                                         tgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAg                    <  1:3828434/59‑1 (MQ=255)
                                                           ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCATAGCAGTATCAGc             <  1:2605134/64‑1 (MQ=255)
                                                           ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAg                    <  1:3365759/57‑1 (MQ=255)
                                                            tttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGc             <  1:2847840/63‑1 (MQ=255)
                                                             ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCTATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg    >  1:757099/1‑71 (MQ=255)
                                                             ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgtg                               >  1:3905905/1‑44 (MQ=255)
                                                             ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgtg                               >  1:3912010/1‑44 (MQ=255)
                                                             ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgtg                               >  1:2279502/1‑44 (MQ=255)
                                                             ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg    >  1:440072/1‑71 (MQ=255)
                                                              tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg    <  1:3221734/70‑1 (MQ=255)
                                                              tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg    <  1:798017/70‑1 (MQ=255)
                                                               cACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGtt  >  1:420133/1‑71 (MQ=255)
                                                                aCTGCCCGTTTTCGCCGCTGAACACTGTATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGtt  <  1:1848761/70‑1 (MQ=255)
                                                                aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGtt  <  1:3885477/70‑1 (MQ=255)
                                                                aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGtt  <  1:2685603/70‑1 (MQ=255)
                                                                aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGAGTTCCAGAGCAGTATCAGCAAGAGCACGtt  <  1:2317800/70‑1 (MQ=255)
                                                                  |                                                                   
TTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTT  >  minE/899606‑899739

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