Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,015,787 T→C 35.8% R264R (CGT→CGC metC → cystathionine beta‑lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,015,7870TC35.8% 47.0 / 46.3 78R264R (CGT→CGCmetCcystathionine beta‑lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (34/16);  new base C (17/11);  total (51/27)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

AGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTT  >  minE/2015720‑2015853
                                                                   |                                                                  
aGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTc                                                                   <  1:1907134/69‑1 (MQ=255)
   tGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCaa                                                                 >  1:3196763/1‑68 (MQ=255)
   tGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCaa                                                                 >  1:941310/1‑68 (MQ=255)
   tGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAAc                                                                >  1:3410918/1‑69 (MQ=255)
   tGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAAc                                                                >  1:548854/1‑69 (MQ=255)
               aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa                                                    <  1:2061819/69‑1 (MQ=255)
               aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa                                                    <  1:2121986/69‑1 (MQ=255)
               aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa                                                    <  1:2456814/69‑1 (MQ=255)
               aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa                                                    <  1:47791/69‑1 (MQ=255)
               aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa                                                    <  1:3034561/69‑1 (MQ=255)
               aTACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCa                                                    <  1:2849034/69‑1 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:2601667/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:2522895/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:555232/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:2113035/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:44302/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:2982367/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:33794/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:3127082/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:637813/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:1323027/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:1483998/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:1669815/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatcat                                                           >  1:324815/1‑60 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatca                                                            >  1:3319813/1‑59 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAAcatca                                                            >  1:3063512/1‑59 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAgata                                                   >  1:2691615/1‑65 (MQ=255)
                 aCCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTAGCGCCAAcatca                                                            >  1:1604232/1‑59 (MQ=255)
                     ccTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGa                                              <  1:2458232/69‑1 (MQ=255)
                           taACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTgg                                        >  1:373631/1‑69 (MQ=255)
                           taACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTgg                                        >  1:2870823/1‑69 (MQ=255)
                           taACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTgg                                        >  1:2868453/1‑69 (MQ=255)
                           taACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTgg                                        >  1:2114939/1‑69 (MQ=255)
                                   cGTGGCCTGCGCACATTGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:17255/69‑1 (MQ=255)
                                    gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:2752883/69‑1 (MQ=255)
                                    gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:687848/69‑1 (MQ=255)
                                    gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:3277069/69‑1 (MQ=255)
                                    gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:452136/69‑1 (MQ=255)
                                    gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:1800645/69‑1 (MQ=255)
                                    gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:920105/69‑1 (MQ=255)
                                    gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:592044/69‑1 (MQ=255)
                                    gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:1738206/69‑1 (MQ=255)
                                    gTGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:3203269/69‑1 (MQ=255)
                                     tGGCCTGCGCACATTAGGTGTGCGTTTGCGCCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggc                               <  1:2578563/68‑1 (MQ=255)
                                              cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa                     >  1:3092125/1‑69 (MQ=255)
                                              cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa                     >  1:995444/1‑69 (MQ=255)
                                              cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa                     >  1:2790791/1‑69 (MQ=255)
                                              cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa                     >  1:2407151/1‑69 (MQ=255)
                                              cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa                     >  1:239642/1‑69 (MQ=255)
                                              cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa                     >  1:219115/1‑69 (MQ=255)
                                              cacaTTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACa                     >  1:2089824/1‑69 (MQ=255)
                                                   tAGGTGTGCGTTTGGGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:1856463/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAAtggct                              >  1:1761841/1‑55 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACGTcc                  >  1:3348394/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:802292/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:949873/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:2850928/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:1043008/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:1490955/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:1825043/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:2024544/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:3172030/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:1383641/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:2687241/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:10284/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:1265260/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:2864416/1‑67 (MQ=255)
                                                   tAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATcc                  >  1:2892792/1‑67 (MQ=255)
                                                    aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa               <  1:2229207/69‑1 (MQ=255)
                                                    aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa               <  1:712547/69‑1 (MQ=255)
                                                    aGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCa               >  1:1704759/1‑69 (MQ=255)
                                                             tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg      <  1:2629259/69‑1 (MQ=255)
                                                             tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg      <  1:60189/69‑1 (MQ=255)
                                                             tttGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGcgcg      <  1:3002648/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)
                                                                   |                                                                  
AGATGGTCGATGCCGATACCGCCTATATAACCAGCCGTGGCCTGCGCACATTAGGTGTGCGTTTGCGTCAACATCATGAAAGCAGTCTGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTT  >  minE/2015720‑2015853

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: