Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,281,121 A→G 12.6% F67S (TTT→TCT)  fabR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,281,1210AG12.6% 309.2 / 31.1 111F67S (TTT→TCT) fabRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (58/39);  new base G (3/11);  total (61/50)
Fisher's exact test for biased strand distribution p-value = 9.23e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 124 total aligned reads displayed.

TAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTGCGTT  >  minE/2281054‑2281187
                                                                   |                                                                  
tAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATaaa                                                                    >  1:362282/1‑68 (MQ=255)
     ttAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGt                                                             <  1:1883399/70‑1 (MQ=255)
     ttAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGt                                                             <  1:1944587/70‑1 (MQ=255)
     ttAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGt                                                             <  1:75722/70‑1 (MQ=255)
     ttAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTg                                                            >  1:1162278/1‑71 (MQ=255)
      tAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTgg                                                           >  1:4056472/1‑71 (MQ=255)
       aaaCCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAaga                                                                >  1:2678681/1‑65 (MQ=255)
       aaaCCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAaga                                                                >  1:3605844/1‑65 (MQ=255)
       aaaCCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAaga                                                                >  1:58302/1‑65 (MQ=255)
       aaaCCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTgg                                                           >  1:2500561/1‑70 (MQ=255)
        aaCCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTggg                                                          <  1:2760253/70‑1 (MQ=255)
          ccGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGc                                                       <  1:2446103/71‑1 (MQ=255)
          ccGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGc                                                       <  1:1428683/71‑1 (MQ=255)
          ccGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGc                                                       <  1:2111345/71‑1 (MQ=255)
            gCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCa                                                      >  1:3690632/1‑70 (MQ=255)
            gCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCa                                                      >  1:2712200/1‑70 (MQ=255)
             ctctCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGTGGGAGCAAt                                                    >  1:3618362/1‑71 (MQ=255)
             ctctCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAaga                                                                >  1:834149/1‑59 (MQ=255)
             ctctCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAaga                                                                >  1:814609/1‑59 (MQ=255)
             ctctCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAaga                                                                >  1:132026/1‑59 (MQ=255)
             ctctCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAaga                                                                >  1:1332810/1‑59 (MQ=255)
              tctcATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAAt                                                    <  1:1719115/70‑1 (MQ=255)
                tcatcaACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATaaa                                                                    <  1:578627/52‑1 (MQ=255)
                  atcaACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGTGGGAGCAATGCCCg                                               >  1:2189916/1‑71 (MQ=255)
                   tcaACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATaaaaa                                                                  <  1:2650614/51‑1 (MQ=255)
                    caACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCt                                             >  1:4444369/1‑71 (MQ=255)
                     aaCCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTg                                                            >  1:2685463/1‑55 (MQ=255)
                     aaCCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCtt                                            >  1:607824/1‑71 (MQ=255)
                      aCCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCtt                                            >  1:1276232/1‑70 (MQ=255)
                      aCCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTc                                           <  1:3045094/71‑1 (MQ=255)
                        cATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGTGGGAGCAATGCCCGCTTCAc                                         <  1:4546479/71‑1 (MQ=255)
                        cATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGTGGGAGCAATGCCCGCTTCAc                                         <  1:4541158/71‑1 (MQ=255)
                        cATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGTGGGAGCAATGCCCGCTTCAc                                         <  1:4642730/71‑1 (MQ=255)
                        cATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCAc                                         <  1:2226269/71‑1 (MQ=255)
                         aTGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGTGGGa                                                         >  1:1124368/1‑54 (MQ=255)
                         aTGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGa                                                         <  1:1127329/54‑1 (MQ=255)
                         aTGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGa                                                         <  1:3437122/54‑1 (MQ=255)
                         aTGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGa                                                         <  1:4089359/54‑1 (MQ=255)
                         aTGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCAcg                                        >  1:3628391/1‑71 (MQ=255)
                          tGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCTGATaaa                                                                    <  1:607806/42‑1 (MQ=255)
                           ggTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGc                                                       >  1:3098197/1‑54 (MQ=255)
                           ggTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGc                                                       >  1:3010820/1‑54 (MQ=255)
                            gTCAGACGCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTc                                           <  1:2716995/65‑1 (MQ=255)
                            gTCAGACCCATTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTc                                           <  1:1043808/65‑1 (MQ=255)
                            gTCAGACCCATTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTc                                           <  1:1049450/65‑1 (MQ=255)
                            gTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATaaaaa                                                                  >  1:2972096/1‑42 (MQ=255)
                            gTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATaaaaa                                                                  >  1:415356/1‑42 (MQ=255)
                            gTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATaaaaa                                                                  >  1:259987/1‑42 (MQ=255)
                            gTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATaaaaa                                                                  >  1:2006761/1‑42 (MQ=255)
                            gTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATaaaaa                                                                  >  1:2635020/1‑42 (MQ=255)
                            gTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTc                                           <  1:226677/65‑1 (MQ=255)
                              cAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCAc                                         >  1:3773836/1‑65 (MQ=255)
                              cAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCAc                                         >  1:3124256/1‑65 (MQ=255)
                                gACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACt                                 >  1:44703/1‑71 (MQ=255)
                                gACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACt                                 >  1:1033887/1‑71 (MQ=255)
                                gACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACt                                 >  1:2701001/1‑71 (MQ=255)
                                 aCCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAag                                                                 <  1:4160218/38‑1 (MQ=255)
                                  cccAGTTCGTCTACGTCGCGGAAATGCCGATAAGAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTc                               <  1:632778/71‑1 (MQ=255)
                                  cccAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGGGGGAGCAATGCCCGCTTCACGCGCCACtt                                <  1:4258837/70‑1 (MQ=255)
                                      gTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACg                            >  1:351109/1‑70 (MQ=255)
                                      gTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGc                           <  1:1202459/71‑1 (MQ=255)
                                        tCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCaa                                                     >  1:2522040/1‑43 (MQ=255)
                                        tCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCa                          >  1:1705212/1‑70 (MQ=255)
                                        tCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCa                          >  1:586299/1‑70 (MQ=255)
                                        tCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCa                          >  1:4098170/1‑70 (MQ=255)
                                        tCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCa                          >  1:2557487/1‑70 (MQ=255)
                                         cGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCaa                                                     >  1:4494562/1‑42 (MQ=255)
                                           tCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACt                      <  1:554570/71‑1 (MQ=255)
                                           tCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACt                      <  1:260422/71‑1 (MQ=255)
                                            cTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCtt                                            <  1:2113926/48‑1 (MQ=255)
                                            cTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCaaa                        >  1:2240063/1‑68 (MQ=255)
                                            cTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACt                      <  1:4700569/70‑1 (MQ=255)
                                             tACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCaaa                        >  1:3211474/1‑67 (MQ=255)
                                             tACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCaaa                        >  1:1431776/1‑67 (MQ=255)
                                               cGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTc                     <  1:4417231/68‑1 (MQ=255)
                                                  cgcgGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTg               >  1:2660571/1‑71 (MQ=255)
                                                  cgcgGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTg               >  1:4227022/1‑71 (MQ=255)
                                                  cgcgGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTg               >  1:349384/1‑71 (MQ=255)
                                                      gAAATGCCGATAAGAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGa           <  1:3731383/71‑1 (MQ=255)
                                                      gAAATGCCGATAAGAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGa           <  1:3645338/71‑1 (MQ=255)
                                                       aaaTGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGc                 >  1:3755620/1‑64 (MQ=255)
                                                       aaaTGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCg            >  1:4223610/1‑69 (MQ=255)
                                                       aaaTGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCg            >  1:472628/1‑69 (MQ=255)
                                                       aaaTGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGaa          <  1:736714/71‑1 (MQ=255)
                                                        aaTGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAg         >  1:2853216/1‑71 (MQ=255)
                                                        aaTGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAg         >  1:3450735/1‑71 (MQ=255)
                                                         aTGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCGCGCAAACTCAGGCTGGCGAAGc        >  1:3000530/1‑71 (MQ=255)
                                                            ccGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCa                          <  1:2087841/50‑1 (MQ=255)
                                                            ccGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCa                    >  1:1754031/1‑56 (MQ=255)
                                                            ccGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCa                    >  1:216548/1‑56 (MQ=255)
                                                            ccGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCa                    >  1:3744219/1‑56 (MQ=255)
                                                            ccGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTg      <  1:560907/70‑1 (MQ=255)
                                                            ccGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTg      <  1:1592815/70‑1 (MQ=255)
                                                            ccGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTg      <  1:2536118/70‑1 (MQ=255)
                                                            ccGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTg      <  1:3749542/70‑1 (MQ=255)
                                                            ccGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTg      <  1:2941596/70‑1 (MQ=255)
                                                             cGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTGCg    <  1:134704/71‑1 (MQ=255)
                                                             cGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTGCg    <  1:781316/71‑1 (MQ=255)
                                                              gATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTGCGt   >  1:1997464/1‑71 (MQ=255)
                                                               aTAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTGCGtt  <  1:3043730/71‑1 (MQ=255)
                                                                   |                                                                  
TAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTACGTCGCGGAAATGCCGATAAAAAGAGGTGGGAGCAATGCCCGCTTCACGCGCCACTTCACGCAAACTCAGGCTGGCGAAGCTGCGTT  >  minE/2281054‑2281187

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: