Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,933 C→T 20.8% S381S (TCC→TCT) ‡ cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,9330CT20.8% 203.0 / 44.2 101S381S (TCC→TCT) ‡cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (41/39);  new base T (14/7);  total (55/46)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

Only 100 of 104 total aligned reads displayed.

GGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCT  >  minE/385863‑385989
                                                                      |                                                        
ggCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTcc                                                          <  1:1768261/71‑1 (MQ=255)
ggCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTcc                                                          <  1:5588546/71‑1 (MQ=255)
  cAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGc                                                        >  1:4717433/1‑71 (MQ=255)
     gCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTg                                                     <  1:5210306/71‑1 (MQ=255)
      agTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGt                                                    <  1:171864/70‑1 (MQ=255)
        tGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTat                                                  >  1:5383460/1‑71 (MQ=255)
          cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatca                                                >  1:796552/1‑71 (MQ=255)
          cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatca                                                <  1:2988911/71‑1 (MQ=255)
          cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatc                                                 <  1:382975/70‑1 (MQ=255)
          cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatc                                                 <  1:2660878/70‑1 (MQ=255)
            tcaTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                                              >  1:1024988/3‑71 (MQ=255)
            tcaTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                                              >  1:2478888/3‑71 (MQ=255)
            tcaTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                                              >  1:4048213/3‑71 (MQ=255)
            tcaTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                                              >  1:5459176/3‑71 (MQ=255)
            tcaTCGGCTACGGAGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                                              >  1:1412700/3‑71 (MQ=255)
             tATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTcc                                                          >  1:5249769/1‑58 (MQ=255)
             caTCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCatcatc                                              >  1:2165164/2‑70 (MQ=255)
              aTCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCg                                             <  1:2290757/70‑1 (MQ=255)
               tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGt                                                    >  1:3036305/1‑62 (MQ=255)
               tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGt                                                    >  1:3774717/1‑62 (MQ=255)
               tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatca                                                <  1:43401/66‑1 (MQ=255)
               tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                                            <  1:564175/70‑1 (MQ=255)
               tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                                            <  1:5618758/70‑1 (MQ=255)
               tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                                            <  1:3811882/70‑1 (MQ=255)
               tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                                            <  1:342877/70‑1 (MQ=255)
               tCGGCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGaa                                           >  1:4965950/1‑70 (MQ=255)
                 ggCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatcatc                                              >  1:3924332/1‑66 (MQ=255)
                 ggCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatcatc                                              >  1:4375191/1‑66 (MQ=255)
                 ggCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatcatc                                              >  1:1404056/1‑66 (MQ=255)
                 ggCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGaac                                          <  1:401342/70‑3 (MQ=255)
                 ggCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAAcc                                         >  1:3009516/1‑71 (MQ=255)
                 ggCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAAcc                                         <  1:1386348/71‑1 (MQ=255)
                  gCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCg                                        <  1:396407/71‑1 (MQ=255)
                  gCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCg                                        <  1:1752776/71‑1 (MQ=255)
                      cGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCggg                                      <  1:1775997/69‑1 (MQ=255)
                         tGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTg                                                     <  1:987350/51‑1 (MQ=255)
                          gggATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTa                                >  1:4061144/1‑71 (MQ=255)
                            gATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatcatc                                              <  1:2527693/55‑1 (MQ=255)
                            gATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAtt                              <  1:3964663/71‑1 (MQ=255)
                             aTGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATGGATCCGGGCGTGTAttt                             >  1:319888/1‑71 (MQ=255)
                             aTGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                             >  1:1417581/1‑71 (MQ=255)
                               gTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTg                           >  1:1918695/1‑71 (MQ=255)
                                ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                                            <  1:229017/53‑1 (MQ=255)
                                ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                                            <  1:2767376/53‑1 (MQ=255)
                                ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGa                                            <  1:782660/53‑1 (MQ=255)
                                ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGAt                                           >  1:3080136/1‑54 (MQ=255)
                                ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGc                          >  1:3382033/1‑71 (MQ=255)
                                 tAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatc                                                 >  1:536425/1‑47 (MQ=255)
                                 tAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTatc                                                 >  1:1563734/1‑47 (MQ=255)
                                  aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGc                                     >  1:5313509/1‑58 (MQ=255)
                                  aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGc                                     >  1:5671578/1‑58 (MQ=255)
                                  aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                             >  1:957108/1‑66 (MQ=255)
                                  aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGc                          >  1:3199603/1‑69 (MQ=255)
                                  aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                        >  1:4313481/1‑71 (MQ=255)
                                  aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                        >  1:442668/1‑71 (MQ=255)
                                  aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                        >  1:2477081/1‑71 (MQ=255)
                                   aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCa                        >  1:2778469/1‑70 (MQ=255)
                                   aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAt                       <  1:3215696/71‑1 (MQ=255)
                                   aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAt                       <  1:4172047/71‑1 (MQ=255)
                                    tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTa                                                   <  1:2489596/42‑1 (MQ=255)
                                    tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTa                                                   <  1:5695117/42‑1 (MQ=255)
                                    tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTa                                                   <  1:690437/42‑1 (MQ=255)
                                    tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTa                                                   <  1:1225190/42‑1 (MQ=255)
                                    tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                             >  1:4449895/1‑64 (MQ=255)
                                    tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                      <  1:4695503/71‑1 (MQ=255)
                                    tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                      <  1:5561103/71‑1 (MQ=255)
                                    tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                      <  1:3588898/71‑1 (MQ=255)
                                     ggAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTCTATTTTGCCATg                      >  1:5022132/1‑70 (MQ=255)
                                      gAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTg                                                     >  1:5454546/1‑38 (MQ=255)
                                      gAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGa                     <  1:1944747/70‑1 (MQ=255)
                                        aTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTg                                                     <  1:696274/36‑1 (MQ=255)
                                         tCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACac                 <  1:1128206/71‑2 (MQ=255)
                                           cTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTat                                                  >  1:4934104/1‑36 (MQ=255)
                                           cTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTat                                                  >  1:140140/1‑36 (MQ=255)
                                           cTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTat                                                  >  1:5364521/1‑36 (MQ=255)
                                           cTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTa                                                   >  1:2975129/1‑35 (MQ=255)
                                           cTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGcc               >  1:144419/1‑71 (MQ=255)
                                           cTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc               >  1:5577617/1‑71 (MQ=255)
                                             tCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCg             >  1:1968933/1‑71 (MQ=255)
                                              cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGAt                                           >  1:2813461/1‑40 (MQ=255)
                                              cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGAt                                           >  1:3225489/1‑40 (MQ=255)
                                              cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGAt                                           >  1:1710372/1‑40 (MQ=255)
                                              cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGAt                                           >  1:3701285/1‑40 (MQ=255)
                                              cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGAt                                           >  1:5564122/1‑40 (MQ=255)
                                              cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc              >  1:4800594/1‑69 (MQ=255)
                                              cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc              <  1:5523451/69‑1 (MQ=255)
                                              cGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc            >  1:1954198/1‑70 (MQ=255)
                                              cGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc            >  1:2396515/1‑70 (MQ=255)
                                              cGTGGCGATTACGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc              <  1:1744653/69‑1 (MQ=255)
                                                  gCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc            <  1:1177082/67‑2 (MQ=255)
                                                  gCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc            <  1:1170488/67‑2 (MQ=255)
                                                     aTTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt     <  1:5356536/71‑1 (MQ=255)
                                                       tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc               >  1:2232338/1‑59 (MQ=255)
                                                       tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc               >  1:122627/1‑59 (MQ=255)
                                                         tGGCGCTGGTTTCCGCCTGTATCATCGAGCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCt  <  1:2370925/70‑1 (MQ=255)
                                                              cTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATg          >  1:4033427/1‑57 (MQ=255)
                                                                ggTTTCCGCCTGTATCATCGATCCTGGCGTGTATTTTGCCATGAACAGccc              <  1:46340/51‑1 (MQ=255)
                                                                ggTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc              <  1:1721911/51‑1 (MQ=255)
                                                                    tCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                      <  1:5483642/39‑1 (MQ=255)
                                                                    tCCGCCTGTATCATCGATCCGGGCGTGTATTTAGCCATg                      <  1:1252400/39‑1 (MQ=255)
                                                                      |                                                        
GGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCT  >  minE/385863‑385989

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: