Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,229 G→A 12.9% F41F (TTC→TTT 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,2290GA12.9% 225.7 / 20.4 85F41F (TTC→TTTcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base A (7/4);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

AACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAG  >  minE/1211160‑1211298
                                                                     |                                                                     
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)
  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)
             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:725742/70‑1 (MQ=255)
                  cTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATa                                                     <  1:3792426/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)
                            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:1420971/38‑1 (MQ=255)
                                    gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                                                   <  1:1179802/38‑1 (MQ=255)
                                    gAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGc                                                                   <  1:1687002/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)
                                                             agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAg          >  1:1540361/1‑70 (MQ=255)
                                                             agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt         <  1:2122453/71‑1 (MQ=255)
                                                               agTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa       >  1:629527/1‑71 (MQ=255)
                                                               agTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTaa       >  1:2564493/1‑71 (MQ=255)
                                                                 ttttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg     <  1:1300803/71‑1 (MQ=255)
                                                                   ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg     >  1:689997/1‑69 (MQ=255)
                                                                   ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg     >  1:3587239/1‑69 (MQ=255)
                                                                   ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa    <  1:3211159/70‑1 (MQ=255)
                                                                     gAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg  <  1:2114138/70‑1 (MQ=255)
                                                                     |                                                                     
AACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAG  >  minE/1211160‑1211298

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: