Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 356,229 C→T 6.4% L196L (CTG→CTA folD ← bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase and 5,10‑methylene‑tetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE356,2290CT 6.4% 293.0 / 6.8 93L196L (CTG→CTAfolDbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase and 5,10‑methylene‑tetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base C (43/44);  new base T (4/2);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGC  >  minE/356161‑356299
                                                                    |                                                                      
tCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAg                                                                      <  1:1766409/71‑1 (MQ=255)
 cACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAg                                                                      <  1:938132/70‑1 (MQ=255)
  aCCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAt                                                                    <  1:1173427/71‑1 (MQ=255)
  aCCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAt                                                                    >  1:1614411/1‑71 (MQ=255)
   ccGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaaa                                                                   <  1:3097243/71‑3 (MQ=38)
   ccGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaaa                                                                   <  1:871102/71‑3 (MQ=38)
     gggAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGc                                                                        <  1:3658331/64‑1 (MQ=255)
     gggAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 >  1:2089001/1‑71 (MQ=255)
     gggAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 <  1:1660627/71‑1 (MQ=255)
      ggAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAttttt                                                                <  1:1541796/71‑1 (MQ=255)
      ggAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 >  1:3431266/1‑70 (MQ=255)
      ggAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 >  1:1149837/1‑70 (MQ=255)
         aTAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAg                                                              <  1:371620/70‑1 (MQ=255)
         aTAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGt                                                             <  1:2377896/71‑1 (MQ=255)
            aaGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTa                                                               >  1:479060/1‑66 (MQ=255)
            aaGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                                            >  1:3521415/1‑69 (MQ=255)
            aaGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGa                                                           <  1:2973903/70‑1 (MQ=255)
            aaGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGa                                                           <  1:70644/70‑1 (MQ=255)
             aGCCTGTCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAg                                                         <  1:2967062/71‑1 (MQ=255)
             aGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAg                                                         <  1:147513/71‑1 (MQ=255)
             aGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAg                                                         <  1:1609524/71‑1 (MQ=255)
                cTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGc                                                        >  1:1214753/1‑69 (MQ=255)
                   gCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 >  1:2416302/1‑57 (MQ=255)
                   gCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGt                                                     >  1:1220474/1‑69 (MQ=255)
                   gCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGt                                                     <  1:3223746/69‑1 (MQ=255)
                   gCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGa                                                   <  1:2655944/71‑1 (MQ=255)
                    cTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaa                                                                    >  1:3724851/1‑52 (MQ=25)
                      tGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGt                                                 >  1:3341468/1‑70 (MQ=255)
                       gCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTc                                                >  1:379525/1‑70 (MQ=255)
                        ccAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGt                                                     <  1:904135/64‑1 (MQ=255)
                        ccAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCa                                               <  1:219120/70‑1 (MQ=255)
                         cAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACt                                             >  1:1783385/1‑71 (MQ=255)
                          aaCGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                                            <  1:2891818/71‑1 (MQ=255)
                           aCGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaa                                                                    >  1:3674418/1‑45 (MQ=21)
                           aCGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaa                                                                    >  1:140988/1‑45 (MQ=21)
                           aCGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaa                                                                    >  1:1398935/1‑45 (MQ=21)
                           aCGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTa                                                               >  1:333005/1‑51 (MQ=255)
                           aCGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTa                                                               >  1:1558075/1‑51 (MQ=255)
                           aCGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCa                                               >  1:1769264/1‑67 (MQ=255)
                           aCGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                                            >  1:3836879/1‑70 (MQ=255)
                           aCGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                                            >  1:1197997/1‑70 (MQ=255)
                           aCGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGt                                           <  1:3110618/71‑1 (MQ=255)
                            cGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                                            >  1:229627/1‑69 (MQ=255)
                            cGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                                            >  1:771315/1‑69 (MQ=255)
                            cGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                                            >  1:530555/1‑69 (MQ=255)
                             ggCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                                         <  1:501356/71‑1 (MQ=255)
                             ggCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                                         <  1:189688/71‑1 (MQ=255)
                               cAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 <  1:1317519/45‑1 (MQ=255)
                               cAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 <  1:2545156/45‑1 (MQ=255)
                               cAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAgtg                                       <  1:1634308/71‑1 (MQ=255)
                                aaCGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAgtgg                                      >  1:3061742/1‑71 (MQ=255)
                                  cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                                            >  1:2097613/1‑47 (MQ=255)
                                  cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                                            >  1:2888654/1‑47 (MQ=255)
                                  cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                                            >  1:2291782/1‑47 (MQ=255)
                                  cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                                            >  1:24567/1‑47 (MQ=255)
                                  cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                                            >  1:955595/1‑47 (MQ=255)
                                  cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAgtggtg                                    >  1:301287/1‑71 (MQ=255)
                                   gATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 <  1:3586129/41‑1 (MQ=255)
                                   gATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 <  1:2508337/41‑1 (MQ=255)
                                   gATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                                 <  1:2470064/41‑1 (MQ=255)
                                   gATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTa                                                               >  1:2601648/1‑43 (MQ=255)
                                   gATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTa                                                               >  1:1961045/1‑43 (MQ=255)
                                   gATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAg                                                         >  1:3136976/1‑49 (MQ=255)
                                   gATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                                         <  1:3511481/65‑1 (MQ=255)
                                   gATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGc                                   <  1:1283179/71‑1 (MQ=255)
                                     tCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTa                                                               >  1:24554/1‑41 (MQ=255)
                                     tCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCaa                                 >  1:891079/1‑71 (MQ=255)
                                       aaTAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAAcc                               <  1:3683357/71‑1 (MQ=255)
                                       aaTAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAAc                                <  1:2640658/70‑1 (MQ=255)
                                           gATCGGCATTTTCTACGTGATGACCCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCt                              >  1:3896844/1‑68 (MQ=255)
                                            aTCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCa                                  >  1:1118347/1‑63 (MQ=255)
                                            aTCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCa                                  >  1:1068672/1‑63 (MQ=255)
                                              cGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcag                         >  1:502660/1‑70 (MQ=255)
                                               ggCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagca                       <  1:2399872/71‑1 (MQ=255)
                                               ggCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagca                       <  1:1812996/71‑1 (MQ=255)
                                               ggCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagca                       <  1:271204/71‑1 (MQ=255)
                                               ggCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagca                       <  1:1300857/71‑1 (MQ=255)
                                               ggCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagca                       <  1:1421071/71‑1 (MQ=255)
                                                gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcag                      <  1:1267342/71‑1 (MQ=255)
                                                gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcag                      >  1:3922411/1‑71 (MQ=255)
                                                gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcag                      <  1:1488086/71‑1 (MQ=255)
                                                 cATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcagc                     <  1:1685698/71‑1 (MQ=255)
                                                  aTTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcagca                    <  1:1672339/71‑1 (MQ=255)
                                                  aTTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcag                      >  1:1810453/1‑69 (MQ=255)
                                                   ttttCTACGTGATGACGCAGATTTTTAGTGAAGCgg                                                      >  1:3148232/1‑36 (MQ=255)
                                                   ttttCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAAc                                <  1:1427141/58‑1 (MQ=255)
                                                     ttCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCa                                               >  1:1340208/1‑41 (MQ=255)
                                                      tCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAgtggt                                     >  1:1441213/1‑50 (MQ=255)
                                                      tCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGtt                 <  1:722935/70‑1 (MQ=255)
                                                        tACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAgtggt                                     >  1:445424/1‑48 (MQ=255)
                                                           gTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTcc               <  1:515493/67‑1 (MQ=255)
                                                              atgaCGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                                         <  1:1188666/38‑1 (MQ=255)
                                                              atgaCGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCa        <  1:3268128/71‑1 (MQ=255)
                                                               tgaCGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCAt       >  1:3137776/1‑71 (MQ=255)
                                                                    cAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGggc  <  1:3852654/71‑1 (MQ=255)
                                                                    cAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGggc  <  1:2938484/71‑1 (MQ=255)
                                                                    |                                                                      
TCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGC  >  minE/356161‑356299

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: