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

GGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGTTT  >  AP009048/3825150‑3825279
                                                                  |                                                               
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:7330531/1‑69 (MQ=255)
              aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                                 >  1:5849133/1‑69 (MQ=255)
              aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                                 >  1:7036496/1‑69 (MQ=255)
              aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                               <  1:5783270/71‑1 (MQ=255)
              aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                               <  1:5910878/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:4790229/70‑1 (MQ=255)
                  cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                            <  1:6926011/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:7676805/67‑1 (MQ=255)
                        cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                         <  1:3258414/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:2667087/1‑51 (MQ=38)
                             tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                                    >  1:5288726/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:6962896/36‑1 (MQ=21)
                                    tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:2956585/36‑1 (MQ=21)
                                    tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:3878447/36‑1 (MQ=21)
                                    tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:5909923/36‑1 (MQ=21)
                                    tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:5738160/36‑1 (MQ=21)
                                    tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                            <  1:1290945/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:2507891/1‑42 (MQ=25)
                                     atcatcACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                                     >  1:5645982/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:1863989/71‑1 (MQ=255)
                                                   aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg          <  1:6358385/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:4204105/1‑51 (MQ=38)
                                                            tGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa                     >  1:1837806/1‑51 (MQ=38)
                                                                  |                                                               
GGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGTTT  >  AP009048/3825150‑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: