Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,919 G→T 100% L91L (CTG→CTT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,9190GT100.0% 215.6 / NA 71L91L (CTG→CTTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (36/35);  total (36/35)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTG  >  AP009048/3824852‑3824983
                                                                   |                                                                
gCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                                               <  1:1032726/71‑1 (MQ=255)
   cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                            >  1:1407351/1‑71 (MQ=255)
   cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacc                                                              <  1:3064459/69‑1 (MQ=255)
    tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                                           >  1:236181/1‑71 (MQ=255)
    tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                                           >  1:1571272/1‑71 (MQ=255)
     gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                            >  1:5343904/1‑69 (MQ=255)
     gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                                          >  1:2630670/1‑71 (MQ=255)
      cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                                               >  1:5163541/1‑65 (MQ=255)
      cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                                           <  1:4735984/69‑1 (MQ=255)
      cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                         <  1:2728394/71‑1 (MQ=255)
      cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                         >  1:2393895/1‑71 (MQ=255)
         aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGt                                                        <  1:2199800/69‑1 (MQ=255)
         aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                                                       <  1:1416437/70‑1 (MQ=255)
         aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                                                       <  1:3888244/70‑1 (MQ=255)
         aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                                                      >  1:2931801/1‑71 (MQ=255)
          cAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACt                                                     <  1:3402405/71‑1 (MQ=255)
           aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACt                                                     >  1:665569/1‑70 (MQ=255)
           aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                                    <  1:1187243/71‑1 (MQ=255)
            aaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc                                                   >  1:249454/1‑71 (MQ=255)
             aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACtt                                                    <  1:4805076/69‑2 (MQ=255)
             aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                  >  1:5325971/1‑71 (MQ=255)
               ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                                    >  1:1938844/1‑67 (MQ=255)
               ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCtt                                                 >  1:1682255/1‑70 (MQ=255)
                gATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                                >  1:914967/1‑70 (MQ=255)
                   cccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                                <  1:3829576/67‑1 (MQ=255)
                   cccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                                            >  1:4560910/1‑71 (MQ=255)
                   cccGATTGCGACCACTACCGCCGTGGCGCTAGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                                            >  1:1972628/1‑71 (MQ=255)
                   cccGATTGCGACCACTACCGCCGTGGCGCTAGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                                            >  1:1966408/1‑71 (MQ=255)
                    ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTa                                               <  1:4371101/67‑1 (MQ=255)
                      gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                          >  1:4813155/1‑70 (MQ=255)
                      gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                          >  1:5027929/1‑70 (MQ=255)
                      gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                          >  1:2878316/1‑70 (MQ=255)
                      gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                          >  1:1331533/1‑70 (MQ=255)
                       aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                                                      <  1:1848147/57‑1 (MQ=255)
                              ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                         <  1:4028579/47‑1 (MQ=25)
                              ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                         <  1:2239891/47‑1 (MQ=25)
                              ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaaga                                     >  1:1097051/1‑67 (MQ=255)
                              ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                  >  1:400672/1‑70 (MQ=255)
                               cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaaga                                     >  1:3048149/1‑66 (MQ=255)
                               cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                <  1:3644079/71‑1 (MQ=255)
                               cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                <  1:1886298/71‑1 (MQ=255)
                                 ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                              <  1:3964653/69‑1 (MQ=255)
                                 ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                <  1:2837012/67‑1 (MQ=255)
                                   aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                                                       >  1:1973744/1‑44 (MQ=25)
                                   aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTaa                                              <  1:1569717/53‑1 (MQ=37)
                                    ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAc                                             >  1:4929748/1‑53 (MQ=37)
                                    ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGa                            >  1:4187187/1‑70 (MQ=255)
                                    ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                           >  1:2466350/1‑71 (MQ=255)
                                    ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                           >  1:1156487/1‑71 (MQ=255)
                                     cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                  <  1:1488623/47‑1 (MQ=25)
                                     cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCCTACTGCTTTAACCGCAAAGAGGCaaaa                              <  1:2610851/67‑1 (MQ=255)
                                      gccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCa                         <  1:4970944/68‑1 (MQ=255)
                                       ccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaaga                                     >  1:3171140/3‑58 (MQ=255)
                                            gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                    >  1:1088947/4‑70 (MQ=255)
                                            gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                    <  1:3004996/67‑1 (MQ=255)
                                            gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGa                   <  1:299146/68‑1 (MQ=255)
                                             cgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                    <  1:580835/67‑1 (MQ=255)
                                              gcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaaga                                     <  1:1779155/50‑1 (MQ=255)
                                              gcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaag                                      <  1:3097260/49‑1 (MQ=255)
                                               cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAgg                <  1:1131829/71‑1 (MQ=255)
                                               cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAgg                <  1:4706368/71‑1 (MQ=255)
                                                tGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCATAGAGGCAAAAGACCACGGTGAAggc               >  1:1031417/1‑71 (MQ=255)
                                                tGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggc               >  1:2513449/1‑71 (MQ=255)
                                                tGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggc               >  1:5191957/1‑71 (MQ=255)
                                                   cGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcggc            <  1:2471079/71‑1 (MQ=255)
                                                    ggAGCATCATGACCTTGACCTGCCGTACTGCTTTaa                                              >  1:3925094/1‑36 (MQ=37)
                                                    ggAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggc               <  1:3801430/67‑1 (MQ=255)
                                                      aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt         >  1:2243833/1‑71 (MQ=255)
                                                      aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt         >  1:2471032/1‑71 (MQ=255)
                                                      aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt         >  1:4499200/1‑71 (MQ=255)
                                                       gCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATc        <  1:818987/71‑1 (MQ=255)
                                                        catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCt       <  1:3798469/71‑1 (MQ=255)
                                                         aTCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGAGGCAATCtg      <  1:3096047/71‑1 (MQ=255)
                                                          tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACgg                      <  1:568617/54‑1 (MQ=255)
                                                          tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtgg     <  1:526134/71‑1 (MQ=255)
                                                          tCATGAACTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAcc                          >  1:691059/1‑50 (MQ=255)
                                                            atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCaaaa                              >  1:2946090/3‑44 (MQ=38)
                                                            atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGaa                  >  1:39058/3‑56 (MQ=255)
                                                            atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt   <  1:3949909/69‑1 (MQ=255)
                                                            atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt   <  1:4933271/69‑1 (MQ=255)
                                                            atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt   >  1:3490259/3‑71 (MQ=255)
                                                            atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggt    <  1:2792591/68‑1 (MQ=255)
                                                             tgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttg  >  1:5217725/2‑71 (MQ=255)
                                                             tgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttg  <  1:4696915/70‑1 (MQ=255)
                                                              gaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACAACGGTGAAGGCGGCAATCtggtt   >  1:840854/1‑69 (MQ=255)
                                                                   |                                                                
GCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTG  >  AP009048/3824852‑3824983

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: