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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,9490AG100.0% 216.8 / NA 62E101E (GAA→GAGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/30);  total (32/30)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGG  >  AP009048/3824882‑3825013
                                                                   |                                                                
ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                                                >  1:1877022/1‑70 (MQ=255)
ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCa                                                               <  1:45416/71‑1 (MQ=255)
   ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                                            <  1:861457/69‑1 (MQ=255)
    taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                                              >  1:699391/2‑68 (MQ=255)
     aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGa                                                          >  1:3531944/1‑71 (MQ=255)
      ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                                                         >  1:301170/1‑71 (MQ=255)
       cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                                                         >  1:2433211/1‑70 (MQ=255)
         ccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                                            <  1:2300501/63‑1 (MQ=255)
          cgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGa                                                          >  1:787647/2‑66 (MQ=255)
          cgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACg                                                     >  1:3542628/2‑71 (MQ=255)
            tgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGt                                                   <  1:1153678/71‑1 (MQ=255)
              gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                                                  >  1:1464760/4‑70 (MQ=255)
              gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                                                  <  1:4145690/67‑1 (MQ=255)
                gcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGaa                                                <  1:1666346/69‑1 (MQ=255)
                 cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAgg                                              <  1:2756953/71‑1 (MQ=255)
                 cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACACGGTg                                                  <  1:4227687/66‑1 (MQ=255)
                   ggCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                                                  <  1:3488427/65‑1 (MQ=255)
                    gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg                                            <  1:4102944/70‑1 (MQ=255)
                    gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg                                            <  1:6348/70‑1 (MQ=255)
                     cGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcggc                                          <  1:4521222/71‑1 (MQ=255)
                       gAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCaa                                        >  1:2592744/1‑71 (MQ=255)
                        aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt                                       <  1:354690/71‑1 (MQ=255)
                           aTCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtg                                    <  1:580753/71‑1 (MQ=255)
                            tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtgg                                   >  1:1549055/1‑71 (MQ=255)
                            tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtgg                                   <  1:3171544/71‑1 (MQ=255)
                              atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                                            >  1:1827565/3‑44 (MQ=38)
                              atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGt                                                   >  1:779067/3‑53 (MQ=255)
                              atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt                                 >  1:3935334/3‑71 (MQ=255)
                              atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt                                 >  1:180533/3‑71 (MQ=255)
                                gaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCa                                                               <  1:4384097/39‑1 (MQ=38)
                                 accttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCaa                                        <  1:4358095/61‑1 (MQ=255)
                                 accttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttgg                               <  1:4422715/70‑1 (MQ=255)
                                 accttgacctgCCGTAATGCTTTAACCGCAAAGAGGCAAAAGACCAc                                                      >  1:419034/1‑47 (MQ=38)
                                     tgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTa                             >  1:4434094/2‑68 (MQ=255)
                                     tgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAgcg                          >  1:3633362/2‑71 (MQ=255)
                                     tgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACAACGGTGAAGGCGGCAATc                                      >  1:4358478/2‑59 (MQ=39)
                                      gacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAgc                           <  1:1223983/69‑1 (MQ=255)
                                        cctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggt                                  >  1:1917802/1‑60 (MQ=255)
                                        cctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGt                       >  1:484505/1‑71 (MQ=255)
                                         ctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGtt                      <  1:4046602/71‑1 (MQ=255)
                                           gCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTa                     >  1:2787669/1‑70 (MQ=255)
                                           gCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTa                     >  1:4372421/1‑70 (MQ=255)
                                            ccGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTAc                    >  1:4409848/1‑70 (MQ=255)
                                             cGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACaa                  <  1:3609967/71‑1 (MQ=255)
                                             cGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACaa                  <  1:1195127/71‑1 (MQ=255)
                                             cGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACaa                  <  1:2786256/71‑1 (MQ=255)
                                             cGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACaa                  <  1:4044171/71‑1 (MQ=255)
                                              gTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCaa                                        <  1:1992013/48‑1 (MQ=255)
                                              gTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAg                 <  1:2575451/71‑1 (MQ=255)
                                              gTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAg                 <  1:74315/71‑1 (MQ=255)
                                                aCTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGa               >  1:3067986/1‑71 (MQ=255)
                                                 cTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAc              <  1:2534745/71‑1 (MQ=255)
                                                   gCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATc                                      >  1:665524/1‑45 (MQ=255)
                                                      ttAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGt          >  1:2692478/1‑70 (MQ=255)
                                                       tAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGa               >  1:455540/1‑64 (MQ=255)
                                                       tAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTa         <  1:1475214/70‑1 (MQ=255)
                                                         aCCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAAt       >  1:1581771/1‑70 (MQ=255)
                                                         aCCGCAAAGAGCCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTa         <  1:1358585/68‑1 (MQ=255)
                                                          ccGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggt                                  >  1:2651042/1‑42 (MQ=255)
                                                             cAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTgg  >  1:4022381/1‑71 (MQ=255)
                                                             cAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTgg  >  1:679173/1‑71 (MQ=255)
                                                              aaaGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttgg                               >  1:1033526/1‑41 (MQ=255)
                                                                agaggCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTgg  >  1:3602123/1‑68 (MQ=255)
                                                                agaggCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTgg  >  1:2251706/1‑68 (MQ=255)
                                                                   |                                                                
CCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGG  >  AP009048/3824882‑3825013

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: