Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,880 C→G 100% A78A (GCC→GCG pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8800CG100.0% 197.0 / NA 58A78A (GCC→GCGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (30/28);  total (30/28)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAA  >  AP009048/3824812‑3824946
                                                                    |                                                                  
gcattGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgac                                                                  <  1:1600689/68‑1 (MQ=255)
 cattGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacc                                                                 <  1:732104/69‑1 (MQ=255)
  attGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacca                                                                >  1:4037447/2‑71 (MQ=255)
   ttGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccac                                                               <  1:3313188/71‑1 (MQ=255)
   ttGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccac                                                               <  1:2469207/71‑1 (MQ=255)
     ggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacca                                                                >  1:4353015/1‑68 (MQ=255)
     ggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccact                                                              <  1:1173327/70‑2 (MQ=255)
      gtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccacta                                                             <  1:2385354/70‑3 (MQ=255)
        ggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacca                                                                >  1:4082331/1‑65 (MQ=255)
          aTTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGcc                                                        >  1:1773308/1‑70 (MQ=255)
                gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCcgtg                                                     <  1:4048882/68‑1 (MQ=255)
                gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcg                                                  >  1:2605870/1‑70 (MQ=255)
                  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)
                              ctgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcg                                                  >  1:1508981/1‑56 (MQ=255)
                              ctgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                    >  1:2036439/1‑71 (MQ=255)
                              ctgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                    >  1:2908663/1‑71 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                       <  1:1666986/65‑1 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCatg                                  <  1:4338451/70‑3 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                    <  1:4464989/68‑1 (MQ=255)
                                 ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                    <  1:4097109/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:402277/1‑62 (MQ=255)
                                      tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATc                                     >  1:1155934/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)
                                           cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGc                                                   <  1:980639/43‑1 (MQ=37)
                                           cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                       >  1:4154892/1‑71 (MQ=255)
                                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGa                                          <  1:3218598/49‑1 (MQ=37)
                                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                    <  1:3952206/71‑1 (MQ=255)
                                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                    >  1:4034216/1‑71 (MQ=255)
                                              cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                    <  1:2477700/71‑1 (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:4614601/71‑1 (MQ=255)
                                                        gATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTa          <  1:1466832/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)
                                                                    |                                                                  
GCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAA  >  AP009048/3824812‑3824946

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: