Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,104,168 C→T 20.4% F56F (TTC→TTT 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,1680CT20.4% 208.1 / 44.6 103F56F (TTC→TTTkbaZtagatose 6‑phosphate aldolase 1, kbaZ subunit
Reads supporting (aligned to +/- strand):  ref base C (38/44);  new base T (16/5);  total (54/49)
Fisher's exact test for biased strand distribution p-value = 1.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

CAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATT  >  minE/2104098‑2104238
                                                                      |                                                                      
cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                        <  1:1891971/71‑1 (MQ=255)
cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                        >  1:4476016/1‑71 (MQ=255)
cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                        >  1:1552048/1‑71 (MQ=255)
cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                        <  1:2613918/71‑1 (MQ=255)
  aTCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcgg                                                                      <  1:5047803/71‑1 (MQ=255)
   tCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggc                                                                     <  1:1224130/71‑1 (MQ=255)
   tCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggc                                                                     >  1:3645520/1‑71 (MQ=255)
     cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcgg                                                                   >  1:3361451/1‑71 (MQ=255)
     cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcg                                                                    <  1:2947930/70‑1 (MQ=255)
     cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcg                                                                    <  1:2347729/70‑1 (MQ=255)
      gCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGt                                                                  <  1:390094/71‑1 (MQ=255)
         actccAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                                >  1:4658344/4‑70 (MQ=255)
              aGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCgg                                                          <  1:4261789/71‑1 (MQ=255)
                tGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGaa                                                        <  1:2162299/71‑1 (MQ=255)
                  caatccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGTTAGACCAGTTCGGCGGTTATACCGGAAt                                                       <  1:2460635/70‑1 (MQ=255)
                  caaaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                                       >  1:1291697/1‑70 (MQ=255)
                   aaaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGa                                                     <  1:961907/71‑1 (MQ=255)
                    aaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                                >  1:1162522/1‑59 (MQ=255)
                     accaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcg                                                                       >  1:268734/1‑51 (MQ=255)
                     accaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcg                                                   >  1:3833679/1‑71 (MQ=255)
                       caaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgcc                                                 >  1:4833883/1‑71 (MQ=255)
                          acGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtata                                                              >  1:5337673/1‑55 (MQ=255)
                           cGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgcccgcc                                             >  1:2849938/1‑71 (MQ=255)
                           cGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgcccgcc                                             >  1:2964781/1‑71 (MQ=255)
                            gCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCccgccg                                            <  1:5108881/71‑1 (MQ=255)
                              cGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                                          >  1:688022/1‑71 (MQ=255)
                               gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                         <  1:4663429/71‑1 (MQ=255)
                               gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                         <  1:5237288/71‑1 (MQ=255)
                                ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGt                                                                  <  1:5040655/45‑1 (MQ=255)
                                ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtata                                                              <  1:5157905/49‑1 (MQ=255)
                                ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                                <  1:3134208/47‑1 (MQ=255)
                                ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAttt                                        <  1:4517395/71‑1 (MQ=255)
                                ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                                         <  1:571206/70‑1 (MQ=255)
                                 tgcTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTta                                                                <  1:5381246/44‑1 (MQ=255)
                                 tgcTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGActtt                                       >  1:1452742/3‑71 (MQ=255)
                                 tACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtata                                                              <  1:4972172/48‑1 (MQ=255)
                                  aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcg                                                                    >  1:3768371/1‑41 (MQ=255)
                                  aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                        >  1:4137225/1‑37 (MQ=255)
                                  aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                                        >  1:5181723/1‑37 (MQ=255)
                                  aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCAGCCGATTTTc                                      >  1:3781488/1‑71 (MQ=255)
                                   cTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTtata                                                              >  1:4144947/1‑46 (MQ=255)
                                   cTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTcg                                     >  1:4115025/1‑71 (MQ=255)
                                    tGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGActt                                        <  1:4642488/67‑1 (MQ=255)
                                     gATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTcgcg                                   >  1:1766241/1‑71 (MQ=255)
                                      aTTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATAcc                                                            <  1:4982771/45‑1 (MQ=255)
                                      aTTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                                       <  1:4570683/50‑1 (MQ=255)
                                      aTTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                                       <  1:1889668/50‑1 (MQ=255)
                                       ttGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaa                                 >  1:2560392/1‑69 (MQ=255)
                                       ttGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaa                                 >  1:2350933/1‑69 (MQ=255)
                                       ttGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                                       <  1:2531918/49‑1 (MQ=255)
                                        tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCgg                                                          >  1:2220899/1‑45 (MQ=255)
                                        tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                                >  1:3330137/1‑68 (MQ=255)
                                        tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                                >  1:4235724/1‑68 (MQ=255)
                                        tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat                                >  1:5112299/1‑68 (MQ=255)
                                        tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtat                                                               <  1:4668275/40‑1 (MQ=255)
                                        tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACTCCCGCCGATTTTCGCgg                                  >  1:189999/1‑69 (MQ=255)
                                         gAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                                >  1:2316795/1‑38 (MQ=255)
                                           aGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGTAATGACACCGGCAGACTTTCGCGAAtttg                             <  1:1525802/71‑1 (MQ=255)
                                           aGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                                >  1:1416292/1‑36 (MQ=255)
                                           aGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                                                >  1:1862434/1‑36 (MQ=255)
                                             cAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGaa                                                        >  1:1464212/1‑42 (MQ=255)
                                             cAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGaa                                                        >  1:2724506/1‑42 (MQ=255)
                                              aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGa                                                     >  1:491660/1‑44 (MQ=255)
                                              aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGa                                                     >  1:3224518/1‑44 (MQ=255)
                                              aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTCGCGAAtttgttt                          <  1:4287668/71‑1 (MQ=255)
                                              aaCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                          <  1:3228111/71‑1 (MQ=255)
                                               aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgc                                                  >  1:982795/1‑46 (MQ=255)
                                               aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                                          <  1:959601/54‑1 (MQ=255)
                                               aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                                          <  1:4117836/54‑1 (MQ=255)
                                               aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgtttg                         <  1:1969733/71‑1 (MQ=255)
                                               aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                            >  1:3450725/1‑68 (MQ=255)
                                               aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgt                            >  1:4174387/1‑68 (MQ=255)
                                                ccTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                                       <  1:5025624/40‑1 (MQ=255)
                                                ccTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                                       <  1:5441795/40‑1 (MQ=255)
                                                 cTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg                             <  1:5583957/65‑1 (MQ=255)
                                                  tCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgccc                                                >  1:5458015/1‑45 (MQ=255)
                                                  tCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCt                                >  1:5653725/1‑61 (MQ=255)
                                                  tCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCt                                >  1:220938/1‑61 (MQ=255)
                                                      aTCAGGTAGACCAGTTCGGCGGTTATACCGGAATGa                                                     >  1:4285320/1‑36 (MQ=255)
                                                      aTCAGGTAGACCAGTTCGGCGGTTATACCGGAATGa                                                     >  1:3025590/1‑36 (MQ=255)
                                                      aTCAGGTAGACCAGTTCGGCGGTTATACCGGAATGa                                                     >  1:953726/1‑36 (MQ=255)
                                                       tCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg                             <  1:1641440/59‑1 (MQ=255)
                                                       tCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg                             <  1:5314207/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:1520900/1‑71 (MQ=255)
                                                         aGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGcc               >  1:4632385/1‑71 (MQ=255)
                                                           gTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                          <  1:1211817/58‑1 (MQ=255)
                                                            tAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGa             >  1:5044543/1‑70 (MQ=255)
                                                              gACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTc          >  1:4728619/1‑71 (MQ=255)
                                                               aCCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGc                                 <  1:5684670/47‑1 (MQ=255)
                                                                  aGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTg      <  1:1085077/71‑1 (MQ=255)
                                                                   gTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                          <  1:3600970/50‑1 (MQ=255)
                                                                   gTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt                          <  1:1050824/50‑1 (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)
                                                                      |                                                                      
CAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATT  >  minE/2104098‑2104238

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: