Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,015,844 T→A 24.6% V283V (GTT→GTA metC → cystathionine beta‑lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,015,8440TA24.6% 197.5 / 65.9 126V283V (GTT→GTAmetCcystathionine beta‑lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (17/78);  new base A (19/12);  total (36/90)
Fisher's exact test for biased strand distribution p-value = 1.43e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 126 total aligned reads displayed.

TTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTT  >  minE/2015781‑2015906
                                                               |                                                              
tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg                                                           <  1:60189/69‑1 (MQ=255)
tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg                                                           <  1:3002648/69‑1 (MQ=255)
tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg                                                           <  1:2629259/69‑1 (MQ=255)
  tGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGa                                                          >  1:2247121/1‑68 (MQ=255)
     gTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt                                                       <  1:500730/68‑1 (MQ=255)
     gTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt                                                       <  1:443798/68‑1 (MQ=255)
     gTCAACATCATGAAAGCAGTATGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGtt                                                       <  1:3399390/68‑1 (MQ=255)
       cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                    <  1:1630664/69‑1 (MQ=255)
       cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                    <  1:79589/69‑1 (MQ=255)
       cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                    <  1:1086066/69‑1 (MQ=255)
       cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                    <  1:1208876/69‑1 (MQ=255)
       cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                    <  1:2361095/69‑1 (MQ=255)
       cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                    <  1:3059581/69‑1 (MQ=255)
       cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                    <  1:136949/69‑1 (MQ=255)
       cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAac                                                    <  1:2509861/69‑1 (MQ=255)
       cAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCCAGTTAac                                                    <  1:1038301/69‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                                 <  1:2424879/69‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                                 <  1:3321296/69‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                                 <  1:708550/69‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                                 <  1:990212/69‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                                 <  1:1544209/69‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                                 <  1:1236079/69‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                                  <  1:359847/68‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                                  <  1:3208546/68‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                                  <  1:2923755/68‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                                  <  1:3429090/68‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                                  <  1:2551652/68‑1 (MQ=255)
          catcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                                  <  1:2176321/68‑1 (MQ=255)
           atcatGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAaccac                                                 <  1:2199582/68‑1 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTTGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:1294104/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                            <  1:2479816/68‑1 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                            <  1:385221/68‑1 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                            <  1:1414412/68‑1 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                            <  1:628064/68‑1 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                            <  1:3323792/68‑1 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                            <  1:3183957/68‑1 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGc                                            <  1:2971152/68‑1 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGtt                                                       >  1:2487544/1‑57 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTg                                             >  1:2303020/1‑67 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:1305925/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:283834/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:250423/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:2484602/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:1340347/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:3011618/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:1252073/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:452796/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:2315036/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:2203668/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:2083284/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:1655359/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:1692589/1‑68 (MQ=255)
                gAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGc                                            >  1:758152/1‑68 (MQ=255)
                         cTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAacca                                                  <  1:1854868/53‑1 (MQ=255)
                              aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                             <  1:1491898/69‑1 (MQ=255)
                              aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                             <  1:546135/69‑1 (MQ=255)
                              aGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                             <  1:2706342/69‑1 (MQ=255)
                               gTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTaa                             <  1:1060870/68‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:2928959/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:3069843/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:982790/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:322791/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:745524/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:637286/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:1162778/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:2751695/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:2266201/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:2037846/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:2911939/49‑1 (MQ=255)
                                 ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCt                                              <  1:2931938/49‑1 (MQ=255)
                                          gctggcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAAtt                 >  1:2927828/1‑69 (MQ=255)
                                          gctggcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAAtt                 >  1:3300975/1‑69 (MQ=255)
                                          gctggcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAAtt                 >  1:3131526/1‑69 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:2641345/67‑1 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:746635/67‑1 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:446231/67‑1 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:1587551/67‑1 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:561828/67‑1 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:2216875/67‑1 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:2038803/67‑1 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:1647734/67‑1 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGAGTTAAACACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:2519435/67‑1 (MQ=255)
                                          gctggcAGAACATCCCCAAGTAGCGCGACTTAACCACCCTGCTCTGCCTGGCAGTAAAGGCCACGaa                   <  1:2296054/67‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           <  1:963995/69‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           <  1:1867486/69‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           <  1:1889297/69‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           <  1:3040452/69‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           <  1:157391/69‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           <  1:1754557/69‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           <  1:306179/69‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGa            <  1:2303998/68‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGa            <  1:277958/68‑1 (MQ=255)
                                                aGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGa            <  1:1366272/68‑1 (MQ=255)
                                                   aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           >  1:444014/1‑66 (MQ=255)
                                                   aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           >  1:1117412/1‑66 (MQ=255)
                                                   aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGaa           >  1:2446180/1‑66 (MQ=255)
                                                          cAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACtt  >  1:201109/1‑68 (MQ=255)
                                                          cAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACtt  >  1:509817/1‑68 (MQ=255)
                                                          cAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACtt  >  1:3342078/1‑68 (MQ=255)
                                                          cAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACtt  >  1:102047/1‑68 (MQ=255)
                                                               |                                                              
TTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTT  >  minE/2015781‑2015906

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: