Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,283 A→T 31.0% A23A (GCT→GCA) ‡ cspC ← stress protein, member of the CspA‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,211,2830AT31.0% 251.4 / 134.4 152A23A (GCT→GCA) ‡cspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base A (47/57);  new base T (25/22);  total (73/79)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 155 total aligned reads displayed.

TCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATT  >  minE/1211215‑1211352
                                                                    |                                                                     
tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                                                                     >  1:3933927/1‑71 (MQ=255)
 cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAg                                                                      <  1:1528557/69‑1 (MQ=255)
 cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                                                                     <  1:1600824/70‑1 (MQ=255)
  aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg                                                                   <  1:2632015/71‑1 (MQ=255)
      agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGact                                                               <  1:2794389/71‑3 (MQ=255)
      agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAg                                                                >  1:5312850/1‑70 (MQ=255)
      agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                                                               <  1:1759384/71‑1 (MQ=255)
        agTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTaa                                                             <  1:4061192/71‑1 (MQ=16)
         gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                                                            <  1:5075619/71‑1 (MQ=255)
         gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                                                            >  1:2661813/1‑71 (MQ=255)
         gTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAAt                                                            <  1:1547175/71‑1 (MQ=16)
            ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa                                                          <  1:2655222/70‑1 (MQ=255)
            ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa                                                          <  1:2205268/70‑1 (MQ=255)
               aaaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAcc                                                      <  1:663283/68‑1 (MQ=12)
                aGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGcc                                                      >  1:431655/1‑70 (MQ=255)
                  ccATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa                                                   <  1:2245700/71‑1 (MQ=255)
                  ccATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa                                                   <  1:2066571/71‑1 (MQ=255)
                   cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGa                                                                 >  1:2481815/1‑56 (MQ=2)
                   cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCaaa                                                   <  1:951577/70‑1 (MQ=9)
                   cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCaaa                                                   <  1:5092732/70‑1 (MQ=9)
                     ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa                                                    <  1:3383326/67‑1 (MQ=255)
                        gtcTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtt                                              <  1:1279319/68‑1 (MQ=12)
                        gtcTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtt                                              <  1:4722556/68‑1 (MQ=12)
                        ctcTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt                                             <  1:2084618/69‑1 (MQ=9)
                        cccTGGATAGCGGAGAAGTGTACGAACACATCTTTGATGCCATCAGCCGGAGTAATGAAGCCAAAACCtt                                              >  1:4280542/1‑70 (MQ=255)
                        cccTGGATAGCGGAGAAGTGTAAGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt                                              <  1:4050024/70‑1 (MQ=255)
                         tcTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATcc                                                <  1:35234/66‑1 (MQ=12)
                          cTGGATAGCGGAGAAGTGTACGAACACATCTTTTCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt                                              >  1:2917133/1‑68 (MQ=255)
                          cTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAg                                           <  1:3836516/71‑1 (MQ=255)
                           tGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGaaac                                                       >  1:5419329/1‑56 (MQ=18)
                           tGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGa                                          >  1:290164/1‑71 (MQ=255)
                            ggATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGAc                                         >  1:1426323/1‑71 (MQ=255)
                             gATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGa                                          >  1:4172713/1‑69 (MQ=255)
                             gATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACt                                        <  1:523556/71‑1 (MQ=12)
                             gATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACt                                        >  1:2843540/1‑71 (MQ=255)
                                 gCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGc                                                       >  1:3857525/1‑52 (MQ=255)
                                 gCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                            <  1:3655965/63‑1 (MQ=255)
                                 gCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGtt                                    <  1:153090/71‑1 (MQ=255)
                                  cagagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtta                                   <  1:1167976/69‑2 (MQ=255)
                                  cagagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtt                                    <  1:3159219/68‑1 (MQ=12)
                                  cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTa                                                              <  1:1732270/44‑1 (MQ=255)
                                   ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                  <  1:1613290/71‑1 (MQ=255)
                                   ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                  <  1:303111/71‑1 (MQ=255)
                                   ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                  <  1:5539536/71‑1 (MQ=255)
                                    gagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGAc                                         >  1:3861295/1‑63 (MQ=12)
                                    gagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGAc                                         >  1:2598203/1‑63 (MQ=12)
                                    gagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCAGGAGTAATGAAACCGAATCCtttg                                            <  1:1869002/60‑2 (MQ=12)
                                    gagaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                                                            >  1:3001490/1‑44 (MQ=255)
                                       aaGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt                                              <  1:3736081/55‑1 (MQ=255)
                                        aGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCAc                             <  1:4487975/71‑1 (MQ=255)
                                          tgtACGAACACATCTTTGCTGCCATCAGCCGGAGTa                                                              <  1:2598985/36‑1 (MQ=255)
                                          tgtAAGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAccg                                                     <  1:1828006/45‑2 (MQ=2)
                                           gtACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATg                                                           >  1:3088576/1‑38 (MQ=18)
                                           gtACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg                                                           <  1:1027758/38‑1 (MQ=255)
                                            tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGaa                                                         <  1:5693580/39‑1 (MQ=17)
                                            tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTaa                         >  1:702603/1‑71 (MQ=255)
                                            tACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                                                            >  1:4421765/1‑36 (MQ=255)
                                            tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                            >  1:4843779/1‑52 (MQ=255)
                                            tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                            >  1:3796598/1‑52 (MQ=255)
                                            tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGt                                     >  1:2812273/1‑59 (MQ=255)
                                            tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                          <  1:2181620/70‑1 (MQ=255)
                                            tACGAACACATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCAtttta                         >  1:4551257/1‑69 (MQ=255)
                                             aCGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                                 >  1:1208086/1‑46 (MQ=255)
                                             aCGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                                 >  1:2164978/1‑46 (MQ=255)
                                              cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATg                                                           >  1:4512650/1‑35 (MQ=12)
                                              cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGa                                                          >  1:2849694/1‑36 (MQ=12)
                                              cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGa                                                          >  1:3804661/1‑36 (MQ=12)
                                              cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTaaa                                 >  1:3629146/1‑61 (MQ=9)
                                              cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAcg                       >  1:192844/1‑70 (MQ=12)
                                              cGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAcc                       <  1:4801092/71‑1 (MQ=255)
                                               gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                      >  1:3011678/1‑71 (MQ=255)
                                               gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                      >  1:4944556/1‑71 (MQ=255)
                                                aaCACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAcgtt                     <  1:4154567/71‑4 (MQ=12)
                                                aaCACATCTTTACTACCATCAACAGGAGAAATAAAGccga                                                    >  1:5206818/1‑38 (MQ=9)
                                                aaCACATCTTTACTACCATCAACAGGAGAAATAAAGccga                                                    >  1:1397473/1‑38 (MQ=9)
                                                 acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                    <  1:3624785/71‑1 (MQ=255)
                                                 acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                    >  1:2078914/1‑71 (MQ=255)
                                                 acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                    <  1:3971008/71‑1 (MQ=255)
                                                  cacaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCAc                             <  1:3855277/61‑1 (MQ=255)
                                                    caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCt                 <  1:4781149/71‑1 (MQ=255)
                                                      tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg                                   <  1:5642545/51‑1 (MQ=255)
                                                      tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                          <  1:4881688/60‑1 (MQ=255)
                                                      tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                          <  1:5158342/60‑1 (MQ=255)
                                                       cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                                 >  1:4205672/1‑36 (MQ=255)
                                                       cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa                         >  1:5689306/1‑60 (MQ=255)
                                                        tttGGTGCCATCAGCCGGAGTAATGAAGCCAAAACCt                                               >  1:5021466/1‑37 (MQ=255)
                                                        tttGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtttgg                                           <  1:1025835/41‑3 (MQ=12)
                                                          tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTaa             <  1:3695205/69‑1 (MQ=255)
                                                          tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa              <  1:3454643/68‑1 (MQ=255)
                                                           gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt          >  1:5137044/1‑71 (MQ=255)
                                                            ctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCtt         <  1:5041135/71‑1 (MQ=255)
                                                               ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGc      >  1:1797797/1‑71 (MQ=255)
                                                                 aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa    <  1:5384946/71‑1 (MQ=255)
                                                                  tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt          <  1:2535722/64‑1 (MQ=255)
                                                                  tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt          <  1:3575/64‑1 (MQ=255)
                                                                  tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt   <  1:1346808/71‑1 (MQ=255)
                                                                  tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt   >  1:5457350/1‑71 (MQ=255)
                                                                    aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa                         >  1:5227390/1‑47 (MQ=255)
                                                                    aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCttt        >  1:4610099/1‑64 (MQ=255)
                                                                    aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt  >  1:3981886/1‑70 (MQ=255)
                                                                    |                                                                     
TCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATT  >  minE/1211215‑1211352

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: