Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,104,192 C→G 12.7% P64P (CCC→CCG kbaZ → tagatose 6‑phosphate aldolase 1, kbaZ subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,104,1920CG12.7% 243.7 / 21.4 95P64P (CCC→CCGkbaZtagatose 6‑phosphate aldolase 1, kbaZ subunit
Reads supporting (aligned to +/- strand):  ref base C (36/46);  new base G (8/4);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 105 total aligned reads displayed.

CGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTG  >  minE/2104125‑2104261
                                                                   |                                                                     
cGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgcccgcc                                                                    >  1:2964781/1‑71 (MQ=255)
 gCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCccgccg                                                                   <  1:5108881/71‑1 (MQ=255)
   cGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                                                                 >  1:688022/1‑71 (MQ=255)
    gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                                                <  1:5237288/71‑1 (MQ=255)
    gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                                                <  1:4670184/71‑1 (MQ=255)
    gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                                                <  1:4663429/71‑1 (MQ=255)
     ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAttt                                                               <  1:4517395/71‑1 (MQ=255)
     ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                                                <  1:571206/70‑1 (MQ=255)
      tgcTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGActtt                                                              >  1:1452742/3‑71 (MQ=255)
       aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCAGCCGATTTTc                                                             >  1:3781488/1‑71 (MQ=255)
        cTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTcg                                                            >  1:4115025/1‑71 (MQ=255)
         tGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGActt                                                               <  1:4642488/67‑1 (MQ=255)
           aTTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcg                                                          >  1:3244603/1‑70 (MQ=255)
            ttGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaa                                                        >  1:2350933/1‑69 (MQ=255)
            ttGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaa                                                        >  1:2560392/1‑69 (MQ=255)
            ttGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGc                                                        >  1:5301465/1‑71 (MQ=255)
             tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                                                       >  1:4235724/1‑68 (MQ=255)
             tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                                                       >  1:3330137/1‑68 (MQ=255)
             tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                                                       >  1:5112299/1‑68 (MQ=255)
             tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                                                       >  1:2493561/1‑68 (MQ=255)
             tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACTCCCGCCGATTTTCGCgg                                                         >  1:189999/1‑69 (MQ=255)
                aGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGTAATGACACCGGCAGACTTTCGCGAAtttg                                                    <  1:1525802/71‑1 (MQ=255)
                   aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTCGCGAAtttgttt                                                 <  1:4287668/71‑1 (MQ=255)
                   aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTCGCGAAtttgttt                                                 <  1:5011078/71‑1 (MQ=255)
                   aaCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                                                 <  1:5089189/71‑1 (MQ=255)
                   aaCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                                                 <  1:3228111/71‑1 (MQ=255)
                    aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                                                                 <  1:4117836/54‑1 (MQ=255)
                    aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                                                                 <  1:959601/54‑1 (MQ=255)
                    aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgtttg                                                <  1:1969733/71‑1 (MQ=255)
                    aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                                                   >  1:4174387/1‑68 (MQ=255)
                    aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                                                   >  1:3450725/1‑68 (MQ=255)
                     ccTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTttgtttgt                                               <  1:4330402/71‑1 (MQ=255)
                      cTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg                                                    <  1:5583957/65‑1 (MQ=255)
                      cTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTc                                              <  1:3011683/71‑1 (MQ=255)
                       tCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgccc                                                                       >  1:5458015/1‑45 (MQ=255)
                       tCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCt                                                       >  1:5653725/1‑61 (MQ=255)
                       tCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCt                                                       >  1:220938/1‑61 (MQ=255)
                            tCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg                                                    <  1:5314207/59‑1 (MQ=255)
                            tCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg                                                    <  1:1641440/59‑1 (MQ=255)
                             cAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTttgtttgt                                               >  1:4656143/1‑63 (MQ=255)
                             cAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGc                                       >  1:3287018/1‑71 (MQ=255)
                              aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTcgc                                                           >  1:4649050/1‑50 (MQ=255)
                              aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGcc                                      >  1:2577266/1‑71 (MQ=255)
                              aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGcc                                      >  1:4632385/1‑71 (MQ=255)
                              aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGcc                                      >  1:1520900/1‑71 (MQ=255)
                                gTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                                                 <  1:1211817/58‑1 (MQ=255)
                                   gACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTc                                 >  1:4728619/1‑71 (MQ=255)
                                    aCCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGc                                                        <  1:4646609/47‑1 (MQ=255)
                                    aCCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGc                                                        <  1:5684670/47‑1 (MQ=255)
                                       aGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTg                             <  1:1085077/71‑1 (MQ=255)
                                        gTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                                                 <  1:1050824/50‑1 (MQ=255)
                                        gTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                                                 <  1:3600970/50‑1 (MQ=255)
                                        gTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCTCCGACTCGTTg                             >  1:5497384/1‑70 (MQ=255)
                                         ttCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                                                 <  1:2550838/49‑1 (MQ=255)
                                           cggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg                                                    <  1:1064807/44‑1 (MQ=255)
                                           cggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCt                                          >  1:2059892/1‑54 (MQ=255)
                                           cggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAAtt                         >  1:24378/1‑71 (MQ=255)
                                            ggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAAttt                        <  1:3279030/71‑1 (MQ=255)
                                            ggcggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAAtt                         <  1:1397652/70‑1 (MQ=255)
                                              cggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTc                                              <  1:1322977/47‑1 (MQ=255)
                                              cggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTc                                              <  1:3511699/47‑1 (MQ=255)
                                              cggTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTc                                              <  1:654665/47‑1 (MQ=255)
                                                gTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCg                    >  1:1562606/1‑71 (MQ=255)
                                                 ttATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCt                                          >  1:1570889/1‑48 (MQ=255)
                                                 ttATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGc                   <  1:110540/71‑1 (MQ=255)
                                                 ttATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGc                   <  1:1796869/71‑1 (MQ=255)
                                                 ttATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGc                   <  1:4609439/71‑1 (MQ=255)
                                                 ttATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGc                   <  1:1803659/71‑1 (MQ=255)
                                                  tataCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCg                    >  1:1695852/1‑69 (MQ=255)
                                                  tataCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCg                    >  1:2805223/1‑69 (MQ=255)
                                                   ataCCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                                                 >  1:875173/1‑39 (MQ=255)
                                                   ataCCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                                                 >  1:217082/1‑39 (MQ=255)
                                                   ataCCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                                                 >  1:3887035/1‑39 (MQ=255)
                                                   ataCCGGAATGACGCCCGCCGATTTTCGCGGCtttgtt                                                  >  1:140867/1‑38 (MQ=255)
                                                   ataCCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                                                   >  1:898032/1‑37 (MQ=255)
                                                   ataCCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                                                   >  1:4929813/1‑37 (MQ=255)
                                                   ataCCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                                                   >  1:1132465/1‑37 (MQ=255)
                                                   ataCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAg                 >  1:3104457/1‑71 (MQ=255)
                                                    taCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAg                 <  1:4967494/70‑1 (MQ=255)
                                                     aCCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                                                   <  1:1490198/35‑1 (MQ=255)
                                                     aCCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                                                   <  1:164171/35‑1 (MQ=255)
                                                     aCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTc                                              <  1:4168486/40‑1 (MQ=255)
                                                     aCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTc                                              <  1:4887608/40‑1 (MQ=255)
                                                      ccGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTccc                     <  1:2348319/64‑1 (MQ=255)
                                                      ccGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGAt              >  1:2578364/1‑71 (MQ=255)
                                                       cGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGa                                    >  1:3034888/1‑48 (MQ=255)
                                                         gAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGc                   <  1:1010246/63‑1 (MQ=255)
                                                           aTGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGt          >  1:5342419/1‑70 (MQ=255)
                                                            tGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTg                             <  1:614634/50‑1 (MQ=255)
                                                            tGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAg                 <  1:4775459/62‑1 (MQ=255)
                                                            tGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAg                 <  1:4053138/62‑1 (MQ=255)
                                                            tGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAg                 <  1:2545375/62‑1 (MQ=255)
                                                              aCGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGa       <  1:2370794/70‑1 (MQ=255)
                                                              aCGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGa       <  1:2223626/70‑1 (MQ=255)
                                                              aCGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGAt      <  1:1293599/71‑1 (MQ=255)
                                                                gcccgccGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTc    <  1:3129834/71‑1 (MQ=255)
                                                                gcccgccGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTc    <  1:3123312/71‑1 (MQ=255)
                                                                 cccgccGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCt   >  1:2485545/1‑71 (MQ=255)
                                                                  ccgccgATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTg  >  1:3719077/1‑71 (MQ=255)
                                                                  ccgccgATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTg  >  1:949989/1‑71 (MQ=255)
                                                                   |                                                                     
CGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTG  >  minE/2104125‑2104261

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: