Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,892 T→C 100% A82A (GCT→GCC pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8920TC100.0% 284.1 / NA 83A82A (GCT→GCCpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (39/44);  total (39/44)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACC  >  AP009048/3824828‑3824959
                                                                |                                                                   
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)
    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)
          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:939231/70‑2 (MQ=255)
                      tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                          <  1:2721762/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:2981782/1‑71 (MQ=255)
                                     aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                          >  1:3138110/1‑71 (MQ=255)
                                     aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                          >  1:2808094/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: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)
                                                         ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                        >  1:252124/3‑53 (MQ=37)
                                                         ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      <  1:501875/69‑1 (MQ=255)
                                                         ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      <  1:2502336/69‑1 (MQ=255)
                                                         ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      <  1:4517601/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)
                                                                |                                                                   
GCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACC  >  AP009048/3824828‑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: