Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,240 C→T 10.8% G38S (GGT→AGT) ‡ 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,2400CT10.8% 214.2 / 13.7 74G38S (GGT→AGT) ‡cspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base T (2/6);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTA  >  minE/1211173‑1211310
                                                                   |                                                                      
ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                                                     <  1:2079596/71‑1 (MQ=255)
ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                                                     >  1:876310/1‑71 (MQ=255)
   ttCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                                   >  1:561217/1‑70 (MQ=255)
    tCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                                   >  1:2545893/1‑69 (MQ=255)
    tCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                                   >  1:226970/1‑69 (MQ=255)
        gCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAg                                                                <  1:2579163/68‑1 (MQ=255)
         ccGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGga                                                            <  1:408955/71‑1 (MQ=255)
         ccGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGga                                                            <  1:3319079/71‑1 (MQ=255)
          cGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                           <  1:832819/71‑1 (MQ=255)
            tCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                                                     <  1:3574691/59‑1 (MQ=255)
            tCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                           >  1:3493370/1‑69 (MQ=255)
              cTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                        >  1:302072/1‑54 (MQ=255)
              cTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                        >  1:762135/1‑54 (MQ=255)
              cTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgt                                                       <  1:3163020/71‑1 (MQ=255)
               gtgaTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAgtg                                                      >  1:591151/4‑71 (MQ=255)
                gAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                                                     <  1:3904064/55‑1 (MQ=255)
                 aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                      <  1:2205189/69‑1 (MQ=255)
                  aTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCagag                                                           <  1:3136989/63‑1 (MQ=255)
                   tttCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACg                                                  <  1:2421040/71‑1 (MQ=255)
                     tCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGa                                                 <  1:803385/70‑1 (MQ=255)
                     tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCATGGATAGCGGAGAAGTGTACGaa                                                >  1:629314/1‑71 (MQ=255)
                        aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac                                             <  1:3285844/71‑1 (MQ=255)
                        aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca                                              >  1:1987656/1‑70 (MQ=255)
                           tCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac                                             <  1:2279969/68‑1 (MQ=255)
                            cAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCt                                         >  1:149205/1‑71 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                        >  1:3603370/1‑39 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGaa                                                         >  1:1486940/1‑54 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac                                             >  1:2877128/1‑66 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca                                              <  1:2092200/65‑1 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                        >  1:3693883/1‑71 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                        >  1:2726902/1‑71 (MQ=255)
                             aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                        <  1:1291169/71‑1 (MQ=255)
                              aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCttt                                       <  1:3964724/71‑1 (MQ=255)
                              aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCttt                                       <  1:818930/71‑1 (MQ=255)
                               cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                        >  1:1666898/1‑37 (MQ=255)
                               cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                        >  1:666624/1‑37 (MQ=255)
                               cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                        >  1:558073/1‑37 (MQ=255)
                                gTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTa                                                    <  1:3536567/56‑1 (MQ=255)
                                gTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                     >  1:3809707/1‑71 (MQ=255)
                                 ttCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgct                                    >  1:2979536/1‑71 (MQ=255)
                                  tCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac                                             >  1:770969/1‑61 (MQ=255)
                                   cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                                   >  1:3330496/1‑38 (MQ=255)
                                   cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctgc                                  >  1:941143/1‑71 (MQ=255)
                                   cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctg                                   >  1:3447212/1‑70 (MQ=255)
                                   cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctg                                   <  1:3052052/70‑1 (MQ=255)
                                    tGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTACGaa                                                <  1:1301994/56‑1 (MQ=255)
                                    tGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTACGaa                                                <  1:1278779/56‑1 (MQ=255)
                                     gACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAg                                                                >  1:1740228/1‑39 (MQ=255)
                                      aCCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                <  1:3725181/70‑1 (MQ=255)
                                      aCCTTCAGCCAGAGTTTTAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                                <  1:2341625/69‑1 (MQ=255)
                                          tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                           >  1:3223482/1‑71 (MQ=255)
                                          tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                           >  1:81415/1‑71 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATc                                          >  1:1257206/1‑55 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                     >  1:191753/1‑60 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                          >  1:2992955/1‑71 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                          >  1:2591574/1‑71 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                          >  1:3032835/1‑71 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                          >  1:3071137/1‑71 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc                          >  1:2272288/1‑71 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                           <  1:864049/70‑1 (MQ=255)
                                           cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc                           <  1:1122496/70‑1 (MQ=255)
                                             gCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                      <  1:873024/41‑1 (MQ=255)
                                               cAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCt                                         >  1:1740751/1‑52 (MQ=255)
                                                agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                     <  1:2435086/71‑1 (MQ=255)
                                                agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                     <  1:436698/71‑1 (MQ=255)
                                                 gagTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGt                     >  1:3538257/1‑70 (MQ=255)
                                                    ttttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATc                                          <  1:1960922/46‑1 (MQ=255)
                                                    ttttAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAacacgt                                           <  1:477609/45‑3 (MQ=255)
                                                    ttttAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCtt                                        <  1:2600963/48‑1 (MQ=255)
                                                      ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgct                                    >  1:2327533/1‑50 (MQ=255)
                                                      ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                     >  1:1325835/1‑49 (MQ=255)
                                                      ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                     >  1:2066364/1‑49 (MQ=255)
                                                       tGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGaa               <  1:4144428/70‑1 (MQ=255)
                                                       tGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg              >  1:1055842/1‑71 (MQ=255)
                                                        gAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCa                             <  1:294587/55‑1 (MQ=255)
                                                               ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCAt                               <  1:4052745/46‑1 (MQ=255)
                                                               ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa                <  1:3547024/61‑1 (MQ=255)
                                                                tACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                                        <  1:3053343/36‑1 (MQ=255)
                                                                 aCCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc                                     >  1:2597496/1‑38 (MQ=255)
                                                                 aCCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAg                            >  1:2499169/1‑47 (MQ=255)
                                                                  cccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt    <  1:519036/70‑1 (MQ=255)
                                                                   ccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCTTCAGCCGGAGTAATGAAGCCAAAACCTTTa  <  1:1024990/71‑1 (MQ=255)
                                                                   ccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa  <  1:174806/71‑1 (MQ=255)
                                                                   |                                                                      
CCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTA  >  minE/1211173‑1211310

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: