Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,072 C→T 100% A142A (GCC→GCT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,0720CT100.0% 231.0 / NA 78A142A (GCC→GCTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/47);  total (31/47)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTT  >  AP009048/3825006‑3825139
                                                                  |                                                                   
aaTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGTTTCGCGAATCGATGGAGATTATTCAGGCTAATg                                                                 >  1:2773843/1‑71 (MQ=255)
aaTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAAt                                                                  <  1:2935818/70‑1 (MQ=255)
 aTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATgg                                                                >  1:975689/1‑71 (MQ=255)
 aTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATgg                                                                <  1:1547425/71‑1 (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)
           atgatgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATg                                                                 <  1:409352/60‑1 (MQ=255)
           atgatgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCtt                                                      <  1:1294167/71‑1 (MQ=255)
               tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                                                  >  1:4554703/1‑71 (MQ=255)
               tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                                                  <  1:3471435/71‑1 (MQ=255)
               tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                                                  >  1:4751754/1‑71 (MQ=255)
                gtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTgg                                                 >  1:3186607/1‑71 (MQ=255)
                 tgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGc                                                    <  1:1002157/67‑1 (MQ=255)
                 tgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTACTCAGGCTAATgg                                                                >  1:3625188/1‑55 (MQ=255)
                    tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgta                                             <  1:1666408/71‑2 (MQ=255)
                    tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgta                                             <  1:5083493/71‑2 (MQ=255)
                    tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgta                                             <  1:5143105/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:1420943/1‑71 (MQ=255)
                        cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                         >  1:3307288/1‑71 (MQ=255)
                         gccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                                        <  1:567019/71‑1 (MQ=255)
                         gccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                                        >  1:3659965/1‑71 (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)
                                           gaatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGa                      <  1:2676720/68‑1 (MQ=255)
                                            aatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGaa                     >  1:486598/3‑71 (MQ=255)
                                             atCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                    <  1:3211332/70‑1 (MQ=255)
                                              tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                            <  1:1901699/62‑1 (MQ=255)
                                              tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                         >  1:4035541/1‑65 (MQ=255)
                                              tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                         >  1:49131/1‑65 (MQ=255)
                                              tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcg                   <  1:4945953/71‑1 (MQ=255)
                                                gATGGAGATTATTCAGGCTAATGGCGGGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                 <  1:4502927/71‑3 (MQ=255)
                                                gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                 >  1:1258178/1‑69 (MQ=255)
                                                gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                  <  1:2912654/70‑2 (MQ=255)
                                                gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                  <  1:4375421/70‑2 (MQ=255)
                                                 aTGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                       <  1:2598413/64‑1 (MQ=255)
                                                  tGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGgggg               <  1:4910219/71‑1 (MQ=255)
                                                   ggAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc              >  1:2954676/1‑71 (MQ=255)
                                                   ggAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc              >  1:1833793/1‑71 (MQ=255)
                                                    gagaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                             >  1:2046396/1‑55 (MQ=255)
                                                     agaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcgc            >  1:5158911/1‑71 (MQ=255)
                                                        ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc         >  1:4599990/1‑71 (MQ=255)
                                                         tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                    <  1:1408462/59‑1 (MQ=255)
                                                          attCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTc                          >  1:2909494/1‑52 (MQ=255)
                                                           ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                         >  1:1824494/1‑52 (MQ=255)
                                                             cAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAt    >  1:4052875/1‑71 (MQ=255)
                                                              aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc                                   <  1:3474160/39‑1 (MQ=37)
                                                              aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAtt   <  1:592841/71‑1 (MQ=255)
                                                               ggCTAATGGCCCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt  <  1:2476450/71‑1 (MQ=255)
                                                                  |                                                                   
AATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTT  >  AP009048/3825006‑3825139

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: