Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,880 C→G 100% A78A (GCC→GCG pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8800CG100.0% 239.3 / NA 71A78A (GCC→GCGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (35/36);  total (35/36)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCA  >  AP009048/3824812‑3824944
                                                                    |                                                                
gcattGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgac                                                                <  1:2632743/68‑1 (MQ=255)
    tggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccac                                                             >  1:4205767/1‑70 (MQ=255)
     ggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccacta                                                           >  1:4402566/1‑69 (MQ=255)
      gtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccacta                                                           <  1:806911/70‑3 (MQ=255)
          aTTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGcc                                                      <  1:3824718/71‑2 (MQ=255)
          aTTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGcc                                                      >  1:3821869/1‑70 (MQ=255)
           ttCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGccg                                                     >  1:3806270/1‑69 (MQ=255)
             ctggCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacca                                                              >  1:3458435/3‑60 (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)
                   ttGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgac                                                                <  1:5140721/52‑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:1964258/70‑1 (MQ=255)
                      aGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGc                                           <  1:4426331/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)
                              ctgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacca                                                              >  1:1285804/1‑43 (MQ=255)
                              ctgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacca                                                              >  1:696581/1‑43 (MQ=255)
                              ctgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                  >  1:4460260/1‑71 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                     <  1:4787212/65‑1 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                    <  1:2418029/66‑1 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                    <  1:1959734/66‑1 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                <  1:4081962/70‑3 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                               <  1:629433/71‑1 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                               >  1:1268068/1‑71 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                               <  1:4508300/71‑1 (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:939231/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)
                                              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:3138110/1‑71 (MQ=255)
                                                     aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt           >  1:2981782/1‑71 (MQ=255)
                                                     aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt           >  1:3842773/1‑71 (MQ=255)
                                                     aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt           >  1:2808094/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:3203776/42‑1 (MQ=16)
                                                               aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                              <  1:4284209/42‑1 (MQ=16)
                                                               aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                              <  1:4383776/42‑1 (MQ=16)
                                                                ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacc                             <  1:2145523/42‑1 (MQ=16)
                                                                 tGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa  >  1:3155027/1‑68 (MQ=255)
                                                                    |                                                                
GCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCA  >  AP009048/3824812‑3824944

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: