Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 356,222 G→C 6.2% H199D (CAC→GAC)  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,2220GC 6.2% 300.9 / 7.0 97H199D (CAC→GAC) folDbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase and 5,10‑methylene‑tetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base G (46/45);  new base C (4/2);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 101 total aligned reads displayed.

ATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCC  >  minE/356155‑356286
                                                                   |                                                                
aTCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGtgat                                                                <  1:647548/70‑1 (MQ=255)
aTCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGtgat                                                                >  1:750706/1‑70 (MQ=255)
aTCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGtgat                                                                >  1:1230065/1‑70 (MQ=255)
 tCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGtgatga                                                              >  1:2841245/1‑71 (MQ=255)
  ccAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGAc                                                             >  1:2418920/1‑71 (MQ=255)
  ccAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGAc                                                             >  1:2186257/1‑71 (MQ=255)
   cAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGAc                                                             <  1:2221831/70‑1 (MQ=255)
   cAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACg                                                            >  1:1161135/1‑71 (MQ=255)
   cAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACg                                                            <  1:3303227/71‑1 (MQ=255)
    aGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACg                                                            >  1:1397475/1‑70 (MQ=255)
      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:871102/71‑3 (MQ=38)
         ccGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaaa                                                      <  1:3097243/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:70644/70‑1 (MQ=255)
                  aaGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGa                                              <  1:2973903/70‑1 (MQ=255)
                   aGCCTGTCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAg                                            <  1:2967062/71‑1 (MQ=255)
                   aGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAg                                            <  1:1609524/71‑1 (MQ=255)
                   aGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAg                                            <  1:147513/71‑1 (MQ=255)
                      cTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGc                                           >  1:1214753/1‑69 (MQ=255)
                         gCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                    >  1:2416302/1‑57 (MQ=255)
                         gCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGt                                        <  1:3223746/69‑1 (MQ=255)
                         gCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGt                                        >  1:1220474/1‑69 (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:1398935/1‑45 (MQ=21)
                                 aCGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaa                                                       >  1:3674418/1‑45 (MQ=21)
                                 aCGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaa                                                       >  1:140988/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:1197997/1‑70 (MQ=255)
                                 aCGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                               >  1:3836879/1‑70 (MQ=255)
                                 aCGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGt                              <  1:3110618/71‑1 (MQ=255)
                                  cGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                               >  1:771315/1‑69 (MQ=255)
                                  cGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                               >  1:530555/1‑69 (MQ=255)
                                  cGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                               >  1:229627/1‑69 (MQ=255)
                                   ggCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                            <  1:189688/71‑1 (MQ=255)
                                   ggCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                            <  1:501356/71‑1 (MQ=255)
                                     cAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                    <  1:2545156/45‑1 (MQ=255)
                                     cAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                    <  1:1317519/45‑1 (MQ=255)
                                     cAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAgtg                          <  1:1634308/71‑1 (MQ=255)
                                      aaCGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAgtgg                         >  1:3061742/1‑71 (MQ=255)
                                        cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                               >  1:2888654/1‑47 (MQ=255)
                                        cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                               >  1:24567/1‑47 (MQ=255)
                                        cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                               >  1:2291782/1‑47 (MQ=255)
                                        cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                               >  1:2097613/1‑47 (MQ=255)
                                        cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAgtggtg                       >  1:301287/1‑71 (MQ=255)
                                         gATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                    <  1:3586129/41‑1 (MQ=255)
                                         gATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                    <  1:2470064/41‑1 (MQ=255)
                                         gATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                    <  1:2508337/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:1421071/71‑1 (MQ=255)
                                                     ggCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagca          <  1:1300857/71‑1 (MQ=255)
                                                      gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcag         >  1:3922411/1‑71 (MQ=255)
                                                      gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcag         <  1:1488086/71‑1 (MQ=255)
                                                      gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcag         <  1:1267342/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)
                                                                   |                                                                
ATCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCC  >  minE/356155‑356286

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: