Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,892 T→C 100% A82A (GCT→GCC pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8920TC100.0% 446.6 / NA 133A82A (GCT→GCCpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (61/72);  total (61/72)
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 148 total aligned reads displayed.

GCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACC  >  AP009048/3824828‑3824959
                                                                |                                                                   
gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcg                                                               >  1:2934536/1‑70 (MQ=255)
gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcg                                                               <  1:3230129/71‑2 (MQ=255)
 cATTGAGTTCGATCTGCTGTTTGGCCATGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcgc                                                              >  1:6920961/1‑69 (MQ=255)
        ttCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCg                                                         <  1:7503128/69‑1 (MQ=255)
          cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                      >  1:7597167/1‑70 (MQ=255)
          cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                      <  1:3417347/70‑1 (MQ=255)
          cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                      >  1:4655049/1‑70 (MQ=255)
          cGATCTGCTGTTTGGACCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                      >  1:1804610/1‑70 (MQ=255)
            aTCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCcgtg                                                                  <  1:4663369/56‑1 (MQ=255)
                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                    <  1:4835809/65‑1 (MQ=255)
                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                               <  1:4278503/70‑3 (MQ=255)
                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                                 <  1:3493566/68‑1 (MQ=255)
                     ttGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCcgt                                                                   >  1:70043/1‑46 (MQ=38)
                     ttGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                          >  1:1452605/1‑70 (MQ=255)
                      tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCcgtg                                                                  >  1:1006467/1‑46 (MQ=38)
                       ggCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                          >  1:3637420/1‑68 (MQ=255)
                        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)
                           cTGCTTACAAAGGGACCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                          <  1:2035985/65‑2 (MQ=255)
                            tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCcgt                                                                   <  1:1921733/39‑1 (MQ=25)
                            tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGc                                                     >  1:5619638/1‑53 (MQ=37)
                            tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                     >  1:5172867/1‑69 (MQ=255)
                            tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctga                                   >  1:7863963/1‑70 (MQ=255)
                            tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                         <  1:4765372/65‑3 (MQ=255)
                            tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                         <  1:1032786/65‑3 (MQ=255)
                             gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctgac                                  >  1:5771561/1‑69 (MQ=255)
                             gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                  >  1:5187653/1‑71 (MQ=255)
                             gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                   >  1:6753898/1‑70 (MQ=255)
                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCcgt                                                                   >  1:1691111/1‑37 (MQ=25)
                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCcgt                                                                   >  1:1065733/1‑37 (MQ=25)
                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                  <  1:4412129/70‑1 (MQ=255)
                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                  <  1:1206086/70‑1 (MQ=255)
                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                  <  1:1323528/70‑1 (MQ=255)
                               ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                    <  1:2236/51‑1 (MQ=37)
                               ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                    <  1:526291/51‑1 (MQ=37)
                               ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                 >  1:3020865/1‑70 (MQ=255)
                                tACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGc                                                                <  1:522571/38‑1 (MQ=25)
                                tACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                               >  1:6568033/1‑71 (MQ=255)
                                 aCAAAGGGATCCCGATTGCGACCACTACCGCcgtg                                                                  >  1:7532440/1‑35 (MQ=25)
                                 aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcg                                                               <  1:7390472/38‑2 (MQ=25)
                                 aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                               <  1:236852/54‑3 (MQ=37)
                                 aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                               <  1:490045/54‑3 (MQ=37)
                                 aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                              >  1:2751745/1‑71 (MQ=255)
                                  cAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                              <  1:5129861/70‑1 (MQ=255)
                                   aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGc                                                          <  1:1748530/41‑1 (MQ=21)
                                   aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                                   >  1:5148113/1‑48 (MQ=37)
                                     aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                      >  1:549977/1‑43 (MQ=25)
                                     aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                          >  1:2793645/1‑71 (MQ=255)
                                       ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                        <  1:6024241/71‑1 (MQ=255)
                                       ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                        <  1:2743880/71‑1 (MQ=255)
                                            ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgaccttgacct                                     >  1:3075587/1‑53 (MQ=25)
                                            ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                    <  1:4300060/70‑1 (MQ=255)
                                            ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCg                   >  1:4198136/1‑71 (MQ=255)
                                              gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                  >  1:6457176/1‑70 (MQ=255)
                                              gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                 <  1:309930/71‑1 (MQ=255)
                                               aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                            >  1:574762/1‑59 (MQ=255)
                                               aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                            >  1:6292552/1‑59 (MQ=255)
                                               aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                            >  1:937073/1‑59 (MQ=255)
                                               aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                            >  1:7086049/1‑59 (MQ=255)
                                               aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                 >  1:7111877/1‑70 (MQ=255)
                                                ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                    <  1:2438601/50‑1 (MQ=19)
                                                ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                    <  1:6134199/50‑1 (MQ=19)
                                                ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                          <  1:2963052/60‑1 (MQ=255)
                                                ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                          <  1:5684386/60‑1 (MQ=255)
                                                ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa               <  1:2375512/71‑1 (MQ=255)
                                                ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa               <  1:2633275/71‑1 (MQ=255)
                                                ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa               <  1:2632940/71‑1 (MQ=255)
                                                  gcgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTTACCGCa                 <  1:15634/64‑1 (MQ=255)
                                                   cgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagag            <  1:6078481/69‑2 (MQ=255)
                                                    gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                        >  1:166434/2‑42 (MQ=16)
                                                    gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                        >  1:6294539/2‑42 (MQ=16)
                                                    gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                 <  1:679610/48‑1 (MQ=25)
                                                    gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                 <  1:4040841/48‑1 (MQ=25)
                                                    gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                 <  1:1621867/48‑1 (MQ=25)
                                                    gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                 <  1:4692585/48‑1 (MQ=25)
                                                    gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg           >  1:3129186/2‑69 (MQ=255)
                                                     accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTaa                      >  1:6637091/1‑59 (MQ=255)
                                                     accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg           <  1:4546941/70‑3 (MQ=255)
                                                     accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc          <  1:1263178/71‑1 (MQ=255)
                                                      ccacTACCGCCGTGGCGCTGGCGGATCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCa         <  1:5470260/71‑1 (MQ=255)
                                                      ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                  <  1:5545807/62‑1 (MQ=255)
                                                      ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg           >  1:4185384/1‑67 (MQ=255)
                                                      ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc          >  1:5945727/1‑70 (MQ=255)
                                                       cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa        <  1:4634213/71‑1 (MQ=255)
                                                       cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa        <  1:4097911/71‑1 (MQ=255)
                                                       cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa        <  1:5558287/71‑1 (MQ=255)
                                                       cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa        <  1:2034067/71‑1 (MQ=255)
                                                         ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      <  1:919819/69‑1 (MQ=255)
                                                         ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      <  1:6433775/69‑1 (MQ=255)
                                                         ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      <  1:2184884/69‑1 (MQ=255)
                                                          taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                        <  1:1885410/51‑1 (MQ=37)
                                                          taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                 <  1:5992662/58‑1 (MQ=255)
                                                          taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa        >  1:4752525/2‑68 (MQ=255)
                                                           aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagag            >  1:7300140/1‑62 (MQ=255)
                                                           aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa      >  1:7438602/1‑69 (MQ=255)
                                                             cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc                           <  1:7795520/46‑1 (MQ=25)
                                                             cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa               <  1:4701630/58‑1 (MQ=255)
                                                             cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaag              <  1:5910025/59‑1 (MQ=255)
                                                                |                                                                   
GCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACC  >  AP009048/3824828‑3824959

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: