Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 356,225 G→C 6.2% H198D (CAT→GAT)  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,2250GC 6.2% 303.8 / 7.0 97H198D (CAT→GAT) 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 = 9.52e-01

Only 100 of 101 total aligned reads displayed.

TCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCAT  >  minE/356156‑356294
                                                                     |                                                                     
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:3097243/71‑3 (MQ=38)
        ccGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTGCGTCATCACGTAGaaa                                                              <  1:871102/71‑3 (MQ=38)
          gggAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGc                                                                   <  1:3658331/64‑1 (MQ=255)
          gggAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                            <  1:1660627/71‑1 (MQ=255)
          gggAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGAtttt                                                            >  1:2089001/1‑71 (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: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: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:530555/1‑69 (MQ=255)
                                 cGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                                       >  1:229627/1‑69 (MQ=255)
                                 cGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTg                                       >  1:771315/1‑69 (MQ=255)
                                  ggCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                                    <  1:189688/71‑1 (MQ=255)
                                  ggCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                                    <  1:501356/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:24567/1‑47 (MQ=255)
                                       cGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTg                                                       >  1:2291782/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: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)
                                                    ggCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagca                  <  1:2399872/71‑1 (MQ=255)
                                                     gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcag                 <  1:1488086/71‑1 (MQ=255)
                                                     gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCcagcag                 >  1:3922411/1‑71 (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)
                                                                   atgaCGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAg                                    <  1:1188666/38‑1 (MQ=255)
                                                                   atgaCGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCa   <  1:3268128/71‑1 (MQ=255)
                                                                    tgaCGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCAt  >  1:3137776/1‑71 (MQ=255)
                                                                     |                                                                     
TCCAGTCACCGGGAATAAAGCCTGGCTTGCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCAT  >  minE/356156‑356294

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: