Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,309,899 G→A 34.4% R31R (CGC→CGT yefM ← antitoxin of the YoeB‑YefM toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,309,8990GA34.4% 84.7 / 61.0 110R31R (CGC→CGTyefMantitoxin of the YoeB‑YefM toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base A (15/23);  total (54/56)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01

Only 100 of 110 total aligned reads displayed.

AAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCC  >  minE/1309835‑1309950
                                                                |                                                   
aaGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCg                                                   >  1:2386704/1‑67 (MQ=255)
aaGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGa                                                  >  1:348990/1‑68 (MQ=255)
aaGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGa                                                  >  1:695816/1‑68 (MQ=255)
aaGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAg                                                 <  1:3022147/69‑1 (MQ=255)
aaGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAg                                                 <  1:2599423/69‑1 (MQ=255)
aaGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAg                                                 <  1:2842412/69‑1 (MQ=255)
aaGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAg                                                 <  1:1315114/69‑1 (MQ=255)
aaGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAg                                                 <  1:1276272/69‑1 (MQ=255)
 aGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAg                                                 <  1:96555/68‑1 (MQ=255)
      tctcTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAg                                                 <  1:208149/63‑1 (MQ=255)
        tctTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGt                                                >  1:2749678/1‑62 (MQ=255)
        tctTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                         >  1:479931/1‑69 (MQ=255)
            ccAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCg                                                   <  1:1453520/55‑1 (MQ=255)
                   gTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:831147/69‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        <  1:2020499/69‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        <  1:513526/69‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        <  1:358893/69‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        <  1:2778940/69‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        <  1:1390640/69‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        <  1:1258091/69‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        <  1:248848/69‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:2481596/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:2117675/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:2360386/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:1993968/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:2486720/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:1975103/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:2856189/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:2896075/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:1791315/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:3033897/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:1708820/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:3136299/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:1639345/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:431734/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:491259/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:1383413/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:618697/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:961112/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:1222331/54‑1 (MQ=255)
                         ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:1132205/54‑1 (MQ=255)
                          tcttcGAGTGACATCAGAACACAAGACTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:2520402/53‑1 (MQ=255)
                           cttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                         >  1:2453772/1‑50 (MQ=255)
                           cttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                         >  1:601417/1‑50 (MQ=255)
                           cttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                         >  1:582899/1‑50 (MQ=255)
                            ttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                       <  1:3110183/51‑1 (MQ=255)
                                  tGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:3332261/54‑1 (MQ=255)
                                  tGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:1739227/54‑1 (MQ=255)
                                   gACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:2430894/53‑1 (MQ=255)
                                   gACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtc              >  1:1588838/1‑69 (MQ=255)
                                      aTCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat            <  1:300450/68‑1 (MQ=255)
                                      aTCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat            <  1:2830946/68‑1 (MQ=255)
                                      aTCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat            <  1:2116336/68‑1 (MQ=255)
                                      aTCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat            <  1:3043489/68‑1 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:2311045/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:480263/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:51856/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACTAGTAATGAGGa                                        >  1:892967/1‑38 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTtcatca          >  1:2607878/1‑68 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:1459608/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:140894/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:1700828/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:618862/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:701266/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:109188/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:244006/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:1024826/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:2161010/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:2472689/1‑69 (MQ=255)
                                        cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat         >  1:1960799/1‑69 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCTGGGCATGATCTTCAACGGCTTTCATCAtt        >  1:2380576/1‑69 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        >  1:955678/1‑53 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        >  1:3147460/1‑53 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        >  1:3303893/1‑53 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        >  1:2713948/1‑53 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        >  1:2695565/1‑53 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        >  1:719827/1‑53 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        >  1:2407530/1‑53 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        >  1:129262/1‑53 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                        >  1:2492800/1‑53 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt        >  1:1224506/1‑69 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt        >  1:1233976/1‑69 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt        >  1:2390317/1‑69 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt        >  1:140158/1‑69 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt        >  1:2620671/1‑69 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt        >  1:388725/1‑69 (MQ=255)
                                         aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt        >  1:1792158/1‑69 (MQ=255)
                                             cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:2944370/43‑1 (MQ=255)
                                             cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:643581/43‑1 (MQ=255)
                                             cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:1513009/43‑1 (MQ=255)
                                             cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:850317/43‑1 (MQ=255)
                                             cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:893917/43‑1 (MQ=255)
                                             cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttgttg    >  1:1350961/1‑69 (MQ=255)
                                             cacaAGCCTCTCCATTCTGGCGAGAAATAAGGATCGGGGCATg                              <  1:139670/43‑1 (MQ=255)
                                               caAGGCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttg       >  1:1854201/1‑64 (MQ=255)
                                               caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                              <  1:3046661/41‑1 (MQ=255)
                                               caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGcc  >  1:2567617/1‑69 (MQ=255)
                                               caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGcc  >  1:3360290/1‑69 (MQ=255)
                                               caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGcc  >  1:212194/1‑69 (MQ=255)
                                               caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGc   >  1:2662316/1‑68 (MQ=255)
                                                                |                                                   
AAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCC  >  minE/1309835‑1309950

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: