Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,948,347 C→A 16.2% intergenic (‑38/+81) gcvP ← / ← gcvH glycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex/glycine cleavage complex lipoylprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,948,3470CA16.2% 261.8 / 22.6 86intergenic (‑38/+81)gcvP/gcvHglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex/glycine cleavage complex lipoylprotein
Reads supporting (aligned to +/- strand):  ref base C (31/41);  new base A (9/5);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

TTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGATCCTCTCCCGCAGAAGAG  >  minE/1948281‑1948416
                                                                  |                                                                     
ttttCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACctctc                                                                   >  1:3841196/1‑71 (MQ=255)
 tttCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACctctct                                                                  >  1:3704325/1‑71 (MQ=255)
 tttCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCATGAAACGTGCAGTGAATTGTGAACatctc                                                                   <  1:2241299/70‑1 (MQ=255)
  ttCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACctctct                                                                  <  1:3546378/70‑1 (MQ=255)
  ttCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACctctct                                                                  >  1:2145737/1‑70 (MQ=255)
  ttCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCtctctc                                                                 >  1:3705117/1‑71 (MQ=255)
  ttCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCtctctc                                                                 >  1:169746/1‑71 (MQ=255)
  ttCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCtctctc                                                                 >  1:3198929/1‑71 (MQ=255)
      aGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAAcc                                                                       <  1:601852/61‑1 (MQ=255)
       gCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTAc                                                            >  1:2091168/1‑71 (MQ=255)
        cTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACg                                                           <  1:3149787/71‑1 (MQ=255)
        cTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTAc                                                            >  1:1442453/1‑70 (MQ=255)
        cTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTAc                                                            >  1:1489819/1‑70 (MQ=255)
        cTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACg                                                           <  1:3819649/71‑1 (MQ=255)
         tGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACctct                                                                    >  1:2929598/1‑61 (MQ=255)
         tGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACg                                                           >  1:3524428/1‑70 (MQ=255)
         tGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACg                                                           >  1:3682409/1‑70 (MQ=255)
         tGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACg                                                           >  1:3667129/1‑70 (MQ=255)
          ggCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACctct                                                                    >  1:3327687/1‑60 (MQ=255)
           gCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGa                                                          <  1:2228151/69‑1 (MQ=255)
             ttAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGaa                                                         >  1:357743/1‑68 (MQ=255)
               aaCGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACctctct                                                                  >  1:1833398/1‑57 (MQ=255)
               aaCGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAagagag                                                    >  1:1294802/1‑71 (MQ=255)
                aCGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAagagag                                                    <  1:2400177/70‑1 (MQ=255)
                aCGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAagagag                                                    <  1:2649869/70‑1 (MQ=255)
                aCGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgt                                                   <  1:90512/71‑1 (MQ=255)
                 cGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtg                                                  <  1:515016/71‑1 (MQ=255)
                 cGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgt                                                   >  1:3010805/1‑70 (MQ=255)
                 cGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgt                                                   >  1:907047/1‑70 (MQ=255)
                 cGTCTGTGTCATGAACGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgt                                                   >  1:988548/1‑70 (MQ=255)
                  gTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtga                                                 <  1:2119530/71‑1 (MQ=255)
                  gTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtga                                                 <  1:3454517/71‑1 (MQ=255)
                  gTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACGAAGAGAgt                                                   >  1:1336981/1‑69 (MQ=255)
                     tctgtCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgagg                                               <  1:2307603/68‑1 (MQ=255)
                      gtgtCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgag                                                <  1:2363651/68‑1 (MQ=255)
                       tgtCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgagggt                                             <  1:2384251/70‑1 (MQ=255)
                       tgtCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgagggt                                             <  1:290045/70‑1 (MQ=255)
                       tgtCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgaggg                                              <  1:1418421/69‑1 (MQ=255)
                       tgtCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACGAAGAGAgtgagggt                                             >  1:2105031/1‑70 (MQ=255)
                       tgtCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACGAAGAGAgtgagggt                                             <  1:3958688/70‑1 (MQ=255)
                       tgtCATGAGCGATGGATCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgagggt                                             <  1:1630481/70‑1 (MQ=255)
                        gtCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgagggt                                             >  1:398845/1‑69 (MQ=255)
                        gtCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgaggg                                              >  1:1461313/1‑68 (MQ=255)
                           aTGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTTAGGc                                        <  1:4093485/71‑1 (MQ=255)
                            tGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCt                                                               <  1:3207411/47‑1 (MQ=255)
                             gAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCa                                       >  1:1631350/1‑70 (MQ=255)
                              aGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTcc                                                                >  1:2612596/1‑44 (MQ=255)
                              aGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAagag                                                      >  1:1106258/1‑54 (MQ=255)
                              aGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgagggtgag                                          >  1:31944/1‑66 (MQ=255)
                              aGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgag                                                >  1:339373/1‑60 (MQ=255)
                              aGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCAt                                      >  1:1529256/1‑70 (MQ=255)
                              aGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCAt                                      >  1:18876/1‑70 (MQ=255)
                               gCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGaa                                                         <  1:2853123/50‑1 (MQ=255)
                               gCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATa                                     >  1:1679191/1‑70 (MQ=255)
                                cGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAagag                                                      <  1:3330272/52‑1 (MQ=255)
                                cGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGTTGAGGCATaaa                                   <  1:3107923/71‑1 (MQ=255)
                                  aTGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgagggt                                             <  1:2166550/59‑1 (MQ=255)
                                   tGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtgagggtga                                           >  1:3142303/1‑60 (MQ=255)
                                   tGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAAttt                                <  1:2486855/71‑1 (MQ=255)
                                    ggTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAAtttt                               >  1:1668607/1‑71 (MQ=255)
                                    ggTTCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAAtttt                               >  1:2279561/1‑71 (MQ=255)
                                     gTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGa                                                          <  1:1059765/43‑1 (MQ=255)
                                     gTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGc                                        >  1:3872905/1‑61 (MQ=255)
                                      ttCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAag                                                        <  1:2769486/44‑1 (MQ=255)
                                      ttCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTcc                             <  1:2374134/71‑1 (MQ=255)
                                      ttCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACGAaga                                                       <  1:2574576/45‑1 (MQ=255)
                                       tCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGaa                                                         >  1:1139357/1‑42 (MQ=255)
                                       tCCTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCt                            >  1:760001/1‑71 (MQ=255)
                                        ccTTAAACGTGCAGTGAATTGTTAACCTCTCTCCTTAAGaa                                                         <  1:3696254/41‑1 (MQ=37)
                                        ccTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCt                            <  1:423090/70‑1 (MQ=255)
                                         cTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTc                           <  1:2559798/70‑1 (MQ=255)
                                         cTGAAACGTGCAGTGAATTGTGAACATCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTc                           <  1:3054796/70‑1 (MQ=255)
                                          tGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAagag                                                      <  1:414163/42‑1 (MQ=255)
                                            aaaCGTGCAGTGAATTGTGAACCTCTCTCCTTACGaa                                                         <  1:1661919/37‑1 (MQ=255)
                                            aaaCGTGCAGTGAATTGTGAACCTCTCTCCTTACGaa                                                         <  1:2748026/37‑1 (MQ=255)
                                            aaaCGTGCAGTGAAGTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCAcc                        <  1:4036963/70‑1 (MQ=255)
                                             aaCGTGCAGTGAATTGTGAACCTCTCTCCTTACGaa                                                         <  1:1388214/36‑1 (MQ=255)
                                             aaCGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCAcc                        <  1:2798044/69‑1 (MQ=255)
                                              aCGTGCAGTGAATTGTGAACCTCTCTCCTTACGAaga                                                       >  1:3299910/1‑37 (MQ=255)
                                               cGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGa                    <  1:3746288/71‑1 (MQ=255)
                                                   cAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAgtga                                                 <  1:2207257/38‑1 (MQ=255)
                                                     gTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATaa                                    <  1:4053705/49‑1 (MQ=255)
                                                        aaTTGTGAACCTCTCTCCTTACGAAGAGAgtgaggg                                              >  1:388102/1‑36 (MQ=255)
                                                        aaTTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAAtttt                               >  1:2845152/1‑51 (MQ=255)
                                                        aaTTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGATCCTCTccc           >  1:1949101/1‑71 (MQ=255)
                                                          ttGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTTATCCTCTCCCg          >  1:3867417/1‑70 (MQ=255)
                                                            gtgAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCAcc                        <  1:3204371/54‑1 (MQ=255)
                                                            gtgAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGATCCTCTCCCg          >  1:2433118/1‑68 (MQ=255)
                                                             tgAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGATCCTCTCCCGCag       <  1:2192582/70‑1 (MQ=255)
                                                             tgAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGATCCTCTCCCGCag       <  1:123836/70‑1 (MQ=255)
                                                             tgAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGATCCTCTCACGCag       <  1:1624438/70‑1 (MQ=255)
                                                              gAACATCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCAcc                        >  1:2876679/1‑52 (MQ=255)
                                                               aaCCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCAcc                        <  1:2831643/51‑1 (MQ=255)
                                                               aaCCTCTCTCCTTACGAAGAGAGTGAGGGAGAGGCAAAAAtttt                               >  1:718174/1‑44 (MQ=38)
                                                                aCCTCTCTCCTTACGAAGAGAGTGAGGGGGAGGCATAAATTTTCCTCAcc                        <  1:392281/50‑1 (MQ=255)
                                                                 ccTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGATCCTCTCCCGCAGAagag  <  1:3844078/71‑1 (MQ=255)
                                                                 ccTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGATCCTCTCCCGCAGAagag  >  1:2152216/1‑71 (MQ=255)
                                                                  |                                                                     
TTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACGTGCAGTGAATTGTGAACCTCTCTCCTTACGAAGAGAGTGAGGGTGAGGCATAAATTTTCCTCACCCTGATCCTCTCCCGCAGAAGAG  >  minE/1948281‑1948416

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: