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% 249.0 / NA 83L154L (CTC→CTTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/48);  total (35/48)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCA  >  AP009048/3825040‑3825175
                                                                    |                                                                   
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)
     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)
              gATGGAGATTATTCAGGCTAATGGCGCGACTCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                                                             >  1:39516/1‑63 (MQ=255)
              gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                                                     >  1:1681475/1‑69 (MQ=255)
              gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                                                     <  1:2709218/71‑3 (MQ=255)
                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)
                    gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGaa                                                         <  1:2286496/61‑1 (MQ=255)
                    gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                                                     >  1:875163/1‑63 (MQ=255)
                    gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGcgcg                                               <  1:2872675/71‑1 (MQ=255)
                    gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGcgcg                                               <  1:1437815/71‑1 (MQ=255)
                    gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGcgcg                                               <  1:818771/71‑1 (MQ=255)
                      ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                             >  1:1477384/1‑71 (MQ=255)
                      ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                             <  1:3381109/71‑1 (MQ=255)
                      ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcgg                                              <  1:3169786/70‑1 (MQ=255)
                       tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                             >  1:3994035/1‑70 (MQ=255)
                       tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                             >  1:1422559/1‑70 (MQ=255)
                         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:2012385/1‑71 (MQ=255)
                                    ggCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAt                               >  1:4317804/1‑71 (MQ=255)
                                     gcgcGTCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                              <  1:1991209/71‑1 (MQ=255)
                                     gcgcGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                              >  1:343675/1‑71 (MQ=255)
                                         gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg                          <  1:5249711/71‑1 (MQ=255)
                                         gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg                          <  1:3307753/71‑1 (MQ=255)
                                         gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg                          <  1:3642596/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)
                                                        attGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGt               <  1:4882202/66‑1 (MQ=255)
                                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTg                    <  1:3523858/61‑1 (MQ=255)
                                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTa          >  1:2742988/1‑71 (MQ=255)
                                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTa          >  1:1903922/1‑71 (MQ=255)
                                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTa          >  1:4895025/1‑71 (MQ=255)
                                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTa          >  1:861699/1‑71 (MQ=255)
                                                          tGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTa          >  1:4395987/1‑70 (MQ=255)
                                                           gATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt                     >  1:3177528/1‑58 (MQ=255)
                                                           gATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTAc         <  1:850844/70‑1 (MQ=255)
                                                            aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGAtt           <  1:2452891/67‑1 (MQ=255)
                                                            aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACaa       <  1:4291460/71‑1 (MQ=255)
                                                              ttCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACt     >  1:449507/1‑71 (MQ=255)
                                                                 gCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGc   >  1:28862/1‑70 (MQ=255)
                                                                 gCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCa  <  1:922962/71‑1 (MQ=255)
                                                                    |                                                                   
GCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCA  >  AP009048/3825040‑3825175

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: