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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8140GA100.0% 259.6 / NA 86A56A (GCG→GCApyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (49/37);  total (49/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGC  >  AP009048/3824746‑3824879
                                                                    |                                                                 
ttcttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCAtt                                                                 >  1:4243934/1‑71 (MQ=255)
         gCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTc                                                        <  1:1859859/71‑1 (MQ=255)
         gCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTc                                                        >  1:178636/1‑71 (MQ=255)
          ccGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTct                                                       >  1:4940603/1‑70 (MQ=255)
              gCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGc                                                    <  1:3154869/70‑1 (MQ=255)
                  tttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                 <  1:3744749/69‑1 (MQ=255)
                    tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg                                              <  1:3847495/70‑1 (MQ=255)
                    tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg                                              >  1:5176121/1‑70 (MQ=255)
                    tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg                                              >  1:1688792/1‑70 (MQ=255)
                    tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                             >  1:311995/1‑71 (MQ=255)
                    tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                             >  1:3877081/1‑71 (MQ=255)
                    tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                             >  1:4968628/1‑71 (MQ=255)
                    tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                             >  1:882993/1‑71 (MQ=255)
                     aaTACCGGGCGCGACCTGGCACTGTTAGGCCGGTTTTACGCTGAAGCATTGGTGGATTc                                                        >  1:2212865/1‑59 (MQ=255)
                       tACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCg                                          <  1:3629610/71‑1 (MQ=255)
                       tACCGGGCGCGAACTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCg                                          <  1:3919408/71‑1 (MQ=255)
                          cGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTct                                                       <  1:4345755/55‑2 (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:2758583/66‑1 (MQ=255)
                                 gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                  <  1:4187777/66‑1 (MQ=255)
                                 gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                  <  1:3752101/66‑1 (MQ=255)
                                 gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGtt                                >  1:4015590/4‑71 (MQ=255)
                                 gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGtt                                >  1:3904078/4‑71 (MQ=255)
                                 gaccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGtt                                >  1:1591259/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:4710533/64‑1 (MQ=255)
                                         aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGttt                               <  1:1117619/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:1416716/1‑42 (MQ=38)
                                            gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAt                                                  >  1:1680935/1‑42 (MQ=38)
                                            gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAt                                                  >  1:4730681/1‑42 (MQ=38)
                                            gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                        >  1:759893/1‑52 (MQ=255)
                                            gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                      >  1:3377891/1‑70 (MQ=255)
                                             ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                        >  1:3610908/1‑51 (MQ=255)
                                             ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                        >  1:1464078/1‑51 (MQ=255)
                                             ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCtt                     <  1:1965041/70‑1 (MQ=255)
                                             ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                    >  1:2082471/1‑71 (MQ=255)
                                             ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                    >  1:4354797/1‑71 (MQ=255)
                                             ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                    >  1:3439326/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:3090463/50‑1 (MQ=255)
                                                    gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                  <  1:535622/50‑1 (MQ=255)
                                                    gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                  <  1:3373010/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:1798391/1‑70 (MQ=255)
                                                                aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc  >  1:2373371/1‑70 (MQ=255)
                                                                aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc  <  1:3447761/70‑1 (MQ=255)
                                                                    |                                                                 
TTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGC  >  AP009048/3824746‑3824879

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: