Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 556,858 A→T 100% L36Q (CTG→CAG)  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
*AP009048556,8580AT100.0% 204.9 / NA 68L36Q (CTG→CAG) folDbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase and 5,10‑methylene‑tetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (32/36);  total (32/36)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCT  >  AP009048/556790‑556924
                                                                    |                                                                  
cacaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGt                                                                  >  1:958239/1‑71 (MQ=255)
cacaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGt                                                                  >  1:681774/1‑71 (MQ=255)
cacaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGt                                                                  >  1:4740563/1‑71 (MQ=255)
cacaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGt                                                                  >  1:3043291/1‑71 (MQ=255)
 acaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGt                                                                  <  1:3669600/70‑1 (MQ=255)
  caAGCCTTGCGTTTGCTTGCTACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTcc                                                                <  1:2714225/71‑1 (MQ=255)
  caAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTcc                                                                <  1:50818/71‑1 (MQ=255)
   aaGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCt                                                               <  1:42059/71‑1 (MQ=255)
    aGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTg                                                              <  1:5068395/71‑1 (MQ=255)
    aGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTg                                                              <  1:2699063/71‑1 (MQ=255)
    aGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTg                                                              <  1:1939917/71‑1 (MQ=255)
         tGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGcc                                                         >  1:3822953/1‑71 (MQ=255)
          gCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGccc                                                        >  1:853071/1‑71 (MQ=255)
          gCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGccc                                                        >  1:1196611/1‑71 (MQ=255)
           cGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCg                                                       <  1:5013821/71‑1 (MQ=255)
           cGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCg                                                       <  1:1531594/71‑1 (MQ=255)
            gTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTg                                                              <  1:4637245/63‑1 (MQ=255)
               tgcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGt                                                   >  1:1121559/1‑71 (MQ=255)
                   tgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCgg                                               >  1:3831435/1‑71 (MQ=255)
                    gcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCt                                                               >  1:253336/1‑54 (MQ=255)
                    gcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGc                                              >  1:657785/1‑71 (MQ=255)
                    gcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGc                                              >  1:4993914/1‑71 (MQ=255)
                       aCATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGc                                           <  1:4934644/71‑1 (MQ=255)
                        cATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTc                                                                 >  1:1313092/1‑48 (MQ=255)
                          tAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTcc                                                 <  1:1781501/62‑1 (MQ=255)
                          tAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTcc                                                 <  1:1787932/62‑1 (MQ=255)
                          tAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCaa                                         >  1:5318479/1‑70 (MQ=255)
                           aaaTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTc                                                                 >  1:2184137/1‑45 (MQ=255)
                              tttGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTg                                                              <  1:5352986/45‑1 (MQ=255)
                              tttGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTg                                                              <  1:4264431/45‑1 (MQ=255)
                              tttGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTg                                                              <  1:3937512/45‑1 (MQ=255)
                                   gATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCt                                >  1:2677559/1‑70 (MQ=255)
                                   gATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTg                               <  1:45456/71‑1 (MQ=255)
                                     tGCAGGGTTACTACCCACCAGCACAACGGCCTGTc                                                                 >  1:4177273/1‑35 (MQ=255)
                                     tGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGaa                             >  1:1279687/1‑71 (MQ=255)
                                      gCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGa                              >  1:1756784/1‑69 (MQ=255)
                                       cAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACt                           <  1:1370866/71‑1 (MQ=255)
                                       cAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACt                           >  1:4701032/1‑71 (MQ=255)
                                         gggTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATAc                                      >  1:4622762/1‑58 (MQ=255)
                                         gggTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttt                         <  1:683029/71‑1 (MQ=255)
                                         gggTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttt                         <  1:1928172/71‑1 (MQ=255)
                                         gggTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGAAGTCCGGCTGCAATACGCGCCTGAACttt                         <  1:2479941/71‑1 (MQ=255)
                                          ggTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACtttt                        >  1:183374/1‑71 (MQ=255)
                                           gTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttttt                       <  1:5025846/71‑1 (MQ=255)
                                            ttACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTg                      <  1:1405737/71‑1 (MQ=255)
                                             tactacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttttt                       >  1:786072/1‑69 (MQ=255)
                                              actacCCACCAGCACAACGGCTTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGa                     <  1:5057389/70‑1 (MQ=255)
                                              actacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGa                     <  1:1087054/70‑1 (MQ=255)
                                              actacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAg                    >  1:1500421/1‑71 (MQ=255)
                                              actacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAg                    >  1:3877394/1‑71 (MQ=255)
                                                tacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttttt                       <  1:811214/66‑1 (MQ=255)
                                                 acCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCaa                 >  1:4462198/1‑71 (MQ=255)
                                                   ccagcaGCACAACGGCCTGGCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttt                         <  1:1410535/61‑1 (MQ=255)
                                                    caccaGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGaa                             >  1:1675478/1‑56 (MQ=255)
                                                     accaGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTc             <  1:1538602/71‑1 (MQ=255)
                                                      ccaGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCa            <  1:4271352/71‑1 (MQ=255)
                                                      ccaGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCa            <  1:3179991/71‑1 (MQ=255)
                                                      ccaGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCa            <  1:2121554/71‑1 (MQ=255)
                                                            caACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAgcgc      <  1:942652/71‑1 (MQ=255)
                                                             aaCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAAc                <  1:2131540/60‑1 (MQ=255)
                                                             aaCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCag           <  1:3606813/65‑1 (MQ=255)
                                                             aaCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCa     >  1:2081082/1‑71 (MQ=255)
                                                             aaCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCa     >  1:5024390/1‑71 (MQ=255)
                                                             aaCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCa     >  1:1796818/1‑71 (MQ=255)
                                                             aaCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCa     >  1:4238634/1‑71 (MQ=255)
                                                              aCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAAt                                        >  1:2302820/1‑35 (MQ=255)
                                                               cGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCag           <  1:3425908/63‑1 (MQ=255)
                                                                ggCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACct  <  1:4363630/71‑1 (MQ=255)
                                                                    |                                                                  
CACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCT  >  AP009048/556790‑556924

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: