Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,051 G→A 100% E135E (GAG→GAA pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,0510GA100.0% 246.8 / NA 83E135E (GAG→GAApyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/49);  total (34/49)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGG  >  AP009048/3824984‑3825118
                                                                   |                                                                   
gtAGCGCGTTACAAGGACGCGAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCg                                                                  <  1:2468413/70‑1 (MQ=255)
   gcgcGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGAt                                                                <  1:408937/70‑1 (MQ=255)
   gcgcGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATg                                                               <  1:4152363/71‑1 (MQ=255)
   gcgcGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATg                                                               >  1:3308056/1‑71 (MQ=255)
   gcgcGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATg                                                               >  1:3302362/1‑71 (MQ=255)
   gcgcGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATg                                                               >  1:2320400/1‑71 (MQ=255)
     gcgTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATg                                                               >  1:644200/1‑69 (MQ=255)
     gcgTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATg                                                               >  1:323184/1‑69 (MQ=255)
     gcgTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATg                                                               >  1:1663620/1‑69 (MQ=255)
       gTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGAt                                                                >  1:1515076/1‑66 (MQ=255)
       gTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGAt                                                                >  1:2490784/1‑66 (MQ=255)
         tACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGAt                                                                <  1:159269/64‑1 (MQ=255)
         tACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGat                                                          >  1:2519876/1‑70 (MQ=255)
         tACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGat                                                          >  1:2458255/1‑70 (MQ=255)
          aCAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGatta                                                        >  1:2739721/1‑71 (MQ=255)
           cAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattat                                                       <  1:3951075/71‑1 (MQ=255)
           cAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGatta                                                        <  1:1280690/70‑1 (MQ=255)
            aaGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattatt                                                      <  1:3279978/71‑1 (MQ=255)
            aaGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattatt                                                      >  1:216977/1‑71 (MQ=255)
              ggACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCa                                                    <  1:3228694/71‑1 (MQ=255)
               gACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAg                                                   <  1:3388867/71‑1 (MQ=255)
                aCGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAgg                                                  >  1:4432560/1‑71 (MQ=255)
                  gcgTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGgag                                                            >  1:2851943/1‑59 (MQ=255)
                  gcgTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGct                                                <  1:148953/71‑2 (MQ=255)
                  gcgTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGct                                                <  1:3758535/71‑2 (MQ=255)
                      aaTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGTTTCGCGAATCGATGGAGATTATTCAGGCTAATg                                            >  1:2773843/1‑71 (MQ=255)
                      aaTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAAt                                             <  1:2935818/70‑1 (MQ=255)
                       aTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATgg                                           <  1:1547425/71‑1 (MQ=255)
                       aTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATgg                                           >  1:975689/1‑71 (MQ=255)
                          cTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGgcg                                         >  1:2789711/1‑70 (MQ=255)
                             gTAGATGATGTGATCAGCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGcgcg                                       <  1:2111416/69‑1 (MQ=255)
                             gTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGcgcg                                       <  1:5110054/69‑1 (MQ=255)
                             gTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGcgcg                                       <  1:3002781/69‑1 (MQ=255)
                             gTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGAc                                     >  1:4288239/1‑71 (MQ=255)
                              tAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACg                                    <  1:4907272/71‑1 (MQ=255)
                                gatgatGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGc                                   <  1:3714681/70‑1 (MQ=255)
                                 atgatgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGct                                                <  1:4073618/56‑2 (MQ=255)
                                 atgatgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATg                                            <  1:409352/60‑1 (MQ=255)
                                 atgatgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCtt                                 <  1:1294167/71‑1 (MQ=255)
                                    atgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTc                                                     <  1:3140092/48‑1 (MQ=255)
                                     tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                             >  1:4554703/1‑71 (MQ=255)
                                     tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                             <  1:3471435/71‑1 (MQ=255)
                                     tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                             >  1:4751754/1‑71 (MQ=255)
                                     tgtgATCACCGCCGGAACGGCGATTCCCGAATCGATg                                                               >  1:1954837/1‑37 (MQ=37)
                                      gtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTgg                            >  1:3186607/1‑71 (MQ=255)
                                       tgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGc                               <  1:1002157/67‑1 (MQ=255)
                                       tgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTACTCAGGCTAATgg                                           >  1:3625188/1‑55 (MQ=255)
                                         aTCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGct                                                <  1:4781716/48‑2 (MQ=39)
                                         aTCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGct                                                <  1:1661184/48‑2 (MQ=39)
                                          tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgta                        <  1:5143105/71‑2 (MQ=255)
                                          tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgta                        <  1:1666408/71‑2 (MQ=255)
                                          tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgta                        <  1:5083493/71‑2 (MQ=255)
                                            aCCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCg                          <  1:990384/67‑1 (MQ=255)
                                            aCCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                      <  1:598691/71‑1 (MQ=255)
                                             ccgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                      <  1:4398871/70‑1 (MQ=255)
                                              cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGGTGGCGTATTGa                    >  1:890219/1‑71 (MQ=255)
                                              cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgt                         >  1:5067526/1‑66 (MQ=255)
                                              cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                    >  1:3307288/1‑71 (MQ=255)
                                              cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                    >  1:1420943/1‑71 (MQ=255)
                                               gccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                   >  1:3659965/1‑71 (MQ=255)
                                               gccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                   <  1:567019/71‑1 (MQ=255)
                                                gcgGTACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATg                                            <  1:2706552/44‑1 (MQ=37)
                                                ccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgt                         <  1:1176174/64‑1 (MQ=255)
                                                ccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgt                         <  1:1234933/64‑1 (MQ=255)
                                                 cgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATgg                                           >  1:3320328/1‑45 (MQ=38)
                                                   gAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCg               >  1:4421202/1‑71 (MQ=255)
                                                    aaCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc              <  1:2876051/71‑1 (MQ=255)
                                                     aCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCg               <  1:454432/69‑1 (MQ=255)
                                                      cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt            <  1:4540550/71‑2 (MQ=255)
                                                       ggCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCttg           <  1:4885902/71‑3 (MQ=255)
                                                       ggCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt            >  1:1591931/1‑69 (MQ=255)
                                                         cGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGAt         >  1:4939862/1‑71 (MQ=255)
                                                          gATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCg                          <  1:900332/53‑1 (MQ=39)
                                                          gATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc        <  1:4188034/71‑1 (MQ=255)
                                                            ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGc                           <  1:4451041/50‑1 (MQ=39)
                                                            ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt      <  1:1307027/71‑1 (MQ=255)
                                                            ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt      <  1:4841797/71‑1 (MQ=255)
                                                             tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt      <  1:2974398/70‑1 (MQ=255)
                                                             tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt      <  1:3794901/70‑1 (MQ=255)
                                                              cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTc     <  1:557728/70‑1 (MQ=255)
                                                               gcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc                >  1:2408887/1‑58 (MQ=255)
                                                                cgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa    >  1:2111461/1‑69 (MQ=255)
                                                                cgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg  <  1:4351201/71‑1 (MQ=255)
                                                                   |                                                                   
GTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGG  >  AP009048/3824984‑3825118

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: