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% 272.8 / NA 89V150V (GTG→GTApyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/56);  total (33/56)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGT  >  AP009048/3825028‑3825162
                                                                    |                                                                  
aCCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                                                                  <  1:2142433/71‑1 (MQ=255)
aCCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                                                                  <  1:2177988/71‑1 (MQ=255)
  cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTg                                                                 >  1:4917043/1‑70 (MQ=255)
  cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                                                >  1:2043124/1‑71 (MQ=255)
  cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                                                <  1:729446/71‑1 (MQ=255)
   gccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                                                               <  1:1471559/71‑1 (MQ=255)
   gccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                                                               <  1:1516390/71‑1 (MQ=255)
    ccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAtt                                                              <  1:1994214/71‑1 (MQ=255)
    ccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAtt                                                              <  1:2972446/71‑1 (MQ=255)
    ccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAtt                                                              <  1:217157/71‑1 (MQ=255)
    ccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAtt                                                              >  1:5151739/1‑71 (MQ=255)
      ggAACGGCGATTCGCGAATCGATCGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc                                                            <  1:3502057/71‑1 (MQ=255)
       gAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc                                                            <  1:1643196/70‑1 (MQ=255)
        aaCGGCGATTCGCGAATCGATGGATATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc                                                          <  1:2550445/71‑1 (MQ=255)
           ggCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCttg                                                       <  1:1073394/71‑3 (MQ=255)
            gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                                                      <  1:617120/71‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  >  1:4122728/1‑71 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  <  1:2572157/71‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  <  1:2181241/71‑1 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAtt                                                              >  1:263286/1‑58 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  >  1:2429780/1‑70 (MQ=255)
                  cgcgATTCGATGGAGATTATTCAGGCTTATGGCGCGTCGCTTGCTGGCGTATTGATTTCGCTTGATCGTc                                                 <  1:5207062/70‑1 (MQ=255)
                   gcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                               >  1:1929920/1‑71 (MQ=255)
                    cgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                              <  1:5131789/71‑1 (MQ=255)
                    cgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                              <  1:1031153/71‑1 (MQ=255)
                     gaatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGa                                             <  1:4606113/68‑1 (MQ=255)
                     gaatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGa                                             <  1:5022526/68‑1 (MQ=255)
                       atCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                                           <  1:281851/70‑1 (MQ=255)
                       atCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                                           <  1:2496723/70‑1 (MQ=255)
                        tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcg                                          >  1:3432087/1‑71 (MQ=255)
                         cGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCg                                                           >  1:654599/1‑53 (MQ=39)
                          gATGGAGATTATTCAGGCTAATGGCGCGACTCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                                                >  1:39516/1‑63 (MQ=255)
                          gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                                        >  1:1681475/1‑69 (MQ=255)
                          gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                                        <  1:2709218/71‑3 (MQ=255)
                           aTGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                                                               <  1:2143935/47‑1 (MQ=38)
                            tGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAg                                               >  1:1448637/1‑62 (MQ=255)
                            tGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGgggg                                      <  1:664616/71‑1 (MQ=255)
                              gagaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcg                                    <  1:3523770/71‑1 (MQ=255)
                                gaTTATTCAGGCTAATGGCGCGACGCTTGCTTGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGcgcg                                  <  1:2299867/71‑1 (MQ=255)
                                gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCttg                                                       <  1:2463265/50‑3 (MQ=255)
                                gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGaa                                            <  1:2286496/61‑1 (MQ=255)
                                gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                                        >  1:875163/1‑63 (MQ=255)
                                gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGcgcg                                  <  1:2872675/71‑1 (MQ=255)
                                gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGcgcg                                  <  1:818771/71‑1 (MQ=255)
                                gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGcgcg                                  <  1:1437815/71‑1 (MQ=255)
                                  ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                <  1:3381109/71‑1 (MQ=255)
                                  ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                >  1:1477384/1‑71 (MQ=255)
                                  ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcgg                                 <  1:3169786/70‑1 (MQ=255)
                                   tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                >  1:1422559/1‑70 (MQ=255)
                                   tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                >  1:3994035/1‑70 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc                                                            >  1:3829640/1‑40 (MQ=38)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                >  1:1515632/1‑68 (MQ=255)
                                     ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCgag                             <  1:45608/71‑1 (MQ=255)
                                      tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  <  1:325361/49‑1 (MQ=38)
                                       cAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  <  1:1360805/48‑1 (MQ=38)
                                        aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcg                                    >  1:1185276/1‑61 (MQ=255)
                                        aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAtt                          <  1:515600/71‑1 (MQ=255)
                                         ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                         <  1:2790021/71‑1 (MQ=255)
                                           ctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCg                       <  1:521735/69‑1 (MQ=255)
                                             aaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGAt                                                     >  1:5361619/1‑39 (MQ=37)
                                              aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAtt                          >  1:3967269/1‑65 (MQ=255)
                                              aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCgg                      <  1:1791900/69‑1 (MQ=255)
                                               tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                                         <  1:2773914/49‑2 (MQ=38)
                                               tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGgggg                                      <  1:385098/52‑1 (MQ=255)
                                               tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGTGGGCGCGGCGAGATTTCGGCg                    <  1:1652882/70‑1 (MQ=255)
                                               tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGa                   >  1:1507987/1‑71 (MQ=255)
                                               tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGCGCGCGGCGAGATTTCGGCGa                   >  1:103649/1‑71 (MQ=255)
                                                ggCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAt                  >  1:4317804/1‑71 (MQ=255)
                                                ggCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAt                  >  1:2012385/1‑71 (MQ=255)
                                                 gcgcGTCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                 <  1:1991209/71‑1 (MQ=255)
                                                 gcgcGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                 >  1:343675/1‑71 (MQ=255)
                                                     gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg             <  1:3307753/71‑1 (MQ=255)
                                                     gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg             <  1:3642596/71‑1 (MQ=255)
                                                     gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg             <  1:5249711/71‑1 (MQ=255)
                                                      aCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCga                              <  1:2857306/53‑1 (MQ=255)
                                                      aCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCg                       <  1:4466852/60‑1 (MQ=255)
                                                      aCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg             >  1:809314/1‑70 (MQ=255)
                                                      aCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg             >  1:803068/1‑70 (MQ=255)
                                                       cGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg             >  1:2926667/1‑69 (MQ=255)
                                                        gCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAg          >  1:2593382/1‑71 (MQ=255)
                                                           tGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt        >  1:4086254/1‑70 (MQ=255)
                                                            gCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt        <  1:4241194/69‑1 (MQ=255)
                                                            gCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt        <  1:4561962/69‑1 (MQ=255)
                                                             cTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                                         >  1:4007665/1‑34 (MQ=37)
                                                             cTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCga                              <  1:44209/46‑1 (MQ=38)
                                                              tGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGcgcg                                  <  1:2678767/41‑1 (MQ=37)
                                                              tGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAg          <  1:1654460/65‑1 (MQ=255)
                                                                gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGt  >  1:1835720/1‑71 (MQ=255)
                                                                gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGt  <  1:3075122/71‑1 (MQ=255)
                                                                    |                                                                  
ACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGT  >  AP009048/3825028‑3825162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: