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% 211.0 / NA 70V150V (GTG→GTApyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/42);  total (28/42)
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:598691/71‑1 (MQ=255)
 ccgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                                                                  <  1:4398871/70‑1 (MQ=255)
  cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGGTGGCGTATTGa                                                                >  1:890219/1‑71 (MQ=255)
  cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                                                >  1:3307288/1‑71 (MQ=255)
  cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                                                >  1:1420943/1‑71 (MQ=255)
   gccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                                                               <  1:567019/71‑1 (MQ=255)
   gccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGAt                                                               >  1:3659965/1‑71 (MQ=255)
       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)
              gATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                    <  1:4188034/71‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  <  1:4841797/71‑1 (MQ=255)
                ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  <  1:1307027/71‑1 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  <  1:3794901/70‑1 (MQ=255)
                 tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                  <  1:2974398/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:49131/1‑65 (MQ=255)
                        tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                                                >  1:4035541/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:1833793/1‑71 (MQ=255)
                             ggAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc                                     >  1:2954676/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)
                                          gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                         <  1:2021689/67‑1 (MQ=255)
                                           ctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCg                       <  1:1901090/69‑1 (MQ=255)
                                           ctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCg                       >  1:1146799/3‑71 (MQ=255)
                                             aaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                        <  1:1762238/68‑1 (MQ=255)
                                             aaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCgg                      >  1:969761/1‑70 (MQ=255)
                                               tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                   <  1:3139225/39‑1 (MQ=37)
                                               tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcg                                          <  1:938154/48‑1 (MQ=39)
                                                ggCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                   <  1:4367227/38‑1 (MQ=38)
                                                  cgcgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                         <  1:2355370/62‑1 (MQ=255)
                                                   gcgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                 >  1:2399176/1‑69 (MQ=255)
                                                   gcgACGCTTGCTGGCGTATTGATTTCGCTTGATCCTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                 >  1:80286/1‑69 (MQ=255)
                                                    cgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAg              <  1:3163753/71‑1 (MQ=255)
                                                     gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg             <  1:4849352/71‑1 (MQ=255)
                                                     gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAg              >  1:353631/1‑70 (MQ=255)
                                                       cGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg             >  1:4973471/1‑69 (MQ=255)
                                                       cGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGa            <  1:2896662/70‑1 (MQ=255)
                                                            gCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt        >  1:1020674/1‑69 (MQ=255)
                                                            gCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTg       <  1:5123487/70‑1 (MQ=255)
                                                             cTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa      <  1:1316577/70‑1 (MQ=255)
                                                              tGGCGTATTGATTTCGCTTGATCGTCAGGAACGgggg                                      <  1:164915/37‑1 (MQ=25)
                                                                gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCTAGATTTCGGCGATTCAGGAAGTTGAGCGt  <  1:4206122/71‑1 (MQ=255)
                                                                gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCg   <  1:1903097/70‑1 (MQ=255)
                                                                    |                                                                  
ACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGT  >  AP009048/3825028‑3825162

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: