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% 319.5 / NA 106L91L (CTG→CTTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (46/60);  total (46/60)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 140 total aligned reads displayed.

GCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTA  >  AP009048/3824852‑3824986
                                                                   |                                                                   
gCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                                                  <  1:2612917/71‑1 (MQ=255)
gCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                                                  <  1:6096595/71‑1 (MQ=255)
  ccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                                                <  1:3231018/71‑1 (MQ=255)
   cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                               >  1:3735880/1‑71 (MQ=255)
    tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                                                >  1:5172867/1‑69 (MQ=255)
    tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctga                                                              >  1:7863963/1‑70 (MQ=255)
     gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctgac                                                             >  1:5771561/1‑69 (MQ=255)
     gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                                             >  1:5187653/1‑71 (MQ=255)
     gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                                              >  1:6753898/1‑70 (MQ=255)
      cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                                             <  1:4412129/70‑1 (MQ=255)
      cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                                             <  1:1323528/70‑1 (MQ=255)
       ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                            >  1:3578159/1‑70 (MQ=255)
        tACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                                                          >  1:6568033/1‑71 (MQ=255)
         aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                                                         >  1:2751745/1‑71 (MQ=255)
          cAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacc                                                                 <  1:6556550/62‑1 (MQ=255)
           aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                                       >  1:4233641/1‑71 (MQ=255)
               ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                                   <  1:2743880/71‑1 (MQ=255)
               ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                                   <  1:6024241/71‑1 (MQ=255)
                gATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTa                                                  <  1:4548058/71‑1 (MQ=255)
                    ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                                                >  1:2182999/1‑53 (MQ=25)
                    ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                                                >  1:18395/1‑53 (MQ=25)
                    ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                                                >  1:314234/1‑53 (MQ=25)
                    ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                                                >  1:3075587/1‑53 (MQ=25)
                    ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                                               <  1:4300060/70‑1 (MQ=255)
                    ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCg                                              >  1:5451746/1‑71 (MQ=255)
                     cGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                     >  1:6355198/1‑63 (MQ=255)
                      gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACt                                                        <  1:6144929/59‑1 (MQ=255)
                      gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                                            <  1:309930/71‑1 (MQ=255)
                       aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                                       >  1:7086049/1‑59 (MQ=255)
                       aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                                       >  1:574762/1‑59 (MQ=255)
                       aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                                       >  1:937073/1‑59 (MQ=255)
                       aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                                            >  1:7111877/1‑70 (MQ=255)
                        ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                               <  1:2438601/50‑1 (MQ=19)
                        ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                     <  1:2963052/60‑1 (MQ=255)
                        ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                     <  1:5684386/60‑1 (MQ=255)
                        ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                                          <  1:2633275/71‑1 (MQ=255)
                        ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                                          <  1:2375512/71‑1 (MQ=255)
                          gcgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTTACCGCa                                            <  1:15634/64‑1 (MQ=255)
                          gcgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                                          >  1:3020275/4‑69 (MQ=255)
                           cgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagag                                       <  1:6078481/69‑2 (MQ=255)
                            gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                                                   >  1:6294539/2‑42 (MQ=16)
                            gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                                                   >  1:7031577/2‑42 (MQ=16)
                            gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                                                   >  1:166434/2‑42 (MQ=16)
                            gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                                                   >  1:1739204/2‑42 (MQ=16)
                            gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                            <  1:1621867/48‑1 (MQ=25)
                             accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                     <  1:4151007/71‑1 (MQ=255)
                             accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                     <  1:1263178/71‑1 (MQ=255)
                              ccacTACCGCCGTGGCGCTGGCGGATCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCa                                    <  1:5470260/71‑1 (MQ=255)
                              ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                                                         >  1:2438393/1‑50 (MQ=25)
                              ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                             <  1:5545807/62‑1 (MQ=255)
                              ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                                      >  1:4185384/1‑67 (MQ=255)
                              ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                     >  1:5945727/1‑70 (MQ=255)
                               cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                   <  1:2034067/71‑1 (MQ=255)
                               cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                   <  1:4634213/71‑1 (MQ=255)
                               cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                   <  1:4097911/71‑1 (MQ=255)
                                actaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaa                                  <  1:5582654/68‑1 (MQ=255)
                                 ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                 <  1:919819/69‑1 (MQ=255)
                                 ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                 <  1:2184884/69‑1 (MQ=255)
                                 ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                 <  1:6433775/69‑1 (MQ=255)
                                 ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                   >  1:6701361/3‑69 (MQ=255)
                                  taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                                   <  1:5973533/51‑1 (MQ=37)
                                  taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                                   <  1:1885410/51‑1 (MQ=37)
                                  taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                   >  1:4752525/2‑68 (MQ=255)
                                  taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                   >  1:7507720/2‑68 (MQ=255)
                                   aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagag                                       >  1:7300140/1‑62 (MQ=255)
                                    ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTaa                                                 >  1:1097965/1‑52 (MQ=37)
                                     cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc                                                      <  1:7795520/46‑1 (MQ=25)
                                     cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTaa                                                 <  1:2856304/51‑1 (MQ=37)
                                     cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                                          <  1:4701630/58‑1 (MQ=255)
                                     cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                                          <  1:3800974/58‑1 (MQ=255)
                                     cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAcc                             <  1:5518582/71‑1 (MQ=255)
                                      gccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCa                            >  1:3013685/4‑71 (MQ=255)
                                       ccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTaa                                                 >  1:3452781/3‑49 (MQ=37)
                                       ccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCAc                           >  1:7249427/3‑71 (MQ=255)
                                        cgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                 >  1:4429544/2‑64 (MQ=255)
                                          tgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGt                        <  1:3660262/71‑1 (MQ=255)
                                             cgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGaa                     <  1:2612353/69‑1 (MQ=255)
                                              gcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACgg                         >  1:7616592/2‑66 (MQ=255)
                                              gcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAg                    >  1:3183348/2‑71 (MQ=255)
                                                 ggCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaaga                                        >  1:6966821/1‑48 (MQ=255)
                                                  gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg                >  1:1143966/1‑71 (MQ=255)
                                                    ggAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaag                                         <  1:1712235/44‑1 (MQ=38)
                                                      aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                                      >  1:3487897/1‑43 (MQ=255)
                                                      aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                                      >  1:7175297/1‑43 (MQ=255)
                                                      aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAgg                   >  1:3823134/1‑64 (MQ=255)
                                                      aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt            >  1:4462122/1‑71 (MQ=255)
                                                      aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt            <  1:5008932/71‑1 (MQ=255)
                                                       gCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACggg                        <  1:512670/58‑2 (MQ=255)
                                                          tCATGACGTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggc                  <  1:7379329/61‑1 (MQ=255)
                                                          tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                 <  1:7348104/46‑1 (MQ=255)
                                                          tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                 <  1:5443309/46‑1 (MQ=255)
                                                          tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                   <  1:693186/44‑1 (MQ=38)
                                                           catgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggt       >  1:6664949/4‑71 (MQ=255)
                                                            atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                 >  1:283488/3‑44 (MQ=38)
                                                             tgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttg     >  1:5243238/2‑71 (MQ=255)
                                                             tgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggt       >  1:6135328/2‑69 (MQ=255)
                                                              gaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttgg    >  1:750048/1‑71 (MQ=255)
                                                               accttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttgg    >  1:1134962/1‑70 (MQ=255)
                                                               accttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt      >  1:5758100/1‑68 (MQ=255)
                                                                ccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTa  <  1:4770409/71‑1 (MQ=255)
                                                                   |                                                                   
GCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTA  >  AP009048/3824852‑3824986

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: