Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,891,678 T→C 27.0% D111D (GAT→GAC ytfQ → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,891,6780TC27.0% 145.7 / 55.6 133D111D (GAT→GACytfQpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (49/48);  new base C (20/16);  total (69/64)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 134 total aligned reads displayed.

AAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAT  >  minE/2891636‑2891737
                                          |                                                           
aaaaGAGGCGAAAGATGCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAAtc                                   >  1:1011213/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAAtc                                   >  1:1064174/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAAtc                                   >  1:2276384/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAAtc                                   >  1:2809079/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAAtc                                   >  1:555131/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAAt                                    >  1:640047/1‑68 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTTCTCGATCGTTCCATTGATGTGAAAGACAAAtc                                   >  1:577194/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTc                                                        <  1:857584/48‑1 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTc                                                        <  1:632026/48‑1 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTc                                                        >  1:2833657/1‑48 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTc                                                        >  1:3071453/1‑48 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTc                                                        >  1:3287877/1‑48 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtc                                   >  1:1267536/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtc                                   >  1:2058382/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtc                                   >  1:3155210/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtc                                   >  1:394000/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtc                                   >  1:1176520/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtc                                   >  1:1088637/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtc                                   >  1:734304/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATCCCGCTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAAtc                                   >  1:2142988/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCCGAAATACCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtc                                   >  1:1249944/1‑69 (MQ=255)
aaaaGAGGCGAAAGATGCAGAAATCCCGGTATTCTTGCTCGAt                                                             >  1:3189473/1‑43 (MQ=255)
aaaaGAGGCGAAAGATGACGAAATCCCGGTATTCTTGCTCGATCGTTCCATAGATGTGAAAGACAAAt                                    >  1:2915436/1‑68 (MQ=255)
    gagGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtctctc                               <  1:2478787/69‑1 (MQ=255)
    gagGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtctctc                               <  1:577110/69‑1 (MQ=255)
    gagGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtctctc                               <  1:1884170/69‑1 (MQ=255)
    gagGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtctctc                               <  1:1823933/69‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:379015/46‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:2327048/46‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:2304811/46‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:2993751/46‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:1795354/46‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:543908/46‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:1898687/46‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:685507/46‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtg                                                    <  1:1967305/46‑2 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCta                             <  1:1486683/69‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCta                             <  1:488498/69‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCta                             <  1:807931/69‑1 (MQ=255)
      ggCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCta                             <  1:1348079/69‑1 (MQ=255)
      ggCGAAAGATGACGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:3197667/46‑1 (MQ=255)
       gCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCAtt                                                    <  1:1565902/45‑1 (MQ=255)
               tGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGaccac                    >  1:1434670/1‑69 (MQ=255)
                gCCGATATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:2289776/69‑1 (MQ=255)
                gCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:3322186/69‑1 (MQ=255)
                gCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:1796171/69‑1 (MQ=255)
                gCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:433441/69‑1 (MQ=255)
                gCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:884555/69‑1 (MQ=255)
                gCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:105863/69‑1 (MQ=255)
                gCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:2581896/69‑1 (MQ=255)
                gCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:1892750/69‑1 (MQ=255)
                gCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:194618/69‑1 (MQ=255)
                gCCGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:2209042/69‑1 (MQ=255)
                gCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAGAGACAAATCTCTCTATATGaccacc                   >  1:2144617/1‑69 (MQ=255)
                 ccGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:1007700/68‑1 (MQ=255)
                 ccGAAATCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGaccacc                   <  1:2178567/68‑1 (MQ=255)
                   gAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGt                 <  1:14747/68‑1 (MQ=255)
                   gAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:1596442/69‑1 (MQ=255)
                    aaaTCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGACCACCGtta               >  1:1024553/1‑67 (MQ=255)
                    aaaTCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGACCACCGtta               >  1:1063370/1‑67 (MQ=255)
                    aaaTCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGACCACCGtta               >  1:301449/1‑67 (MQ=255)
                    aaaTCCCGGTCTTCTTGCTTGACCGTTCCATCGATGTGAAAGACAAATCTCTCTATATGACCACCGtta               >  1:2193802/1‑67 (MQ=255)
                            gTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGt                 <  1:2229149/59‑1 (MQ=255)
                               ttCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGacaac    >  1:80789/1‑69 (MQ=255)
                               ttCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGacaac    >  1:147307/1‑69 (MQ=255)
                               ttCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGacaac    >  1:3125140/1‑69 (MQ=255)
                               ttCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGacaac    >  1:3126726/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTTCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:2434118/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACTAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:2499892/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAAtctctc                               >  1:2758008/1‑40 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:1221334/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:2662962/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:2365352/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:121427/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:3200434/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:501516/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:2118818/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:473807/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:1283920/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:365997/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:135363/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:422150/55‑1 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCa               >  1:2997748/1‑56 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGacaaca   >  1:2425988/1‑68 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGacaaca   >  1:1940302/1‑68 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGacaaca   >  1:2595487/1‑68 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:621220/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:3252916/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:811947/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:2828904/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:1997068/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:1088101/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:624381/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:593150/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:2008276/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:512736/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:3059372/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:414379/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAt  >  1:1630022/1‑69 (MQ=255)
                                 cTTGCTCGATCGTTCCATTGATATGAAAGACAAATCTCTCTATATGACCACCGTc                <  1:109484/55‑1 (MQ=255)
                                          |                                                           
AAAAGAGGCGAAAGATGCCGAAATCCCGGTATTCTTGCTCGATCGTTCCATTGATGTGAAAGACAAATCTCTCTATATGACCACCGTCACTGCCGACAACAT  >  minE/2891636‑2891737

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: