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% 249.9 / NA 74A84A (GCA→GCGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/37);  total (37/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGT  >  AP009048/3824830‑3824964
                                                                    |                                                                  
aTTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCt                                                                  <  1:3406179/71‑1 (MQ=255)
   gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTgg                                                                <  1:4170334/70‑1 (MQ=255)
   gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTgg                                                                <  1:884798/70‑1 (MQ=255)
   gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGc                                                               >  1:3684528/1‑71 (MQ=255)
    aGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGct                                                              <  1:1566706/71‑2 (MQ=255)
     gTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCgg                                                             <  1:1119110/71‑1 (MQ=255)
     gTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCgg                                                             <  1:369087/71‑1 (MQ=255)
        cGATCTGCTGTTTGGCCCTGATTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                           >  1:4340154/1‑70 (MQ=255)
        cGATCTGCTGTTTGGCCCTGATTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                           >  1:4345456/1‑70 (MQ=255)
            ctgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                                      >  1:2908663/1‑71 (MQ=255)
            ctgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                                      >  1:2036439/1‑71 (MQ=255)
               ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                         <  1:1666986/65‑1 (MQ=255)
               ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                                    <  1:4338451/70‑3 (MQ=255)
               ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                                      <  1:4097109/68‑1 (MQ=255)
               ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                                      <  1:4464989/68‑1 (MQ=255)
               ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                                   >  1:2549689/1‑71 (MQ=255)
                tgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                           >  1:4157981/1‑62 (MQ=255)
                tgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                         <  1:3294156/64‑1 (MQ=255)
                tgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                                                  >  1:535883/1‑71 (MQ=255)
                tgTTTCGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                                                  >  1:4275619/1‑71 (MQ=255)
                   ttGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                               >  1:914818/1‑70 (MQ=255)
                    tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATc                                                       >  1:1155934/1‑62 (MQ=255)
                    tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATc                                                       >  1:402277/1‑62 (MQ=255)
                    tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                               <  1:723396/70‑2 (MQ=255)
                     ggCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                              <  1:946871/70‑3 (MQ=255)
                      gCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGa                                                            <  1:1834914/55‑1 (MQ=255)
                       cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                                   <  1:1173242/63‑1 (MQ=255)
                         cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                         >  1:4154892/1‑71 (MQ=255)
                            cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGa                                                            <  1:3218598/49‑1 (MQ=37)
                            cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                      <  1:3952206/71‑1 (MQ=255)
                            cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                      <  1:2477700/71‑1 (MQ=255)
                            cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                      >  1:4034216/1‑71 (MQ=255)
                               aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                             >  1:4554861/1‑61 (MQ=255)
                               aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                              <  1:1502046/60‑3 (MQ=255)
                                cAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                         <  1:3650222/48‑1 (MQ=37)
                                 aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                                    >  1:3571818/1‑50 (MQ=37)
                                  aaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                                        >  1:3408831/1‑47 (MQ=25)
                                   aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCt                                                                  >  1:2958867/1‑36 (MQ=16)
                                   aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                               >  1:1022278/1‑71 (MQ=255)
                                   aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                               >  1:2749567/1‑71 (MQ=255)
                                     ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                               >  1:3127028/1‑52 (MQ=25)
                                      gATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTa                            <  1:1466832/71‑1 (MQ=255)
                                      gATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTa                            <  1:4614601/71‑1 (MQ=255)
                                         cccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                         >  1:4509904/1‑71 (MQ=255)
                                         cccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                         >  1:440643/1‑71 (MQ=255)
                                          ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCg                        >  1:962894/1‑71 (MQ=255)
                                            gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCat                                                     <  1:65220/40‑2 (MQ=21)
                                            gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                                   >  1:3427009/1‑58 (MQ=255)
                                            gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                       >  1:1161898/1‑70 (MQ=255)
                                             aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                         <  1:2744902/35‑1 (MQ=16)
                                             aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaa                     <  1:4526239/71‑1 (MQ=255)
                                             aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                      >  1:4159049/1‑70 (MQ=255)
                                              ttGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                    <  1:1225498/71‑1 (MQ=255)
                                                gcgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaag                   <  1:3971688/67‑1 (MQ=255)
                                                 cgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaaga                  <  1:2718398/68‑1 (MQ=255)
                                                   accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                >  1:2062596/1‑68 (MQ=255)
                                                    ccacTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                               <  1:3581350/38‑2 (MQ=16)
                                                    ccacTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                               <  1:595433/38‑2 (MQ=16)
                                                    ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc               >  1:1877022/1‑70 (MQ=255)
                                                    ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCa              <  1:45416/71‑1 (MQ=255)
                                                     cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagag                 <  1:1000237/67‑2 (MQ=255)
                                                       ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc                                <  1:2567848/48‑1 (MQ=25)
                                                       ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa           <  1:861457/69‑1 (MQ=255)
                                                        taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa             >  1:699391/2‑68 (MQ=255)
                                                         aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                         >  1:2805475/1‑55 (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)
                                                               gtgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                         >  1:1717083/1‑49 (MQ=38)
                                                                tgGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGt  <  1:1153678/71‑1 (MQ=255)
                                                                tgGCGCTGGCAGAGCATCATGACCTTGACCTGCCGTACTGCtt                              <  1:4030354/43‑1 (MQ=25)
                                                                    |                                                                  
ATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGT  >  AP009048/3824830‑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: