Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,211,232 G→A 17.6% G40G (GGC→GGT 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,2320GA17.6% 197.1 / 29.9 85G40G (GGC→GGTcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base A (7/8);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCC  >  minE/1211165‑1211300
                                                                   |                                                                    
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)
        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)
                             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)
                                           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:77269/71‑1 (MQ=255)
                                                    aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg               >  1:2785260/1‑71 (MQ=255)
                                                    aGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg               <  1:1982487/71‑1 (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:3587239/1‑69 (MQ=255)
                                                              ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATg       >  1:689997/1‑69 (MQ=255)
                                                              ttGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa      <  1:3211159/70‑1 (MQ=255)
                                                                gAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg    <  1:2114138/70‑1 (MQ=255)
                                                                  aGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCGGCCGGAGTAATGAAGcc  >  1:2443707/1‑70 (MQ=255)
                                                                  aGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGcc  >  1:217784/1‑70 (MQ=255)
                                                                   gCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCtt                              <  1:2521295/41‑1 (MQ=255)
                                                                   gCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt           <  1:3169696/60‑1 (MQ=255)
                                                                   |                                                                    
CTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCC  >  minE/1211165‑1211300

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: