Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,207 T→G 8.0% N49H (AAC→CAC) ‡ 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,2070TG 8.0% 183.8 / 6.6 74N49H (AAC→CAC) ‡cspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base G (3/3);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGC  >  minE/1211140‑1211275
                                                                   |                                                                    
gATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCt                                                                    <  1:879768/70‑1 (MQ=255)
gATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCt                                                                    <  1:748010/70‑1 (MQ=255)
gATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCt                                                                    <  1:1105075/70‑1 (MQ=255)
    aGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCt                                                               <  1:3007819/71‑1 (MQ=255)
     gATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCtt                                                              <  1:935294/71‑1 (MQ=255)
       tAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCtt                                                              >  1:2449383/1‑69 (MQ=255)
        aGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGa                                                                  >  1:3587162/1‑64 (MQ=255)
          cTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGa                                                                  <  1:853926/62‑1 (MQ=255)
               aCGTTAACAGCTGCCGGACCTTTGTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGt                                                    >  1:3797890/1‑71 (MQ=255)
               aCGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGt                                                    >  1:2123489/1‑71 (MQ=255)
                   tAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCa                                                            <  1:1509052/59‑1 (MQ=255)
                    aaCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTg                                                <  1:873703/70‑1 (MQ=255)
                    aaCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTg                                                <  1:1399248/70‑1 (MQ=255)
                     aCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTg                                                >  1:1160875/1‑69 (MQ=255)
                     aCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGaa                                              <  1:359650/71‑1 (MQ=255)
                      cAGCTGCCGGACCTTTCTGGCCTTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGtt                                                   <  1:3585270/65‑1 (MQ=255)
                      cAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGa                                               >  1:1860093/1‑69 (MQ=255)
                         cTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCa                                          >  1:1055219/1‑71 (MQ=255)
                          tGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAt                                         >  1:886128/1‑71 (MQ=255)
                          tGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAt                                         >  1:2308363/1‑71 (MQ=255)
                            ccGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAt                                         <  1:3759316/69‑1 (MQ=255)
                            ccGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTa                                       >  1:2699652/1‑71 (MQ=255)
                            ccGGACCTTTCTGGCCGTCCTGAATCTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAtt                                        <  1:531293/70‑1 (MQ=255)
                            ccGGACCTTTCTGGACGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTa                                       <  1:2236277/71‑1 (MQ=255)
                              ggACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCa                                          >  1:2800256/1‑66 (MQ=255)
                               gACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCaga                                                      <  1:1283729/53‑1 (MQ=255)
                                 ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGaa                                              >  1:2902704/1‑59 (MQ=255)
                                 ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                  >  1:778458/1‑71 (MQ=255)
                                 ccTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                  >  1:899542/1‑71 (MQ=255)
                                   tttCTGGCCGTCCTGAATTTCGAACTCAACGTTCt                                                                    <  1:1051872/35‑1 (MQ=255)
                                    ttCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCa                                          >  1:1354117/1‑60 (MQ=255)
                                     tCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTc                                                             >  1:3051221/1‑40 (MQ=255)
                                     tCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                <  1:690552/69‑1 (MQ=255)
                                      cTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATa                              <  1:3792426/70‑1 (MQ=255)
                                      cTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATa                              <  1:725742/70‑1 (MQ=255)
                                         gCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGc                            <  1:3019907/69‑1 (MQ=255)
                                         gCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCg                           <  1:2478583/70‑1 (MQ=255)
                                            gTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                        >  1:1478149/1‑70 (MQ=255)
                                               cTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGcc                                           <  1:3159318/48‑1 (MQ=255)
                                                tGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGttt                                                  >  1:186154/1‑40 (MQ=255)
                                                tGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGttt                                                  >  1:2751592/1‑40 (MQ=255)
                                                gtgaTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAgtg                   <  1:477945/68‑1 (MQ=255)
                                                   aTTTCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTAc                >  1:1059989/1‑71 (MQ=255)
                                                    tttCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                  >  1:217490/1‑52 (MQ=255)
                                                     ttCGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTa                 <  1:706655/68‑1 (MQ=255)
                                                      tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTg                                                <  1:1915725/36‑1 (MQ=255)
                                                      tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGaa                      >  1:2571572/1‑62 (MQ=255)
                                                      tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGaa                      >  1:156016/1‑62 (MQ=255)
                                                      tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGa              <  1:176980/70‑1 (MQ=255)
                                                       cGAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGaa             >  1:1637821/1‑70 (MQ=255)
                                                        gAACTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAac            <  1:1652031/70‑1 (MQ=255)
                                                        gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                            <  1:1179802/38‑1 (MQ=255)
                                                        gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                            <  1:1687002/38‑1 (MQ=255)
                                                        gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                            <  1:1420971/38‑1 (MQ=255)
                                                         aaCTCTACGCGCTGACCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAacac          >  1:372179/1‑71 (MQ=255)
                                                         aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                  >  1:1054772/1‑47 (MQ=255)
                                                         aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac          <  1:24290/71‑1 (MQ=255)
                                                         aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca           >  1:2495753/1‑70 (MQ=255)
                                                          aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAacac          <  1:3755859/70‑1 (MQ=255)
                                                          aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca         >  1:3192446/1‑71 (MQ=255)
                                                          aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca         >  1:2657771/1‑71 (MQ=255)
                                                           cTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                     <  1:2087680/42‑1 (MQ=255)
                                                           cTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCg                           <  1:933161/52‑1 (MQ=255)
                                                           cTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca         >  1:3527977/1‑70 (MQ=255)
                                                            tCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGaa                      <  1:605378/56‑1 (MQ=255)
                                                              aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                     >  1:1478301/1‑39 (MQ=255)
                                                              aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGTAGAAGTGTACGAACACATCt      <  1:2715024/70‑1 (MQ=255)
                                                              aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt     >  1:3407707/1‑71 (MQ=255)
                                                              aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGAGTACGAACACATCt      >  1:103504/1‑70 (MQ=255)
                                                               aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAtt                                        <  1:269154/35‑1 (MQ=255)
                                                               aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAGGTGTACGAACACATCttt    <  1:947240/71‑1 (MQ=255)
                                                                cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACTCATCTTtg   >  1:1761059/1‑71 (MQ=255)
                                                                cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtg   >  1:3923645/1‑71 (MQ=255)
                                                                 tttCTGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTa                 <  1:2721220/55‑1 (MQ=255)
                                                                 gTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgc  >  1:637407/1‑71 (MQ=255)
                                                                  ttCTGACCTTCGTCCAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTACGAacactga       <  1:621917/65‑4 (MQ=255)
                                                                  ttCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCg                           <  1:2999891/45‑1 (MQ=255)
                                                                  ttCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCg                           <  1:1548692/45‑1 (MQ=255)
                                                                  ttCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCg                           <  1:1544182/45‑1 (MQ=255)
                                                                   |                                                                    
GATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGC  >  minE/1211140‑1211275

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: