Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,919 G→T 100% L91L (CTG→CTT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,9190GT100.0% 155.9 / NA 53L91L (CTG→CTTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (23/30);  total (23/30)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGG  >  AP009048/3824855‑3824984
                                                                |                                                                 
cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                             >  1:4154892/1‑71 (MQ=255)
   cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                          >  1:4034216/1‑71 (MQ=255)
   cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                          <  1:3952206/71‑1 (MQ=255)
   cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                          <  1:2477700/71‑1 (MQ=255)
      aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                                                 >  1:4554861/1‑61 (MQ=255)
          aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                   >  1:1022278/1‑71 (MQ=255)
          aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                   >  1:2749567/1‑71 (MQ=255)
             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)
                   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)
                           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:1464760/4‑70 (MQ=255)
                                         gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                     <  1:4145690/67‑1 (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:6348/70‑1 (MQ=255)
                                               gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg               <  1:4102944/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)
                                                         atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCaaaa                               >  1:1827565/3‑44 (MQ=38)
                                                         atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGt                      >  1:779067/3‑53 (MQ=255)
                                                         atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt    >  1:180533/3‑71 (MQ=255)
                                                         atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt    >  1:3935334/3‑71 (MQ=255)
                                                           gaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCa                                  <  1:4384097/39‑1 (MQ=38)
                                                            accttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCaa           <  1:4358095/61‑1 (MQ=255)
                                                            accttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttgg  <  1:4422715/70‑1 (MQ=255)
                                                            accttgacctgCCGTAATGCTTTAACCGCAAAGAGGCAAAAGACCAc                         >  1:419034/1‑47 (MQ=38)
                                                                |                                                                 
CTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGG  >  AP009048/3824855‑3824984

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: