Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,886 A→T 100% T80T (ACA→ACT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8860AT100.0% 235.7 / NA 80T80T (ACA→ACTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (41/39);  total (41/39)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAA  >  AP009048/3824822‑3824953
                                                                |                                                                   
aTTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGcc                                                               <  1:3824718/71‑2 (MQ=255)
aTTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGcc                                                               >  1:3821869/1‑70 (MQ=255)
 ttCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGccg                                                              >  1:3806270/1‑69 (MQ=255)
      gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTAc                                                                   >  1:3499562/1‑61 (MQ=255)
      gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCcgtg                                                            >  1:4585173/1‑68 (MQ=255)
      gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTgg                                                           <  1:1808581/69‑1 (MQ=255)
       cATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcgc                                                        <  1:1165925/71‑3 (MQ=255)
        aTTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCt                                                       <  1:3988880/71‑1 (MQ=255)
         ttGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTAcc                                                                  >  1:4975323/1‑59 (MQ=255)
          tGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGccg                                                              <  1:3550778/62‑3 (MQ=255)
          tGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcg                                                         <  1:4321181/67‑2 (MQ=255)
           gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTgg                                                     <  1:3220473/70‑1 (MQ=255)
           gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGc                                                    >  1:1200662/1‑71 (MQ=255)
            aGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGc                                                    <  1:4426331/70‑1 (MQ=255)
            aGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGc                                                    <  1:1964258/70‑1 (MQ=255)
                cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTAc                                                                   <  1:3679308/51‑1 (MQ=255)
                cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                >  1:99006/1‑70 (MQ=255)
                cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                >  1:1682498/1‑70 (MQ=255)
                  aTCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                             >  1:3231265/1‑71 (MQ=255)
                    ctgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                           >  1:4460260/1‑71 (MQ=255)
                       ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                              <  1:4787212/65‑1 (MQ=255)
                       ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                             <  1:1959734/66‑1 (MQ=255)
                       ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                             <  1:2418029/66‑1 (MQ=255)
                       ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                         <  1:4081962/70‑3 (MQ=255)
                       ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                        <  1:4508300/71‑1 (MQ=255)
                       ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                        <  1:629433/71‑1 (MQ=255)
                       ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                        >  1:1268068/1‑71 (MQ=255)
                       ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                        <  1:3829400/71‑1 (MQ=255)
                         gTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCat                                          >  1:1156158/1‑66 (MQ=255)
                            tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                    <  1:2589201/70‑2 (MQ=255)
                            tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                    <  1:2721762/70‑2 (MQ=255)
                            tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                    <  1:939231/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)
                                    cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTgg                                                     >  1:4819124/1‑45 (MQ=25)
                                    cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                           >  1:3724200/1‑55 (MQ=255)
                                    cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                           >  1:3303739/1‑71 (MQ=255)
                                     ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                           >  1:2299084/1‑70 (MQ=255)
                                       aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTgg                                                     <  1:4673419/42‑1 (MQ=25)
                                       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:3842773/1‑71 (MQ=255)
                                           aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                    >  1:2808094/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)
                                                tCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                             <  1:4650532/41‑1 (MQ=21)
                                                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:4383776/42‑1 (MQ=16)
                                                     aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                                       <  1:4284209/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:221514/71‑1 (MQ=255)
                                                           accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc    <  1:2627195/71‑1 (MQ=255)
                                                            ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc    >  1:2139405/1‑70 (MQ=255)
                                                            ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc    >  1:1800212/1‑70 (MQ=255)
                                                            ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc    >  1:1566546/1‑70 (MQ=255)
                                                             cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa  <  1:2069658/71‑1 (MQ=255)
                                                                |                                                                   
ATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAA  >  AP009048/3824822‑3824953

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