Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,108 C→T 100% L154L (CTC→CTT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,1080CT100.0% 199.9 / NA 67L154L (CTC→CTTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/46);  total (21/46)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAA  >  AP009048/3825040‑3825170
                                                                    |                                                              
gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                                                              >  1:3304369/1‑71 (MQ=255)
gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                                                              >  1:2788483/1‑71 (MQ=255)
 cGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                                                              >  1:4403973/1‑70 (MQ=255)
  gATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                            >  1:2647793/1‑71 (MQ=255)
    ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                          <  1:3742535/71‑1 (MQ=255)
    ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                          <  1:2687081/71‑1 (MQ=255)
    ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                          <  1:517736/71‑1 (MQ=255)
     tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                          >  1:4349746/1‑70 (MQ=255)
      cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                           <  1:1977238/68‑1 (MQ=255)
      cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTc                                                         <  1:2946924/70‑1 (MQ=255)
         gaatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                                      >  1:4062546/4‑70 (MQ=255)
         gaatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGa                                                     <  1:4452387/68‑1 (MQ=255)
           atCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                                                   <  1:1278153/70‑1 (MQ=255)
            tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcg                                                  <  1:4420731/71‑1 (MQ=255)
             cGATTGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                                                   <  1:715837/69‑1 (MQ=255)
              gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                                                 <  1:629530/70‑2 (MQ=255)
              gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                                                 <  1:4599494/70‑2 (MQ=255)
                 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)
                      ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                          <  1:3819662/53‑1 (MQ=255)
                      ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                          <  1:2147224/53‑1 (MQ=255)
                      ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                        <  1:2763703/71‑1 (MQ=255)
                       tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                           <  1:425581/51‑1 (MQ=255)
                         ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCgag                                     <  1:1555226/71‑1 (MQ=255)
                           cAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                                                        <  1:4210280/50‑1 (MQ=255)
                            aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAtt                                  <  1:1617157/71‑1 (MQ=255)
                             ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                 <  1:3627377/71‑1 (MQ=255)
                             ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                 <  1:4522319/71‑1 (MQ=255)
                             ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                 <  1:3087504/71‑1 (MQ=255)
                             ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGGCAGGATCGGGGGCGCGGCGAGAt                                   <  1:620104/69‑1 (MQ=255)
                              gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                <  1:2935298/68‑1 (MQ=255)
                              gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                <  1:1221193/68‑1 (MQ=255)
                              gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                <  1:3418387/68‑1 (MQ=255)
                                  aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                          <  1:3363218/41‑1 (MQ=38)
                                  aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                          <  1:1484898/41‑1 (MQ=38)
                                  aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                             <  1:4265230/70‑1 (MQ=255)
                                  aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                             >  1:4119517/1‑70 (MQ=255)
                                   tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                                      <  1:127545/44‑1 (MQ=38)
                                     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)
                                          aCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                                      >  1:986794/1‑37 (MQ=37)
                                          aCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                         <  1:207328/66‑1 (MQ=255)
                                              ttGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGa                           <  1:3587500/60‑1 (MQ=255)
                                               tGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                                                   >  1:1448843/1‑35 (MQ=37)
                                                 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)
                                                                    |                                                              
GCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAA  >  AP009048/3825040‑3825170

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: