Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,220 C→A 11.8% L44L (CTG→CTT 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,2200CA11.8% 259.2 / 12.9 84L44L (CTG→CTTcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base A (6/4);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGA  >  minE/1211155‑1211289
                                                                 |                                                                     
aCGTTAACAGCTGCCGGACCTTTGTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGt                                                                  >  1:3797890/1‑71 (MQ=255)
aCGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGt                                                                  >  1:2123489/1‑71 (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:2308363/1‑71 (MQ=255)
           tGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCAt                                                       >  1:886128/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)
                     ttCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCa                                                        >  1:1354117/1‑60 (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:2751592/1‑40 (MQ=255)
                                 tGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGttt                                                                >  1:186154/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:156016/1‑62 (MQ=255)
                                       tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGaa                                    >  1:2571572/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:1687002/38‑1 (MQ=255)
                                         gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                                          <  1:1420971/38‑1 (MQ=255)
                                         gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                                          <  1:1179802/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:2657771/1‑71 (MQ=255)
                                           aCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca                       >  1:3192446/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:1544182/45‑1 (MQ=255)
                                                   ttCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCg                                         <  1:1548692/45‑1 (MQ=255)
                                                     cTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctg              <  1:1901325/70‑1 (MQ=255)
                                                      tGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACAt                      >  1:2409280/1‑61 (MQ=255)
                                                       gACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa           >  1:582459/1‑71 (MQ=255)
                                                        aCCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa           <  1:543206/70‑1 (MQ=255)
                                                        aCCTTCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGcc            <  1:1559666/69‑1 (MQ=25)
                                                         ccTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCAt          >  1:3633317/1‑70 (MQ=255)
                                                          cTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGaa                           <  1:3806519/52‑1 (MQ=255)
                                                          cTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCAt          >  1:1249309/1‑69 (MQ=255)
                                                            tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca                         >  1:223531/1‑52 (MQ=255)
                                                            tCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAac                          >  1:111609/1‑51 (MQ=255)
                                                            tCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGcc            >  1:440163/1‑65 (MQ=19)
                                                            tCAGCAAGAGTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTcttc      >  1:3353712/1‑68 (MQ=16)
                                                              aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg    <  1:1982487/71‑1 (MQ=255)
                                                              aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg    <  1:77269/71‑1 (MQ=255)
                                                              aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg    >  1:2785260/1‑71 (MQ=255)
                                                               gCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCgg   <  1:1518876/71‑1 (MQ=255)
                                                                ccAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGa  <  1:559992/71‑1 (MQ=255)
                                                                ccAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTACGAacact                       <  1:3563071/50‑2 (MQ=255)
                                                                ccAGAGATTTGTAACCATCGTTCTGGATAGCAGAGAAGTGTACGAacact                       <  1:2276615/50‑2 (MQ=255)
                                                                 |                                                                     
ACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGA  >  minE/1211155‑1211289

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: