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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE899,6720TC21.9% 195.3 / 55.3 96L15P (CTG→CCG) pspEthiosulfate:cyanide sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (41/34);  new base C (10/11);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTC  >  minE/899602‑899740
                                                                      |                                                                    
aaTTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACt                                                                      >  1:2853038/1‑71 (MQ=255)
aaTTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACt                                                                      >  1:1674951/1‑71 (MQ=255)
 aTTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTg                                                                     >  1:4110586/1‑71 (MQ=255)
  ttttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGc                                                                    >  1:2253094/1‑71 (MQ=255)
    ttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGccc                                                                  >  1:1187676/1‑71 (MQ=255)
     tCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCg                                                                 >  1:1339472/1‑71 (MQ=255)
     tCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCg                                                                 >  1:2802077/1‑71 (MQ=255)
     tCATACATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc                                                                  >  1:1258939/1‑70 (MQ=255)
       aTCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGtt                                                               <  1:1831720/71‑1 (MQ=255)
          cATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGtttt                                                             <  1:282394/70‑1 (MQ=255)
              gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTTCCCGTTTTcgcc                                                         <  1:999460/70‑1 (MQ=255)
              gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTTCCCGTTTTcgcc                                                         <  1:3750951/70‑1 (MQ=255)
              gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgcc                                                         <  1:3867457/70‑1 (MQ=255)
              gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgcc                                                         <  1:3643922/70‑1 (MQ=255)
                 ggACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTg                                                     <  1:1385180/71‑1 (MQ=255)
                 ggACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTg                                                     <  1:30263/71‑1 (MQ=255)
                  gACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGa                                                    <  1:2393291/71‑1 (MQ=255)
                    cGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAac                                                  <  1:656659/71‑1 (MQ=255)
                              tttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCg                                         <  1:2500917/70‑1 (MQ=255)
                               ttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGaa                                                   <  1:3667855/59‑1 (MQ=255)
                               ttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGAt                                       >  1:3879977/1‑71 (MQ=255)
                                tAAAAAGGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                      >  1:542243/1‑71 (MQ=255)
                                tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                      >  1:3945242/1‑71 (MQ=255)
                                tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                      >  1:3184012/1‑71 (MQ=255)
                                tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                      >  1:1705725/1‑71 (MQ=255)
                                 aaaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgc                                                          >  1:3875342/1‑50 (MQ=255)
                                   aaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGc                                                                    <  1:4214504/38‑1 (MQ=255)
                                   aaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGc                                                                    <  1:2452887/38‑1 (MQ=255)
                                   aaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTc                                   >  1:3929391/1‑71 (MQ=255)
                                   aaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTc                                   >  1:3121831/1‑71 (MQ=255)
                                     aaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTc                                                            >  1:31702/1‑44 (MQ=255)
                                     aaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgt                                 <  1:3589835/71‑1 (MQ=255)
                                         cTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGtt                              >  1:2112625/1‑70 (MQ=255)
                                            aCTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAaca                                                 >  1:1905512/1‑48 (MQ=255)
                                             cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGaa                                                   >  1:1446593/1‑45 (MQ=255)
                                             cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGaa                                                   >  1:1242015/1‑45 (MQ=255)
                                             cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGaa                                                   >  1:1488700/1‑45 (MQ=255)
                                             cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCag                          <  1:3740876/70‑1 (MQ=255)
                                             cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATCTTCGTGTTCCag                          <  1:58908/70‑1 (MQ=255)
                                             cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTCTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCag                          <  1:865147/70‑1 (MQ=255)
                                              ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc                                          <  1:1694324/53‑1 (MQ=255)
                                              ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc                                          <  1:2951732/53‑1 (MQ=255)
                                              ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        >  1:2692740/1‑71 (MQ=255)
                                              ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        >  1:2547711/1‑71 (MQ=255)
                                              ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        >  1:3402059/1‑71 (MQ=255)
                                              ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        >  1:2666025/1‑71 (MQ=255)
                                              ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        >  1:964083/1‑71 (MQ=255)
                                              ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCaga                         >  1:869599/1‑70 (MQ=255)
                                              ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCag                          >  1:4222468/1‑69 (MQ=255)
                                               tGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        >  1:2592626/1‑70 (MQ=255)
                                                gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        <  1:1770664/69‑1 (MQ=255)
                                                gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        >  1:1460377/1‑69 (MQ=255)
                                                gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        >  1:4155046/1‑69 (MQ=255)
                                                gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        <  1:2695875/69‑1 (MQ=255)
                                                gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag                        >  1:589137/1‑69 (MQ=255)
                                                gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGc                       <  1:2110549/70‑1 (MQ=255)
                                                gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGc                       <  1:585265/70‑1 (MQ=255)
                                                gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGc                       <  1:3085268/70‑1 (MQ=255)
                                                  tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc                            >  1:1588973/1‑63 (MQ=255)
                                                  tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc                            >  1:2221842/1‑63 (MQ=255)
                                                  tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc                            >  1:3868698/1‑63 (MQ=255)
                                                  tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc                            >  1:1780545/1‑63 (MQ=255)
                                                  tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc                            >  1:1354717/1‑63 (MQ=255)
                                                      gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                  <  1:3951591/69‑1 (MQ=255)
                                                      gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                  >  1:467954/1‑69 (MQ=255)
                                                      gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                  >  1:2904664/1‑69 (MQ=255)
                                                      gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                  <  1:577329/69‑1 (MQ=255)
                                                      gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt                  <  1:2805925/69‑1 (MQ=255)
                                                        aCTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCg                                  <  1:3267146/51‑1 (MQ=255)
                                                         cTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCa             <  1:1891399/71‑1 (MQ=255)
                                                         cTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCa             <  1:3899413/71‑1 (MQ=255)
                                                          tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                      >  1:1829885/1‑45 (MQ=255)
                                                          tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                      >  1:3142369/1‑45 (MQ=255)
                                                          tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                      >  1:1959078/1‑45 (MQ=255)
                                                          tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg                                      >  1:2682987/1‑45 (MQ=255)
                                                          tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgtg                                >  1:21085/1‑51 (MQ=255)
                                                            gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAag           <  1:3288631/70‑1 (MQ=255)
                                                               ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgt                                 <  1:1706569/45‑1 (MQ=255)
                                                               ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgt                                 <  1:3290527/45‑1 (MQ=255)
                                                               ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgt                                 <  1:407799/45‑1 (MQ=255)
                                                               ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCa       <  1:3032676/71‑1 (MQ=255)
                                                                 ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc      >  1:2937631/1‑70 (MQ=255)
                                                                 ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg     >  1:146253/1‑71 (MQ=255)
                                                                  tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc      >  1:206264/1‑69 (MQ=255)
                                                                  tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc      >  1:593895/1‑69 (MQ=255)
                                                                  tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc      >  1:3675365/1‑69 (MQ=255)
                                                                  tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc      >  1:3200008/1‑69 (MQ=255)
                                                                  tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc      >  1:4022009/1‑69 (MQ=255)
                                                                  tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc      >  1:3354023/1‑69 (MQ=255)
                                                                  tCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg     <  1:3169548/70‑1 (MQ=255)
                                                                  tCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg     <  1:395561/70‑1 (MQ=255)
                                                                  tCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg     <  1:3279131/70‑1 (MQ=255)
                                                                  tCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg     <  1:3353920/70‑1 (MQ=255)
                                                                    aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc  <  1:1450144/71‑1 (MQ=255)
                                                                    aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc  <  1:710735/71‑1 (MQ=255)
                                                                    aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc  <  1:3495540/71‑1 (MQ=255)
                                                                    aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc  <  1:3706235/71‑1 (MQ=255)
                                                                    aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc  <  1:116511/71‑1 (MQ=255)
                                                                      |                                                                    
AATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTC  >  minE/899602‑899740

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: