Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,096,686 T→C 8.4% intergenic (‑310/‑79) ydiJ ← / → ydiK predicted FAD‑linked oxidoreductase/predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,096,6860TC 8.4% 237.2 / 11.5 83intergenic (‑310/‑79)ydiJ/ydiKpredicted FAD‑linked oxidoreductase/predicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (49/27);  new base C (4/3);  total (53/30)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

ACCACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGA  >  minE/1096618‑1096755
                                                                    |                                                                     
aCCACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTg                                                                      <  1:2633305/70‑1 (MQ=255)
 ccACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCt                                                                    >  1:3689681/1‑71 (MQ=255)
    caTTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATGGCCTGTTAAAAACCGCTGAAATTGCtcatc                                                                <  1:3643387/71‑1 (MQ=255)
     aTTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATCGCtcat                                                                 >  1:2182080/1‑70 (MQ=255)
      tttttttCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCtcatca                                                               >  1:2392273/1‑71 (MQ=255)
        tttttCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTcatcat                                                              >  1:1434131/1‑70 (MQ=255)
          tttCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTa                                                            <  1:2810135/70‑1 (MQ=255)
                acacGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATCGCTCATCATTa                                                            <  1:2369528/64‑1 (MQ=255)
                  acGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGa                                                   <  1:1534098/71‑1 (MQ=255)
                        aTCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTcg                                             >  1:3666133/1‑71 (MQ=255)
                        aTCAGAATTAAGGGCAGAATCGGACTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTcg                                             >  1:3006557/1‑71 (MQ=255)
                         tCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCtc                                                                   >  1:3956783/1‑48 (MQ=255)
                         tCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGt                                                     <  1:3204929/62‑1 (MQ=255)
                         tCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGtt                                                <  1:939697/67‑1 (MQ=255)
                          cAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCa                                                        <  1:3693088/58‑1 (MQ=255)
                          cAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTcgcg                                           >  1:2175540/1‑71 (MQ=255)
                           aGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTcatcat                                                              >  1:900432/1‑51 (MQ=255)
                           aGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAg                                                  >  1:1459204/1‑63 (MQ=255)
                           aGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAg                                                  >  1:170007/1‑63 (MQ=255)
                           aGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCgt                                          >  1:1364505/1‑71 (MQ=255)
                           aGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCgt                                          >  1:2289513/1‑71 (MQ=255)
                            gAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCgt                                          <  1:1783363/70‑1 (MQ=255)
                             aaTTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCgtgt                                        >  1:2940130/1‑71 (MQ=255)
                              aTTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGtt                                       <  1:1185480/71‑1 (MQ=255)
                              aTTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGtt                                       <  1:520627/71‑1 (MQ=255)
                              aTTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGtt                                       >  1:2150485/1‑71 (MQ=255)
                              aTTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGtt                                       >  1:2348814/1‑71 (MQ=255)
                               ttAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTa                                                            >  1:3672938/1‑49 (MQ=255)
                               ttAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGttt                                               >  1:2753352/1‑62 (MQ=255)
                               ttAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGttt                                               >  1:2835383/1‑62 (MQ=255)
                               ttAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGtt                                       >  1:3998820/1‑70 (MQ=255)
                               ttAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGtt                                       >  1:2325232/1‑70 (MQ=255)
                                 aaGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCtc                                                                   >  1:2288808/1‑40 (MQ=255)
                                 aaGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTa                                                            >  1:1701449/1‑47 (MQ=255)
                                 aaGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAc                                    >  1:3859333/1‑71 (MQ=255)
                                  aGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTc                                      <  1:2004112/68‑1 (MQ=255)
                                   gggCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAc                                    >  1:3470784/1‑69 (MQ=255)
                                   gggCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAc                                    >  1:833164/1‑69 (MQ=255)
                                   gggCAGAATCCGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTCAGTTTCGCGTGTTCAcg                                   >  1:3466147/1‑70 (MQ=255)
                                     gCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAcgtcg                                >  1:1429849/1‑71 (MQ=255)
                                        gAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCAtt                                                             <  1:1232887/39‑1 (MQ=255)
                                        gAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCAtt                                                             <  1:1413752/39‑1 (MQ=255)
                                        gAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCAtt                                                             <  1:2376641/39‑1 (MQ=255)
                                         aaTCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAcgtcgcgtc                            >  1:3296515/1‑71 (MQ=255)
                                          aTCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTAt                                                           >  1:1991636/1‑39 (MQ=255)
                                          aTCGGCCTGTTAAAAACCGCTGAAATTGATCATCATTa                                                            >  1:1764621/1‑38 (MQ=255)
                                           tCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGc                                                         <  1:736662/40‑1 (MQ=255)
                                           tCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTcga                          >  1:3156833/1‑71 (MQ=255)
                                             ggCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTcga                          <  1:2006853/69‑1 (MQ=255)
                                               ccTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGAt                      >  1:89271/1‑71 (MQ=255)
                                               ccTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGAt                      >  1:851596/1‑71 (MQ=255)
                                               ccTGTTAAAAACCGCTGAAATCGCTCATCATTATGCAGGTg                                                    <  1:2050086/41‑1 (MQ=255)
                                               ccTGTTAAAAACCGCTGAAATCGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGAt                      <  1:3675517/71‑1 (MQ=255)
                                                 tGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTc                                      <  1:3005358/53‑1 (MQ=255)
                                                   ttAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTcgcg                                           <  1:2564916/46‑1 (MQ=255)
                                                   ttAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAcgtcgcgtc                            <  1:371599/61‑1 (MQ=255)
                                                   ttAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAcgtcgcgtc                            <  1:895814/61‑1 (MQ=255)
                                                   ttAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAcgtcgcgtc                            <  1:876325/61‑1 (MQ=255)
                                                   ttAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAcgtcgcg                              <  1:2386680/59‑1 (MQ=255)
                                                   ttAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGAttt                    <  1:3463587/69‑1 (MQ=255)
                                                    tAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTcgacga                       <  1:224998/65‑1 (MQ=255)
                                                    tAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGAtt                     >  1:4120988/1‑67 (MQ=255)
                                                    tAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGa                  >  1:183363/1‑70 (MQ=255)
                                                     aaaaaCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGtt                                       >  1:1536185/1‑48 (MQ=255)
                                                     aaaaaCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAc                                    >  1:1682499/1‑51 (MQ=255)
                                                     aaaaaCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAc                                    >  1:1471996/1‑51 (MQ=255)
                                                     aaaaaCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGAc                 >  1:3419782/1‑70 (MQ=255)
                                                     aaaaaCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACg                >  1:3125460/1‑71 (MQ=255)
                                                     aaaaaCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACg                >  1:3548008/1‑71 (MQ=255)
                                                      aaaaCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGAtt                     >  1:2681639/1‑65 (MQ=255)
                                                       aaaCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGAttt                    <  1:518090/65‑1 (MQ=255)
                                                         aCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGAt                      >  1:869979/1‑61 (MQ=255)
                                                         aCCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGcac             >  1:3097940/1‑70 (MQ=255)
                                                          ccGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGcaca            >  1:3014475/1‑70 (MQ=255)
                                                          ccGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGcaca            >  1:2512984/1‑70 (MQ=255)
                                                          ccGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTAGCGTCGACGAtt                     >  1:1459336/1‑61 (MQ=255)
                                                          ccGCTGAAATCGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAc                                    >  1:1323414/1‑46 (MQ=255)
                                                          ccGCTGAAATCGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCAc                                    >  1:1752636/1‑46 (MQ=255)
                                                           agCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCa                                     >  1:3691560/2‑44 (MQ=255)
                                                            gCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACaaaa         <  1:1977755/71‑1 (MQ=255)
                                                             cTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACaaaaa        <  1:857243/71‑1 (MQ=255)
                                                             cTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACaaaa         >  1:2666516/1‑70 (MQ=255)
                                                               gAAATCGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAg      >  1:853150/1‑71 (MQ=255)
                                                                    tGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGa  >  1:3324919/1‑70 (MQ=255)
                                                                    |                                                                     
ACCACATTTTTTTCTAACACGCCCATCAGAATTAAGGGCAGAATCGGCCTGTTAAAAACCGCTGAAATTGCTCATCATTATGCAGGTGAGTTTCGCGTGTTCACGTCGCGTCGACGATTTGACGCACAAAAAAGGTGA  >  minE/1096618‑1096755

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: