Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,913 C→T 100% H89H (CAC→CAT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,9130CT100.0% 193.2 / NA 65H89H (CAC→CATpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/35);  total (30/35)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGG  >  AP009048/3824845‑3824980
                                                                    |                                                                   
ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                                                   >  1:2549689/1‑71 (MQ=255)
 tgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                                                                  >  1:535883/1‑71 (MQ=255)
 tgTTTCGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgac                                                                  >  1:4275619/1‑71 (MQ=255)
    ttGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                                               >  1:914818/1‑70 (MQ=255)
     tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                                               <  1:723396/70‑2 (MQ=255)
      ggCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                                              <  1:946871/70‑3 (MQ=255)
        cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                                                   <  1:1173242/63‑1 (MQ=255)
          cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                         >  1:4154892/1‑71 (MQ=255)
             cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                      <  1:2477700/71‑1 (MQ=255)
             cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                      <  1:3952206/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)
                    aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                               >  1:2749567/1‑71 (MQ=255)
                    aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                               >  1:1022278/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:440643/1‑71 (MQ=255)
                          cccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                                         >  1:4509904/1‑71 (MQ=255)
                           ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCg                                        >  1:962894/1‑71 (MQ=255)
                             gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                                                   >  1:3427009/1‑58 (MQ=255)
                             gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                       >  1:1161898/1‑70 (MQ=255)
                              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)
                                                   gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                 <  1:4145690/67‑1 (MQ=255)
                                                   gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                 >  1:1464760/4‑70 (MQ=255)
                                                     gcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGaa               <  1:1666346/69‑1 (MQ=255)
                                                      cTGGCGGAGCATCATGTCCTTGACCTGCCGTACTGCTTTAACCGCaaa                                    <  1:3357981/48‑1 (MQ=255)
                                                      cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAgg             <  1:2756953/71‑1 (MQ=255)
                                                      cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACACGGTg                 <  1:4227687/66‑1 (MQ=255)
                                                        ggCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                 <  1:3488427/65‑1 (MQ=255)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg           <  1:4102944/70‑1 (MQ=255)
                                                         gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg           <  1:6348/70‑1 (MQ=255)
                                                          cGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcggc         <  1:4521222/71‑1 (MQ=255)
                                                            gAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                                >  1:2255971/1‑44 (MQ=255)
                                                            gAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCaa       >  1:2592744/1‑71 (MQ=255)
                                                             aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                                    <  1:1047800/41‑1 (MQ=38)
                                                             aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt      <  1:354690/71‑1 (MQ=255)
                                                               catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAaaga                                  <  1:3790835/41‑1 (MQ=38)
                                                                aTCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtg   <  1:580753/71‑1 (MQ=255)
                                                                 tCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                                      <  1:842187/35‑1 (MQ=37)
                                                                 tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtgg  >  1:1549055/1‑71 (MQ=255)
                                                                 tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtgg  <  1:3171544/71‑1 (MQ=255)
                                                                    |                                                                   
CTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGG  >  AP009048/3824845‑3824980

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: