Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,123 C→G 100% R159R (CGC→CGG pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,1230CG100.0% 197.7 / NA 56R159R (CGC→CGGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (19/37);  total (19/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

GGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATC  >  AP009048/3825057‑3825189
                                                                  |                                                                  
ggAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc                                                                >  1:3189050/1‑71 (MQ=255)
 gagaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTTATCGTCAGGAACGGGGgcg                                                               <  1:32197/71‑1 (MQ=255)
 gagaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc                                                                <  1:927143/70‑1 (MQ=255)
  agaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcgc                                                              >  1:1673456/1‑71 (MQ=255)
  agaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcgc                                                              <  1:3545027/71‑1 (MQ=255)
   gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACgggg                                                                  >  1:1208203/1‑66 (MQ=255)
     ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                                           <  1:2763703/71‑1 (MQ=255)
        ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCgag                                                        <  1:1555226/71‑1 (MQ=255)
           aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAtt                                                     <  1:1617157/71‑1 (MQ=255)
            ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                                    <  1:4522319/71‑1 (MQ=255)
            ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                                    <  1:3087504/71‑1 (MQ=255)
            ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                                    <  1:3627377/71‑1 (MQ=255)
            ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGGCAGGATCGGGGGCGCGGCGAGAt                                                      <  1:620104/69‑1 (MQ=255)
             gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                                   <  1:1221193/68‑1 (MQ=255)
             gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                                   <  1:2935298/68‑1 (MQ=255)
             gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                                   <  1:3418387/68‑1 (MQ=255)
                 aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                                                >  1:4119517/1‑70 (MQ=255)
                 aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                                                <  1:4265230/70‑1 (MQ=255)
                    gcgcGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                                                <  1:305839/67‑1 (MQ=255)
                    gcgcGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                                            <  1:72610/71‑1 (MQ=255)
                     cgcgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTc                                           <  1:1016784/71‑1 (MQ=255)
                      gcgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                                            >  1:1325275/1‑69 (MQ=255)
                        gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg                                        <  1:1107995/71‑1 (MQ=255)
                         aCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                                            <  1:207328/66‑1 (MQ=255)
                             ttGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGa                                              <  1:3587500/60‑1 (MQ=255)
                                cTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCa                                          <  1:894915/61‑1 (MQ=255)
                                  ggCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGc                               >  1:4297877/1‑70 (MQ=255)
                                   gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCg                              <  1:2333377/70‑1 (MQ=255)
                                      tattGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCgaga                                                       <  1:2541920/40‑1 (MQ=37)
                                      tattGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTg                                  <  1:2278894/61‑1 (MQ=255)
                                       attGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGAt                          >  1:3945307/2‑70 (MQ=255)
                                        ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCgaga                                                       >  1:3202154/1‑40 (MQ=38)
                                        ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCgg                                                 >  1:384887/1‑46 (MQ=255)
                                        ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTa                        >  1:4006316/1‑71 (MQ=255)
                                         tGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAg                                     >  1:820145/1‑57 (MQ=255)
                                          gATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTAc                       >  1:2426494/1‑70 (MQ=255)
                                           aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCgg                                                 <  1:2537616/43‑1 (MQ=38)
                                           aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACaa                     <  1:748236/71‑1 (MQ=255)
                                           aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACaa                     <  1:3483607/71‑1 (MQ=255)
                                            tttCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                                                >  1:2375551/1‑43 (MQ=38)
                                             ttCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                                            <  1:1624203/46‑1 (MQ=255)
                                                    gATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCaaa              <  1:1799145/69‑1 (MQ=255)
                                                     aTCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTg           <  1:3599589/71‑1 (MQ=255)
                                                      tCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCa                <  1:3544199/65‑1 (MQ=255)
                                                        gTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACaa                     >  1:4069892/1‑58 (MQ=255)
                                                          cAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCa                <  1:784673/61‑1 (MQ=255)
                                                          cAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGAtctc      >  1:1534330/1‑71 (MQ=255)
                                                           aGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATctct     <  1:3954948/71‑1 (MQ=255)
                                                            ggAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTa    <  1:4268683/71‑1 (MQ=255)
                                                             gAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTa    <  1:2984781/70‑1 (MQ=255)
                                                             gAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTa    >  1:4384461/1‑70 (MQ=255)
                                                              aaCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGAtctc      >  1:2452095/1‑67 (MQ=255)
                                                              aaCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTatc  <  1:1195985/71‑1 (MQ=255)
                                                              aaCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTatc  <  1:223702/71‑1 (MQ=255)
                                                               aCGGGGGCGCGGGGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTat   >  1:2843365/1‑69 (MQ=255)
                                                               aCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCg                              >  1:70704/1‑42 (MQ=255)
                                                               aCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACaa                     >  1:2580791/1‑51 (MQ=255)
                                                                  |                                                                  
GGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATC  >  AP009048/3825057‑3825189

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: