Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,913 C→T 100% H89H (CAC→CAT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,9130CT100.0% 267.2 / NA 90H89H (CAC→CATpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (41/49);  total (41/49)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGG  >  AP009048/3824845‑3824980
                                                                    |                                                                   
ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                                                   <  1:629433/71‑1 (MQ=255)
ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                                                   <  1:4508300/71‑1 (MQ=255)
ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                                                   >  1:1268068/1‑71 (MQ=255)
ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                                                   <  1:3829400/71‑1 (MQ=255)
     tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                                               <  1:939231/70‑2 (MQ=255)
     tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                                               <  1:2589201/70‑2 (MQ=255)
     tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                                               <  1:2721762/70‑2 (MQ=255)
      ggCCCAGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                                              <  1:3362203/70‑3 (MQ=255)
        cccTGCTTACAAAGTGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacc                                                           >  1:1593928/1‑71 (MQ=255)
        cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacc                                                           >  1:3472502/1‑71 (MQ=255)
        cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                                             >  1:3717598/1‑69 (MQ=255)
         ccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgacctgacctg                                                         <  1:2464907/71‑1 (MQ=255)
          cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                         >  1:1797215/1‑71 (MQ=255)
          cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                                               <  1:1986340/65‑2 (MQ=255)
           tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                                        >  1:3930200/1‑71 (MQ=255)
             cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                      >  1:3303739/1‑71 (MQ=255)
              ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                      >  1:2299084/1‑70 (MQ=255)
                aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGt                                                     <  1:766191/69‑1 (MQ=255)
                aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                                                    <  1:1022044/70‑1 (MQ=255)
                  aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                                              <  1:3780966/58‑3 (MQ=255)
                  aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACt                                                  >  1:730349/1‑70 (MQ=255)
                  aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                                 >  1:5003074/1‑71 (MQ=255)
                    aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                                            >  1:4270914/1‑58 (MQ=255)
                    aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                                            >  1:1026169/1‑58 (MQ=255)
                    aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                               >  1:2808094/1‑71 (MQ=255)
                    aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                               >  1:3842773/1‑71 (MQ=255)
                    aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                               >  1:3138110/1‑71 (MQ=255)
                    aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                               >  1:2981782/1‑71 (MQ=255)
                      ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                             <  1:2266032/71‑1 (MQ=255)
                         tCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAc                                          <  1:4258885/71‑1 (MQ=255)
                          cccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTTCCGTACTGCTTTAAcc                                         >  1:5101386/1‑71 (MQ=255)
                           ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                         >  1:1268162/1‑54 (MQ=25)
                           ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCg                                        <  1:2970445/71‑1 (MQ=255)
                            cGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                                         >  1:616703/1‑69 (MQ=255)
                             gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                       >  1:1465747/1‑70 (MQ=255)
                             gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                                      <  1:3379486/71‑1 (MQ=255)
                             gATTGCGACCACTACCGCCGTGCCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                       >  1:4450107/1‑70 (MQ=255)
                              aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                                                                  <  1:4284209/42‑1 (MQ=16)
                              aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                                                                  <  1:4383776/42‑1 (MQ=16)
                              aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                                                                  <  1:3203776/42‑1 (MQ=16)
                               ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacc                                                                 <  1:2145523/42‑1 (MQ=16)
                                tGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                                      >  1:3155027/1‑68 (MQ=255)
                                 gcgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                                    >  1:2370290/4‑69 (MQ=255)
                                   gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                                <  1:3842048/70‑3 (MQ=255)
                                    accacTACCGCCGTGGCGCTGGCGGAGCATCCTGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                                >  1:4199562/1‑68 (MQ=255)
                                    accacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                                             >  1:2749798/1‑41 (MQ=16)
                                    accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                               <  1:2627195/71‑1 (MQ=255)
                                    accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                               <  1:221514/71‑1 (MQ=255)
                                     ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                               >  1:1566546/1‑70 (MQ=255)
                                     ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                               >  1:2139405/1‑70 (MQ=255)
                                     ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                               >  1:1800212/1‑70 (MQ=255)
                                      cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                             <  1:2069658/71‑1 (MQ=255)
                                        ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                             >  1:252124/3‑53 (MQ=37)
                                        ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                           <  1:4517601/69‑1 (MQ=255)
                                        ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                           <  1:501875/69‑1 (MQ=255)
                                        ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                           <  1:2502336/69‑1 (MQ=255)
                                        ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                             >  1:1627547/3‑69 (MQ=255)
                                         taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                             <  1:4923786/67‑1 (MQ=255)
                                         taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                             <  1:2607802/67‑1 (MQ=255)
                                           ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAATAGGCAAAAGa                         >  1:1178309/1‑70 (MQ=255)
                                           ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAg                          <  1:1888605/69‑1 (MQ=255)
                                           ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                        >  1:1450498/1‑71 (MQ=255)
                                           ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                        >  1:140365/1‑71 (MQ=255)
                                            cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                                                    <  1:1004759/42‑1 (MQ=25)
                                            cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAcc                       <  1:916455/71‑1 (MQ=255)
                                            cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                        >  1:3291487/1‑70 (MQ=255)
                                                   gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                 <  1:4913456/67‑1 (MQ=255)
                                                   gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGa                <  1:4902635/68‑1 (MQ=255)
                                                    cgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGaa               <  1:4943401/69‑1 (MQ=255)
                                                     gcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                             <  1:3198096/39‑1 (MQ=38)
                                                      cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                             <  1:3149869/55‑1 (MQ=255)
                                                      cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAgg             >  1:5121200/1‑71 (MQ=255)
                                                       tGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACgg                   >  1:1463398/1‑64 (MQ=255)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTa                                            <  1:2216593/37‑1 (MQ=38)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggggg          <  1:2219535/71‑1 (MQ=255)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg          <  1:3402987/71‑1 (MQ=255)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg          <  1:900494/71‑1 (MQ=255)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg          >  1:4742852/1‑71 (MQ=255)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg          <  1:1043984/71‑1 (MQ=255)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg           <  1:3770226/70‑1 (MQ=255)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggc            <  1:1582570/69‑1 (MQ=255)
                                                         gCGGAGCATCATGACCTAGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggc            <  1:2548262/69‑1 (MQ=255)
                                                           ggAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                             <  1:2627005/50‑1 (MQ=255)
                                                            gAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCaa       >  1:2876614/1‑71 (MQ=255)
                                                             aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGa                         >  1:4535494/1‑52 (MQ=255)
                                                             aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt      >  1:4546732/1‑71 (MQ=255)
                                                             aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt      <  1:446042/71‑1 (MQ=255)
                                                               catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCa        >  1:4457915/1‑67 (MQ=255)
                                                               catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCt    <  1:1691485/71‑1 (MQ=255)
                                                                 tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                 >  1:4838218/1‑56 (MQ=255)
                                                                 tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                 <  1:1041366/56‑1 (MQ=255)
                                                                 tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGa                <  1:180695/57‑1 (MQ=255)
                                                                 tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtgg  <  1:227643/71‑1 (MQ=255)
                                                                    |                                                                   
CTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGG  >  AP009048/3824845‑3824980

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: