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% 317.2 / NA 93A189A (GCT→GCCpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/52);  total (41/52)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGTTT  >  AP009048/3825147‑3825279
                                                                  |                                                                  
tCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTatc                                                                >  1:6846132/1‑69 (MQ=255)
   ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                             <  1:1550208/71‑1 (MQ=255)
   ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTg                                                              >  1:7330236/1‑70 (MQ=255)
     aaGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                             >  1:6722452/1‑69 (MQ=255)
      aGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAag                                                          <  1:2094922/71‑1 (MQ=255)
       gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                             <  1:5848457/67‑1 (MQ=255)
       gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGaa                                                           >  1:802478/1‑69 (MQ=255)
       gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                         <  1:3594807/71‑1 (MQ=255)
       gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                         <  1:3269838/71‑1 (MQ=255)
       gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                         <  1:251676/71‑1 (MQ=255)
       gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                         <  1:145593/71‑1 (MQ=255)
        ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                        <  1:832222/71‑1 (MQ=255)
         tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                        >  1:7545473/1‑70 (MQ=255)
          gAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                         >  1:404399/1‑68 (MQ=255)
           aGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                                     >  1:2813475/1‑71 (MQ=255)
            gCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAgaga                                                       >  1:3147174/1‑68 (MQ=255)
            gCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                                    <  1:786453/71‑1 (MQ=255)
               tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                   <  1:7353509/69‑1 (MQ=255)
               tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                   <  1:5879748/69‑1 (MQ=255)
               tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                                 <  1:1606917/71‑1 (MQ=255)
                 aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                                 >  1:5849133/1‑69 (MQ=255)
                 aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                                 >  1:7036496/1‑69 (MQ=255)
                 aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                                 >  1:7330531/1‑69 (MQ=255)
                 aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                               <  1:5910878/71‑1 (MQ=255)
                 aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                               <  1:5783270/71‑1 (MQ=255)
                 aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                                <  1:5559905/70‑1 (MQ=255)
                   tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                             >  1:3717963/1‑71 (MQ=255)
                    aCAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                        >  1:1338930/1‑59 (MQ=38)
                     cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                           >  1:3776052/1‑71 (MQ=255)
                     cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                            <  1:6926011/70‑1 (MQ=255)
                     cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                            <  1:4790229/70‑1 (MQ=255)
                     cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                            <  1:987778/70‑1 (MQ=255)
                      aaCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                                    <  1:554199/61‑1 (MQ=255)
                      aaCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGc                                          <  1:1412897/71‑1 (MQ=255)
                        cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                         >  1:3510420/1‑70 (MQ=255)
                         tGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                        >  1:5284581/1‑70 (MQ=255)
                           cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                         <  1:3258414/67‑1 (MQ=255)
                           cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                         <  1:7676805/67‑1 (MQ=255)
                           cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGaa                                      >  1:3668861/1‑70 (MQ=255)
                           cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGa                                       <  1:2694892/69‑1 (MQ=255)
                           cAAAGTGATCTCTATCATCACCCTGAAAGATATGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                        >  1:6504479/1‑68 (MQ=255)
                            aaaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACa                                    <  1:2316765/71‑1 (MQ=255)
                               gTGATCTCTATCATCACCCTGAAAGATCTGATGGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCt                                 >  1:6465145/1‑71 (MQ=255)
                                tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                                    >  1:5288726/1‑51 (MQ=38)
                                tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                                    >  1:2667087/1‑51 (MQ=38)
                                tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                                    >  1:4897106/1‑51 (MQ=38)
                                  aTCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                             <  1:1772639/56‑1 (MQ=38)
                                     tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGaa                                                      <  1:5017695/44‑1 (MQ=25)
                                     tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                           >  1:3900297/1‑71 (MQ=255)
                                      ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGgcggcg                          >  1:2687971/1‑71 (MQ=255)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:5738160/36‑1 (MQ=21)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:1290945/36‑1 (MQ=21)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:5909923/36‑1 (MQ=21)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:3878447/36‑1 (MQ=21)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:2956585/36‑1 (MQ=21)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:6962896/36‑1 (MQ=21)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                           <  1:2485798/69‑1 (MQ=39)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                           <  1:2654715/69‑1 (MQ=39)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                           <  1:7011387/69‑1 (MQ=39)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGcggcgg                         >  1:6820695/1‑71 (MQ=255)
                                        atcatcACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:6046047/35‑1 (MQ=21)
                                        atcatcACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                                     >  1:5645982/1‑42 (MQ=25)
                                        atcatcACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                                     >  1:2507891/1‑42 (MQ=25)
                                         tcatcaCCCTGAAAGATCTGATTGCCTATCTGTAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                       <  1:2948899/71‑1 (MQ=255)
                                          catcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGc                                          >  1:3481267/1‑51 (MQ=38)
                                           atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                       >  1:1660342/1‑69 (MQ=39)
                                           atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa                     <  1:4819452/71‑1 (MQ=255)
                                            tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAg                    <  1:5303332/71‑1 (MQ=255)
                                            tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAg                    <  1:6820222/71‑1 (MQ=255)
                                             caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                             <  1:4128047/45‑1 (MQ=37)
                                             caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcgg                            <  1:3603088/62‑1 (MQ=38)
                                                ccTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAAc                                     <  1:2717707/50‑1 (MQ=38)
                                                 cTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAgg                   <  1:3597620/67‑1 (MQ=39)
                                                  tGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACAt                                   >  1:3660055/1‑50 (MQ=38)
                                                  tGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTAt              >  1:5741753/1‑71 (MQ=255)
                                                  tGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTAt              >  1:7101017/1‑71 (MQ=255)
                                                  agAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGCCGGTTaa                     >  1:5323970/2‑64 (MQ=255)
                                                   gAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                       <  1:371193/61‑1 (MQ=38)
                                                    aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAAc                                     >  1:417316/1‑46 (MQ=37)
                                                    aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcg            >  1:7479453/1‑71 (MQ=255)
                                                      aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg          <  1:6358385/71‑1 (MQ=255)
                                                      aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg          <  1:1863989/71‑1 (MQ=255)
                                                       gatcTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg          <  1:5757150/67‑1 (MQ=255)
                                                         tcTGATTGCCTATCTGGAAGAGAAGCCGGAAATTGCGGAACATCTggcg                             >  1:7514/2‑49 (MQ=37)
                                                         tcTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGCCCTATcgcg          <  1:1933089/67‑1 (MQ=38)
                                                            gATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTa                      <  1:2029221/53‑1 (MQ=38)
                                                             aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGcc                 >  1:2063292/1‑57 (MQ=38)
                                                             aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGcc                 >  1:4939391/1‑57 (MQ=38)
                                                             aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGcc                 >  1:5685929/1‑57 (MQ=38)
                                                             aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGtt   <  1:2957130/71‑1 (MQ=255)
                                                              ttGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGttt  <  1:6031585/71‑1 (MQ=255)
                                                               tGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa                     >  1:1837806/1‑51 (MQ=38)
                                                               tGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa                     >  1:4204105/1‑51 (MQ=38)
                                                                  |                                                                  
TCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGTTT  >  AP009048/3825147‑3825279

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: