Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,213 T→C 100% A189A (GCT→GCC pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,2130TC100.0% 233.6 / NA 67A189A (GCT→GCCpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/37);  total (30/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGT  >  AP009048/3825148‑3825277
                                                                 |                                                                
cAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTatc                                                              >  1:3468556/1‑68 (MQ=255)
  ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                           <  1:1694183/71‑1 (MQ=255)
  ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTg                                                            >  1:4322960/1‑70 (MQ=255)
     aGTTGAGCGTGTTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAag                                                        <  1:1135267/71‑1 (MQ=255)
       ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                      <  1:2282409/71‑1 (MQ=255)
        tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                      >  1:2661754/1‑70 (MQ=255)
        tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAag                                                        >  1:2354713/1‑68 (MQ=255)
         gAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                      >  1:4377842/1‑69 (MQ=255)
         gAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                      >  1:2241810/1‑69 (MQ=255)
            cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                 <  1:67620/71‑1 (MQ=255)
            cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                 <  1:1646192/71‑1 (MQ=255)
             gTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAgaga                                                     >  1:4258595/1‑66 (MQ=255)
                 ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                            <  1:2615440/71‑1 (MQ=255)
                 ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                             <  1:3539447/70‑1 (MQ=255)
                  tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAgaga                                                     >  1:2594038/1‑61 (MQ=255)
                  tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                           >  1:3902058/1‑71 (MQ=255)
                  tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                           >  1:4200201/1‑71 (MQ=255)
                   aCAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                              >  1:2874556/1‑67 (MQ=255)
                   aCAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                           <  1:1318402/70‑1 (MQ=255)
                    cAACTGCTAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                          <  1:3585312/70‑1 (MQ=255)
                    cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                         <  1:1022560/71‑1 (MQ=255)
                    cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                         <  1:1742235/71‑1 (MQ=255)
                    cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                         <  1:4744873/71‑1 (MQ=255)
                     aaCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                               <  1:1021976/64‑1 (MQ=255)
                      aCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                       <  1:1868467/71‑1 (MQ=255)
                       cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                       >  1:4104293/1‑70 (MQ=255)
                        tGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                           <  1:4575471/49‑1 (MQ=38)
                        tGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                           <  1:44249/49‑1 (MQ=38)
                         gCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTat                                                               >  1:4789264/1‑43 (MQ=38)
                          cataGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGa                                     <  1:756402/66‑1 (MQ=255)
                          cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGa                                     <  1:1369398/69‑1 (MQ=255)
                           aaaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                              <  1:3925388/59‑1 (MQ=38)
                           aaaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAAc                                   <  1:4320331/70‑1 (MQ=255)
                             aGTGATCTCTATCATCACCCTGAATGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATc                                >  1:5670/1‑71 (MQ=255)
                             aGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCTGGAAATGGCGGAACAt                                 <  1:2347012/70‑1 (MQ=255)
                               tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACAt                                 >  1:1097074/1‑68 (MQ=255)
                               tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTg                              >  1:1179415/1‑71 (MQ=255)
                                  tctcTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCCGAACAt                                 <  1:12102/65‑1 (MQ=255)
                                    tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcgg                          >  1:3988191/1‑70 (MQ=255)
                                     ctATCATCACCATGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACAt                                 >  1:1806508/1‑62 (MQ=255)
                                      tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCt                               <  1:1672203/63‑1 (MQ=255)
                                       atcatcACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCTGt                      >  1:3085364/1‑71 (MQ=255)
                                        tcatcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                     <  1:2607877/71‑1 (MQ=255)
                                          atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTa                    <  1:5146238/70‑1 (MQ=255)
                                          atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTa                    <  1:4129509/70‑1 (MQ=255)
                                           tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                         >  1:3511601/1‑64 (MQ=38)
                                             aCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                             >  1:2558676/1‑42 (MQ=25)
                                               ccTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                          <  1:2986861/43‑1 (MQ=25)
                                                cTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                       >  1:4448692/1‑45 (MQ=37)
                                                cTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTa             <  1:3513343/71‑1 (MQ=255)
                                                  gAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATc           >  1:3041606/1‑71 (MQ=255)
                                                   aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGTAAATGGCGGAAc                                   >  1:448469/1‑46 (MQ=25)
                                                   aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa                   >  1:96517/1‑62 (MQ=38)
                                                   aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcg          >  1:1456536/1‑71 (MQ=255)
                                                   aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcg          >  1:1414763/1‑71 (MQ=255)
                                                   aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATc           <  1:4844396/70‑1 (MQ=255)
                                                    aaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAg                  >  1:3734219/1‑62 (MQ=38)
                                                    aaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgc         <  1:2182828/71‑1 (MQ=255)
                                                    aaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgc         <  1:62046/71‑1 (MQ=255)
                                                     aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCt              >  1:4268967/1‑65 (MQ=39)
                                                     aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg        <  1:3339955/71‑1 (MQ=255)
                                                       atcTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg        >  1:3445111/3‑69 (MQ=39)
                                                           gATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                      <  1:3291003/35‑1 (MQ=25)
                                                            aTTGCCTATCTGGAAGAGAAGCCGGGAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGt  <  1:3955526/70‑1 (MQ=255)
                                                            aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGaa      <  1:5081628/66‑1 (MQ=255)
                                                            aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAagag   <  1:3410005/69‑1 (MQ=255)
                                                            aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAag     <  1:1421428/67‑1 (MQ=255)
                                                                 |                                                                
CAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGT  >  AP009048/3825148‑3825277

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: