Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,309,908 A→G 39.6% L28L (CTT→CTC 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,9080AG39.6% 41.1 / 66.8 96L28L (CTT→CTCyefMantitoxin of the YoeB‑YefM toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (33/25);  new base G (15/23);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

TCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGCTGT  >  minE/1309843‑1309975
                                                                 |                                                                   
tctTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                                                  >  1:479931/1‑69 (MQ=255)
           gTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                                                       <  1:831147/69‑1 (MQ=255)
           gTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                                                       <  1:3076446/69‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 <  1:2020499/69‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 <  1:1258091/69‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 <  1:513526/69‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 <  1:2778940/69‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 <  1:1390640/69‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 <  1:248848/69‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 <  1:358893/69‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:1708820/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:3033897/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:2856189/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:1791315/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:2486720/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:1975103/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:2481596/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:1993968/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:2360386/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:2117675/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:2896075/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:1639345/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:1622221/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:3136299/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:431734/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:1383413/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:491259/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:618697/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:1222331/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:1132205/54‑1 (MQ=255)
                 ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:961112/54‑1 (MQ=255)
                  tcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 <  1:2004000/68‑1 (MQ=255)
                  tcttcGAGTGACATCAGAACACAAGACTCTCCATTCTGACGAGTAATGAGGAt                                                                <  1:2520402/53‑1 (MQ=255)
                   cttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                                                  >  1:2453772/1‑50 (MQ=255)
                   cttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                                                  >  1:582899/1‑50 (MQ=255)
                   cttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAgg                                                                  >  1:601417/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:2830946/68‑1 (MQ=255)
                              aTCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat                                     <  1:2116336/68‑1 (MQ=255)
                              aTCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat                                     <  1:300450/68‑1 (MQ=255)
                              aTCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat                                     <  1:3043489/68‑1 (MQ=255)
                              aTCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat                                     <  1:2787354/68‑1 (MQ=255)
                                cAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                                  >  1:51856/1‑69 (MQ=255)
                                cAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                                  >  1:480263/1‑69 (MQ=255)
                                cAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                                  >  1:2311045/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:453807/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:2161010/1‑69 (MQ=255)
                                cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                                  >  1:244006/1‑69 (MQ=255)
                                cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                                  >  1:2388687/1‑69 (MQ=255)
                                cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                                  >  1:2472689/1‑69 (MQ=255)
                                cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                                  >  1:1960799/1‑69 (MQ=255)
                                cAGAACACAAGCCTCTCCATTCTGACGAGTAATGAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                                  >  1:1024826/1‑69 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCTGGGCATGATCTTCAACGGCTTTCATCAtt                                 >  1:2380576/1‑69 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 >  1:3303893/1‑53 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 >  1:3147460/1‑53 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 >  1:2492800/1‑53 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 >  1:129262/1‑53 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 >  1:2407530/1‑53 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 >  1:719827/1‑53 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 >  1:2695565/1‑53 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 >  1:2713948/1‑53 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                                 >  1:955678/1‑53 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt                                 >  1:1792158/1‑69 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt                                 >  1:140158/1‑69 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt                                 >  1:2390317/1‑69 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt                                 >  1:1233976/1‑69 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt                                 >  1:1224506/1‑69 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt                                 >  1:388725/1‑69 (MQ=255)
                                 aGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAtt                                 >  1:2620671/1‑69 (MQ=255)
                                     cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                                                       <  1:1513009/43‑1 (MQ=255)
                                     cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                                                       <  1:893917/43‑1 (MQ=255)
                                     cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                                                       <  1:643581/43‑1 (MQ=255)
                                     cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                                                       <  1:2944370/43‑1 (MQ=255)
                                     cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATg                                                       <  1:850317/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:3360290/1‑69 (MQ=255)
                                       caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGcc                           >  1:212194/1‑69 (MQ=255)
                                       caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGcc                           >  1:2567617/1‑69 (MQ=255)
                                       caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGc                            >  1:2662316/1‑68 (MQ=255)
                                                               ataaGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGc     >  1:1987313/1‑67 (MQ=255)
                                                               ataaGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGc     >  1:2939327/1‑67 (MQ=255)
                                                               ataaGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGc     >  1:2431368/1‑67 (MQ=255)
                                                               ataaGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAAATCTGACGCGCTTc       >  1:2729914/1‑65 (MQ=255)
                                                                 aaGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGCTGt  <  1:301622/68‑1 (MQ=255)
                                                                 aaGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCAGACGCGCTTCGCTGt  <  1:2726915/68‑1 (MQ=255)
                                                                 |                                                                   
TCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTTCGCTGT  >  minE/1309843‑1309975

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: