Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,826 C→T 100% S60S (TCC→TCT pyrE → orotate phosphoribosyltransferase

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

GCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTG  >  AP009048/3824760‑3824893
                                                                  |                                                                   
gCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGc                                                                  <  1:3154869/70‑1 (MQ=255)
    tttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                               <  1:3744749/69‑1 (MQ=255)
      tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg                                                            >  1:5176121/1‑70 (MQ=255)
      tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg                                                            <  1:3847495/70‑1 (MQ=255)
      tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg                                                            >  1:1688792/1‑70 (MQ=255)
      tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                                           >  1:311995/1‑71 (MQ=255)
      tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                                           >  1:882993/1‑71 (MQ=255)
      tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                                           >  1:4968628/1‑71 (MQ=255)
      tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                                           >  1:3877081/1‑71 (MQ=255)
         tACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCg                                                        <  1:3629610/71‑1 (MQ=255)
         tACCGGGCGCGAACTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCg                                                        <  1:3919408/71‑1 (MQ=255)
             gggCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg                                                            >  1:1735975/1‑63 (MQ=255)
               gcgcGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgc                                                  >  1:383580/1‑71 (MQ=255)
               gcgcGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctg                                                   >  1:1015085/1‑70 (MQ=255)
                cgcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgct                                                 >  1:2785846/1‑71 (MQ=255)
                 gcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTCGATTCTGGCATTGAGTTCGATctgctg                                                <  1:4800625/71‑1 (MQ=255)
                   gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                                <  1:4187777/66‑1 (MQ=255)
                   gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                                <  1:3752101/66‑1 (MQ=255)
                   gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                                <  1:2758583/66‑1 (MQ=255)
                   gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGtt                                              >  1:3904078/4‑71 (MQ=255)
                   gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGtt                                              >  1:1591259/4‑71 (MQ=255)
                   gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGtt                                              >  1:4015590/4‑71 (MQ=255)
                     ccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                                >  1:4956668/2‑67 (MQ=255)
                      cTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTgg                                           >  1:1269617/1‑71 (MQ=255)
                         gCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTgg                                           <  1:984147/68‑1 (MQ=255)
                         gCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGccc                                        >  1:2423547/1‑71 (MQ=255)
                         gCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGc                                          <  1:4571858/69‑1 (MQ=255)
                          cACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                               <  1:4712211/47‑1 (MQ=39)
                           aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGttt                                             <  1:1117619/64‑1 (MQ=255)
                           aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGttt                                             <  1:4710533/64‑1 (MQ=255)
                             tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTg                                                              <  1:2106426/45‑1 (MQ=38)
                             tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTgg                                           >  1:4642167/1‑64 (MQ=255)
                             tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTgg                                           <  1:3607521/64‑1 (MQ=255)
                             tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                                    >  1:3533853/1‑71 (MQ=255)
                             tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                                    >  1:3533850/1‑71 (MQ=255)
                              gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                               <  1:1791598/43‑1 (MQ=38)
                              gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAt                                                                >  1:909545/1‑42 (MQ=38)
                              gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAt                                                                >  1:4730681/1‑42 (MQ=38)
                              gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAt                                                                >  1:1416716/1‑42 (MQ=38)
                              gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAt                                                                >  1:1680935/1‑42 (MQ=38)
                              gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                      >  1:759893/1‑52 (MQ=255)
                              gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                                    >  1:3377891/1‑70 (MQ=255)
                               ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                      >  1:1464078/1‑51 (MQ=255)
                               ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                      >  1:3610908/1‑51 (MQ=255)
                               ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCtt                                   <  1:1965041/70‑1 (MQ=255)
                               ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                                  >  1:4354797/1‑71 (MQ=255)
                               ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                                  >  1:3439326/1‑71 (MQ=255)
                               ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                                  >  1:2082471/1‑71 (MQ=255)
                                tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTg                                                              >  1:1875243/1‑42 (MQ=38)
                                tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTg                                                              >  1:3478397/1‑42 (MQ=38)
                                tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGTGTTCGATctgctg                                                <  1:4618510/56‑1 (MQ=255)
                                tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                      >  1:1091863/1‑50 (MQ=255)
                                tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTg                                            >  1:481313/1‑60 (MQ=255)
                                 aGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                                           >  1:3938672/1‑44 (MQ=38)
                                 aGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTAc                                 >  1:3894684/1‑70 (MQ=255)
                                  ggCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACa                                <  1:3828044/70‑1 (MQ=255)
                                   gCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGccc                                        >  1:910664/1‑61 (MQ=255)
                                   gCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACaaa                              <  1:1225803/71‑1 (MQ=255)
                                   gCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACaaa                              <  1:1941044/71‑1 (MQ=255)
                                   gCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACa                                >  1:2577642/1‑69 (MQ=255)
                                   gCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACa                                <  1:3698390/69‑1 (MQ=255)
                                      gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                                <  1:3373010/50‑1 (MQ=255)
                                      gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                                <  1:3090463/50‑1 (MQ=255)
                                      gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                                <  1:535622/50‑1 (MQ=255)
                                      gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAgg                            >  1:3651200/1‑70 (MQ=255)
                                       tttttACGCTGAAGCATTTGTGGATTCTGGCATTGAGTTCGATct                                                    <  1:2434028/45‑1 (MQ=38)
                                       tttttACGCTGAAGCATTGGTGGATTCTGGCATTg                                                              <  1:1276960/35‑1 (MQ=37)
                                       tttttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGa                          <  1:1565237/71‑1 (MQ=255)
                                         tttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGa                                                       >  1:3293349/1‑40 (MQ=38)
                                         tttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATc                        >  1:946229/1‑71 (MQ=255)
                                         tttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATc                        >  1:3860668/1‑71 (MQ=255)
                                           tACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATccc                      >  1:3646043/1‑71 (MQ=255)
                                             cGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGa                    <  1:3783638/71‑1 (MQ=255)
                                              gCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGa                    <  1:2002136/70‑1 (MQ=255)
                                               cTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGAtt                  <  1:2337360/71‑1 (MQ=255)
                                               cTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGAtt                  <  1:3348861/71‑1 (MQ=255)
                                                  aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTg                 <  1:2020680/69‑1 (MQ=255)
                                                  aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc                >  1:2373371/1‑70 (MQ=255)
                                                  aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc                <  1:3447761/70‑1 (MQ=255)
                                                  aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc                >  1:1798391/1‑70 (MQ=255)
                                                    gcattGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgac             <  1:2632743/68‑1 (MQ=255)
                                                      attGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGa                          <  1:1320094/55‑1 (MQ=255)
                                                        tggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccac          >  1:4205767/1‑70 (MQ=255)
                                                         ggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccacta        >  1:4402566/1‑69 (MQ=255)
                                                          gtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccacta        <  1:806911/70‑3 (MQ=255)
                                                          ctggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCg                     <  1:650920/56‑1 (MQ=255)
                                                             gATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACa                                <  1:3835061/43‑1 (MQ=255)
                                                              aTTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGcc   >  1:3821869/1‑70 (MQ=255)
                                                              aTTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGcc   <  1:3824718/71‑2 (MQ=255)
                                                               ttCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGccg  >  1:3806270/1‑69 (MQ=255)
                                                                  |                                                                   
GCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTG  >  AP009048/3824760‑3824893

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