Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,898 A→G 100% A84A (GCA→GCG pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8980AG100.0% 287.3 / NA 84A84A (GCA→GCGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/42);  total (42/42)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACC  >  AP009048/3824830‑3824959
                                                                    |                                                             
aTTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCt                                                             <  1:3988880/71‑1 (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)
        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:3829400/71‑1 (MQ=255)
               ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                              <  1:629433/71‑1 (MQ=255)
               ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                              >  1:1268068/1‑71 (MQ=255)
                 gTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCat                                                >  1:1156158/1‑66 (MQ=255)
                    tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                          <  1:2721762/70‑2 (MQ=255)
                    tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                          <  1:939231/70‑2 (MQ=255)
                    tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                          <  1:2589201/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:1026169/1‑58 (MQ=255)
                                   aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                       >  1:4270914/1‑58 (MQ=255)
                                   aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                          >  1:2808094/1‑71 (MQ=255)
                                   aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                          >  1:2981782/1‑71 (MQ=255)
                                   aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                          >  1:3138110/1‑71 (MQ=255)
                                   aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                          >  1:3842773/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:4383776/42‑1 (MQ=16)
                                             aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                                             <  1:4284209/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: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)
                                                       ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                        >  1:252124/3‑53 (MQ=37)
                                                       ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      <  1:2502336/69‑1 (MQ=255)
                                                       ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      <  1:4517601/69‑1 (MQ=255)
                                                       ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      <  1:501875/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:140365/1‑71 (MQ=255)
                                                          ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc   >  1:1450498/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)
                                                                    |                                                             
ATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACC  >  AP009048/3824830‑3824959

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: