Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,898 A→G 100% A84A (GCA→GCG pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8980AG100.0% 451.2 / NA 135A84A (GCA→GCGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (66/69);  total (66/69)
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 135 total aligned reads displayed.

TTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGT  >  AP009048/3824831‑3824964
                                                                   |                                                                  
ttGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                                                 >  1:1648294/1‑71 (MQ=255)
  gAGTTCGACCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCt                                                                  >  1:5155731/1‑68 (MQ=255)
   aGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCg                                                              <  1:7393854/71‑1 (MQ=255)
   aGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCg                                                              >  1:4732563/1‑71 (MQ=255)
     ttCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCg                                                              <  1:7503128/69‑1 (MQ=255)
       cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                           <  1:3417347/70‑1 (MQ=255)
       cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                           >  1:7597167/1‑70 (MQ=255)
       cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                           >  1:4655049/1‑70 (MQ=255)
       cGATCTGCTGTTTGGACCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                           >  1:1804610/1‑70 (MQ=255)
              ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                         <  1:4835809/65‑1 (MQ=255)
              ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                                    <  1:4278503/70‑3 (MQ=255)
              ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                                      <  1:3493566/68‑1 (MQ=255)
                  ttGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                               >  1:1452605/1‑70 (MQ=255)
                   tGGCCCTGCTTACAAAGGGATCCCGAATGCGACCACAACCGCCGTGGCGCTgg                                                                >  1:226497/1‑53 (MQ=255)
                    ggCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                               >  1:3637420/1‑68 (MQ=255)
                     gCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                            <  1:2612917/71‑1 (MQ=255)
                     gCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                            <  1:3284182/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)
                        cTGCTTACAAAGGGACCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                               <  1:2035985/65‑2 (MQ=255)
                         tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGc                                                          >  1:5619638/1‑53 (MQ=37)
                         tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                          >  1:5172867/1‑69 (MQ=255)
                         tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                              <  1:1032786/65‑3 (MQ=255)
                         tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                              <  1:4765372/65‑3 (MQ=255)
                          gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                                      >  1:5711953/1‑56 (MQ=255)
                          gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctgac                                       >  1:5771561/1‑69 (MQ=255)
                          gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                        >  1:6753898/1‑70 (MQ=255)
                           cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                       <  1:4412129/70‑1 (MQ=255)
                           cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                       <  1:1323528/70‑1 (MQ=255)
                           cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                       <  1:1206086/70‑1 (MQ=255)
                            ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                         <  1:526291/51‑1 (MQ=37)
                            ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                         <  1:2236/51‑1 (MQ=37)
                            ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                      >  1:3020865/1‑70 (MQ=255)
                            ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                      >  1:3578159/1‑70 (MQ=255)
                             tACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                                    >  1:6568033/1‑71 (MQ=255)
                              aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                                    <  1:490045/54‑3 (MQ=37)
                              aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                                    <  1:606600/54‑3 (MQ=37)
                              aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                                   >  1:2751745/1‑71 (MQ=255)
                               cAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacc                                           <  1:6556550/62‑1 (MQ=255)
                               cAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                                   <  1:5129861/70‑1 (MQ=255)
                                aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGc                                                               <  1:1748530/41‑1 (MQ=21)
                                aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATc                                                       >  1:2007642/1‑49 (MQ=37)
                                aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                 >  1:4233641/1‑71 (MQ=255)
                                  aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                           >  1:1274278/1‑43 (MQ=25)
                                  aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                               >  1:2793645/1‑71 (MQ=255)
                                    ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                             <  1:6024241/71‑1 (MQ=255)
                                      aTCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                         <  1:3094763/41‑1 (MQ=21)
                                         ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                          >  1:18395/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)
                                           gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                       >  1:6457176/1‑70 (MQ=255)
                                           gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                      <  1:309930/71‑1 (MQ=255)
                                            aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                 >  1:7086049/1‑59 (MQ=255)
                                             ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                         <  1:6134199/50‑1 (MQ=19)
                                             ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                               <  1:5684386/60‑1 (MQ=255)
                                             ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                    <  1:2633275/71‑1 (MQ=255)
                                             ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                    <  1:2632940/71‑1 (MQ=255)
                                             ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                    <  1:2375512/71‑1 (MQ=255)
                                               gcgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                    >  1:3020275/4‑69 (MQ=255)
                                                 gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                             >  1:6294539/2‑42 (MQ=16)
                                                 gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                             >  1:7031577/2‑42 (MQ=16)
                                                 gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                             >  1:166434/2‑42 (MQ=16)
                                                 gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                      <  1:4692585/48‑1 (MQ=25)
                                                 gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                      <  1:679610/48‑1 (MQ=25)
                                                 gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                      <  1:1621867/48‑1 (MQ=25)
                                                 gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                >  1:3129186/2‑69 (MQ=255)
                                                  accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                <  1:4546941/70‑3 (MQ=255)
                                                  accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc               <  1:4151007/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:5558287/71‑1 (MQ=255)
                                                    cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa             <  1:7888624/71‑1 (MQ=255)
                                                    cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa             <  1:4634213/71‑1 (MQ=255)
                                                    cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa             <  1:4097911/71‑1 (MQ=255)
                                                      ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa           <  1:919819/69‑1 (MQ=255)
                                                      ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa           <  1:2184884/69‑1 (MQ=255)
                                                      ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa             >  1:6701361/3‑69 (MQ=255)
                                                       taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                      <  1:5992662/58‑1 (MQ=255)
                                                       taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa             >  1:4752525/2‑68 (MQ=255)
                                                        aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa           >  1:7438602/1‑69 (MQ=255)
                                                         ccgccgTGGCGCTGGCGGAGCATCATgaccttgac                                            >  1:7144412/1‑35 (MQ=16)
                                                         ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTaa                           >  1:1097965/1‑52 (MQ=37)
                                                          cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc                                <  1:7795520/46‑1 (MQ=25)
                                                          cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTaa                           <  1:2856304/51‑1 (MQ=37)
                                                          cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                    <  1:3800974/58‑1 (MQ=255)
                                                          cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaag                   <  1:5910025/59‑1 (MQ=255)
                                                          cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAcc       <  1:5518582/71‑1 (MQ=255)
                                                           gccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCa      >  1:72968/4‑71 (MQ=255)
                                                           gccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCa      <  1:7023213/68‑1 (MQ=255)
                                                           gccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCa      >  1:3013685/4‑71 (MQ=255)
                                                            ccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCAc     >  1:7249427/3‑71 (MQ=255)
                                                            ccgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCAc     >  1:7585875/3‑71 (MQ=255)
                                                             cgtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa           >  1:4429544/2‑64 (MQ=255)
                                                                   |                                                                  
TTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGT  >  AP009048/3824831‑3824964

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: