Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,286 C→A 7.9% P22P (CCG→CCT 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,2860CA 7.9% 422.1 / 14.4 150P22P (CCG→CCTcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base C (66/72);  new base A (8/4);  total (74/76)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 161 total aligned reads displayed.

AGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTG  >  minE/1211217‑1211354
                                                                     |                                                                    
   aagagTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAg                                                                  >  1:2641007/2‑71 (MQ=255)
    agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGact                                                                 <  1:2794389/71‑3 (MQ=255)
    agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAg                                                                  >  1:5312850/1‑70 (MQ=255)
    agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                                                                 <  1:1759384/71‑1 (MQ=255)
      agTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa                                                               <  1:5371059/71‑1 (MQ=255)
       gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                                                              <  1:5394219/71‑1 (MQ=255)
       gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                                                              >  1:2186337/1‑71 (MQ=255)
       gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                                                              <  1:5075619/71‑1 (MQ=255)
       gTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAAt                                                              <  1:1547175/71‑1 (MQ=16)
          ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa                                                            <  1:2655222/70‑1 (MQ=255)
             aaaccATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAcc                                                        <  1:663283/68‑1 (MQ=12)
                ccATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa                                                     <  1:2245700/71‑1 (MQ=255)
                ccATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa                                                     <  1:2066571/71‑1 (MQ=255)
                 cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGa                                                                   >  1:3652710/1‑56 (MQ=2)
                 cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGa                                                                   >  1:2481815/1‑56 (MQ=2)
                 cATTGCTCTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCaaa                                                     <  1:951577/70‑1 (MQ=9)
                 cATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGa                                                                   >  1:2581149/1‑56 (MQ=12)
                   ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa                                                      <  1:3383326/67‑1 (MQ=255)
                      gtcTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtt                                                <  1:4722556/68‑1 (MQ=12)
                      gtcTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtt                                                <  1:1279319/68‑1 (MQ=12)
                      cccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt                                                <  1:3886370/70‑1 (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)
                         tGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtt                                                >  1:4297971/1‑67 (MQ=12)
                            aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCgg                                                                    >  1:5179639/1‑44 (MQ=255)
                               gCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGtt                                      <  1:153090/71‑1 (MQ=255)
                                cagagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtta                                     <  1:1167976/69‑2 (MQ=255)
                                cagagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtt                                      >  1:4202839/3‑70 (MQ=12)
                                cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                              <  1:5060615/62‑1 (MQ=255)
                                 ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                    <  1:1613290/71‑1 (MQ=255)
                                 ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                    <  1:303111/71‑1 (MQ=255)
                                  gagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGAc                                           >  1:2598203/1‑63 (MQ=12)
                                  gagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGAc                                           >  1:3861295/1‑63 (MQ=12)
                                  gagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCAGGAGTAATGAAACCGAATCCtttg                                              <  1:1869002/60‑2 (MQ=12)
                                    gaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAAc                                  >  1:5097356/1‑70 (MQ=255)
                                     aaGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt                                                <  1:3736081/55‑1 (MQ=255)
                                      aGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCAc                               <  1:4487975/71‑1 (MQ=255)
                                        tgtACGAACACATCTTTGCTGCCATCAGCCGGAGTa                                                                <  1:5235387/36‑1 (MQ=255)
                                        tgtACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg                                                             <  1:5460608/39‑1 (MQ=255)
                                        tgtAAGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAAccg                                                       <  1:1828006/45‑2 (MQ=2)
                                        tgtAAGAACACGTCATTGCTGCCGTCTTCCGGAGTAATGATACCGAATCCtt                                                <  1:3246410/52‑1 (MQ=255)
                                         gtACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATg                                                             >  1:1775226/1‑38 (MQ=18)
                                         gtACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATg                                                             >  1:1348652/1‑38 (MQ=18)
                                         gtACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg                                                             <  1:1027758/38‑1 (MQ=255)
                                          tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGaaa                                                          <  1:4490316/40‑2 (MQ=17)
                                          tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGaa                                                           <  1:5693580/39‑1 (MQ=17)
                                          tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTaa                           >  1:702603/1‑71 (MQ=255)
                                          tACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTaa                           >  1:2964801/1‑71 (MQ=255)
                                          tACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                                                              >  1:4421765/1‑36 (MQ=255)
                                          tACGAACACATCTTTGCTGCCATCAGCCGGAGTAAt                                                              >  1:3011408/1‑36 (MQ=255)
                                          tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                              >  1:4843779/1‑52 (MQ=255)
                                          tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                              >  1:3796598/1‑52 (MQ=255)
                                          tACGAACACATCTTTAGAGCCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCAtttta                           >  1:4551257/1‑69 (MQ=255)
                                           aCGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                                   >  1:2164978/1‑46 (MQ=255)
                                           aCGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                                   >  1:3856104/1‑46 (MQ=255)
                                           aCGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                                   >  1:1208086/1‑46 (MQ=255)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGa                                                            >  1:2823655/1‑36 (MQ=12)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGa                                                            >  1:2849694/1‑36 (MQ=12)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGa                                                            >  1:2873384/1‑36 (MQ=12)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGa                                                            >  1:3804661/1‑36 (MQ=12)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTaaa                                   >  1:3629146/1‑61 (MQ=9)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAAc                                  >  1:2731558/1‑62 (MQ=9)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAAc                                  >  1:2622983/1‑62 (MQ=9)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAcg                         >  1:4127071/1‑70 (MQ=12)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAcg                         >  1:192844/1‑70 (MQ=12)
                                            cGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAc                          <  1:3081630/70‑1 (MQ=12)
                                            cGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAcc                         <  1:4801092/71‑1 (MQ=255)
                                             gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                        >  1:4944556/1‑71 (MQ=255)
                                              aaCACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAcgtt                       <  1:4154567/71‑4 (MQ=12)
                                              aaCACATCTTTACTACCATCAACAGGAGAAATAAAGccga                                                      >  1:1525315/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)
                                                cacaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCAc                               <  1:3855277/61‑1 (MQ=255)
                                                  caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCt                   <  1:4781149/71‑1 (MQ=255)
                                                    tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg                                     <  1:5642545/51‑1 (MQ=255)
                                                     cTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc                                                   >  1:4205672/1‑36 (MQ=255)
                                                      tttGGTGCCATCAGCCGGAGTAATGAAGCCAAAACCt                                                 >  1:5021466/1‑37 (MQ=255)
                                                      tttGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtttgg                                             <  1:1025835/41‑3 (MQ=12)
                                                      tttGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtttg                                              <  1:4885998/40‑2 (MQ=12)
                                                      tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTg                       >  1:5654221/1‑63 (MQ=255)
                                                        tgctgcCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTaa               <  1:3695205/69‑1 (MQ=255)
                                                        tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                <  1:2164328/68‑1 (MQ=255)
                                                        tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                <  1:3678572/68‑1 (MQ=255)
                                                         gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt            >  1:5137044/1‑71 (MQ=255)
                                                         gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt            >  1:3595011/1‑71 (MQ=255)
                                                           tgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                <  1:1711364/65‑1 (MQ=255)
                                                                tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt     >  1:5457350/1‑71 (MQ=255)
                                                                  ttccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCa                                >  1:4018360/3‑42 (MQ=255)
                                                                  ttccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCa                                >  1:3968665/3‑42 (MQ=255)
                                                                  aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa                           >  1:5227390/1‑47 (MQ=255)
                                                                  aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCttt          >  1:4610099/1‑64 (MQ=255)
                                                                  aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt    >  1:3981886/1‑70 (MQ=255)
                                                                    ccGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTa                            <  1:3480573/44‑1 (MQ=255)
                                                                    ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAAc                          >  1:2053871/1‑46 (MQ=16)
                                                                    ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtcaa  <  1:4366361/70‑4 (MQ=255)
                                                                    ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtcaa  <  1:3476175/70‑4 (MQ=255)
                                                                    ccGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCTTTAATCTTAGACAtca   <  1:3136418/69‑3 (MQ=255)
                                                                     |                                                                    
AGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTG  >  minE/1211217‑1211354

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: