Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,096 G→A 100% V150V (GTG→GTA pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,0960GA100.0% 222.2 / NA 72V150V (GTG→GTApyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/47);  total (25/47)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCG  >  AP009048/3825028‑3825161
                                                                    |                                                                 
aCCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                                                                 >  1:3431158/1‑71 (MQ=255)
  cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                                               >  1:212467/1‑71 (MQ=255)
      ggAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc                                                           <  1:2888825/71‑1 (MQ=255)
        aaCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc                                                         <  1:1641416/71‑1 (MQ=255)
         aCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt                                                        >  1:3202809/1‑71 (MQ=255)
          cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                       <  1:4522355/71‑2 (MQ=255)
          cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                       >  1:1341393/1‑70 (MQ=255)
          cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                       <  1:4295501/71‑2 (MQ=255)
           ggCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc                                                         <  1:162912/68‑1 (MQ=255)
           ggCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                                       >  1:1754405/1‑69 (MQ=255)
            gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt                                                        >  1:646567/1‑68 (MQ=255)
            gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                                                     >  1:3304369/1‑71 (MQ=255)
            gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                                                     >  1:2788483/1‑71 (MQ=255)
             cGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt                                                        >  1:3030462/1‑67 (MQ=255)
             cGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                                                     >  1:4403973/1‑70 (MQ=255)
              gATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                   >  1:2647793/1‑71 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                 <  1:2687081/71‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                 <  1:517736/71‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                 <  1:3742535/71‑1 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                 >  1:4349746/1‑70 (MQ=255)
                  cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc                                                           >  1:866565/1‑59 (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)
                        tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                                               >  1:2281312/1‑49 (MQ=38)
                        tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcg                                         <  1:4420731/71‑1 (MQ=255)
                         cGATTGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                                          <  1:715837/69‑1 (MQ=255)
                          gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                                        <  1:4599494/70‑2 (MQ=255)
                          gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                                        <  1:629530/70‑2 (MQ=255)
                           aTGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                                                                 <  1:3686640/44‑1 (MQ=38)
                             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)
                                gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCttg                                                      >  1:1876353/1‑48 (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)
                                   tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                                               >  1:1484997/1‑38 (MQ=37)
                                   tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                  <  1:425581/51‑1 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCgag                            <  1:1555226/71‑1 (MQ=255)
                                       cAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                                               <  1:4210280/50‑1 (MQ=255)
                                        aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCg                                                          <  1:1783954/38‑1 (MQ=37)
                                        aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAtt                         <  1:1617157/71‑1 (MQ=255)
                                         ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                        <  1:3627377/71‑1 (MQ=255)
                                         ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                        <  1:3087504/71‑1 (MQ=255)
                                         ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                        <  1:4522319/71‑1 (MQ=255)
                                         ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGGCAGGATCGGGGGCGCGGCGAGAt                          <  1:620104/69‑1 (MQ=255)
                                          gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                       <  1:2935298/68‑1 (MQ=255)
                                          gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                       <  1:3418387/68‑1 (MQ=255)
                                          gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                       <  1:1221193/68‑1 (MQ=255)
                                              aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                 <  1:3363218/41‑1 (MQ=38)
                                              aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                 <  1:1484898/41‑1 (MQ=38)
                                              aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                    >  1:4119517/1‑70 (MQ=255)
                                              aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                    <  1:4265230/70‑1 (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)
                                                                    |                                                                 
ACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCG  >  AP009048/3825028‑3825161

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: