Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,216 C→T 100% Y190Y (TAC→TAT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,2160CT100.0% 194.0 / NA 65Y190Y (TAC→TATpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/37);  total (28/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGT  >  AP009048/3825150‑3825277
                                                                  |                                                             
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:4200201/1‑71 (MQ=255)
                tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                           >  1:3902058/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:44249/49‑1 (MQ=38)
                      tGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                           <  1:4575471/49‑1 (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)
                                                                  |                                                             
GGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGT  >  AP009048/3825150‑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: