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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,2130TC100.0% 303.7 / NA 87A189A (GCT→GCCpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/51);  total (36/51)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGTTTG  >  AP009048/3825147‑3825280
                                                                  |                                                                   
tCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTatc                                                                 >  1:1380045/1‑69 (MQ=255)
 cAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCt                                                                >  1:723800/1‑71 (MQ=255)
  aGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCt                                                                >  1:225641/1‑70 (MQ=255)
   ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                              <  1:5199028/71‑1 (MQ=255)
   ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACACTGAAAGATCTGATTGCCTATCTg                                                               >  1:3569649/1‑70 (MQ=255)
     aaGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTatc                                                                 >  1:3972867/1‑64 (MQ=255)
     aaGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                              >  1:1136426/1‑69 (MQ=255)
     aaGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                              >  1:4552958/1‑69 (MQ=255)
     aaGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGaa                                                            <  1:1989433/71‑1 (MQ=255)
     aaGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGCaa                                                            <  1:530275/71‑1 (MQ=255)
      aGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                              >  1:3746999/1‑68 (MQ=255)
      aGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                              >  1:2727574/1‑68 (MQ=255)
        ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                         <  1:656157/71‑1 (MQ=255)
        ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                         <  1:1824974/71‑1 (MQ=255)
        ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                         <  1:4766633/71‑1 (MQ=255)
        ttGAGCGTGATTACAACTGCAAAGTGATCTCTATAATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                         <  1:207557/71‑1 (MQ=255)
           aGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                         <  1:4908944/68‑1 (MQ=255)
           aGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGaa                                                       <  1:3992676/70‑1 (MQ=255)
           aGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGaa                                                       <  1:2769300/70‑1 (MQ=255)
           aGCGTCATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTTGAagag                                                         <  1:2199811/68‑1 (MQ=255)
             cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                    >  1:3836769/1‑71 (MQ=255)
             cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                    <  1:1583298/71‑1 (MQ=255)
             cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                    <  1:3579771/71‑1 (MQ=255)
             cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                    >  1:2159176/1‑71 (MQ=255)
               tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAATGATCTGATTGCCTATCTGGAaga                                                          <  1:4077854/63‑1 (MQ=255)
               tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTat                                                                  >  1:3024279/1‑54 (MQ=255)
               tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                              <  1:2849197/59‑1 (MQ=255)
               tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                          >  1:1487725/1‑63 (MQ=255)
                gATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                         <  1:4516012/63‑1 (MQ=255)
                 aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                                 <  1:5199474/70‑1 (MQ=255)
                 aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                                 >  1:4544600/1‑70 (MQ=255)
                 aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                                 <  1:1628120/70‑1 (MQ=255)
                  ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGaa                                                            <  1:253477/58‑1 (MQ=38)
                  ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                               <  1:2217108/71‑1 (MQ=255)
                  ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                               >  1:3635434/1‑71 (MQ=255)
                  ttACAACTGCAAAGTGATCTCTATCATCACACTGAAAGATCTGATTGCCTATCTGGaa                                                            <  1:4911055/58‑1 (MQ=38)
                    aCATCTGCAATGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGTAGCCGGAAAt                                              <  1:1480923/70‑1 (MQ=255)
                    aCAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTg                                                               <  1:123844/53‑1 (MQ=255)
                    aCAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                             >  1:1658584/1‑71 (MQ=255)
                     cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                         <  1:1583524/58‑1 (MQ=38)
                     cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                             <  1:4913651/70‑1 (MQ=255)
                     cAACTGCAAAGTATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                            <  1:4058534/70‑1 (MQ=255)
                      aaCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGATAAGCCGGAAATGGc                                           <  1:2026270/71‑1 (MQ=255)
                        cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                         <  1:233654/71‑1 (MQ=255)
                        cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                         <  1:2966120/71‑1 (MQ=255)
                        cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                          >  1:4700676/1‑70 (MQ=255)
                          gCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGc                                           <  1:573974/67‑1 (MQ=255)
                           cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGc                                           <  1:749294/66‑1 (MQ=255)
                           cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGa                                        <  1:1551747/69‑1 (MQ=255)
                            aaaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                             <  1:2053783/63‑1 (MQ=255)
                            aaaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACa                                     <  1:3771984/71‑1 (MQ=255)
                             aaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                          >  1:1761266/1‑49 (MQ=37)
                             aaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACAt                                    <  1:2255514/71‑1 (MQ=255)
                                tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                               >  1:145823/1‑57 (MQ=38)
                                tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTg                                 <  1:3027874/71‑1 (MQ=255)
                                  aTCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                               >  1:517287/1‑55 (MQ=38)
                                  aTCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                                >  1:523616/1‑54 (MQ=38)
                                     tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                                      >  1:4099755/1‑45 (MQ=37)
                                     tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                              >  1:2108268/1‑53 (MQ=38)
                                      ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGgcggcg                           >  1:70968/1‑71 (MQ=255)
                                      ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGgcggcg                           >  1:64746/1‑71 (MQ=255)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcgg                             <  1:2623498/68‑1 (MQ=255)
                                       tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGcggcgg                          >  1:2702396/1‑71 (MQ=255)
                                        atcatcACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                          <  1:2836143/54‑1 (MQ=38)
                                         tcatcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                        <  1:2710251/71‑1 (MQ=255)
                                         tcatcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                        <  1:3833124/71‑1 (MQ=255)
                                         tcatcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                        <  1:1115469/71‑1 (MQ=255)
                                         tcatcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                        <  1:1797598/71‑1 (MQ=255)
                                           atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                        >  1:4324868/1‑69 (MQ=39)
                                           atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa                      <  1:1685945/71‑1 (MQ=255)
                                           atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa                      <  1:1045871/71‑1 (MQ=255)
                                           atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTa                       <  1:4748167/70‑1 (MQ=255)
                                           atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTa                       <  1:3237067/70‑1 (MQ=255)
                                            tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa                      <  1:3568076/70‑1 (MQ=255)
                                              aCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGc                   >  1:2270666/1‑71 (MQ=255)
                                                ccTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                         <  1:2969039/47‑1 (MQ=37)
                                                ccTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGcc                  >  1:1203659/1‑70 (MQ=255)
                                                 cTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTa                <  1:437362/71‑1 (MQ=255)
                                                  tGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTAt               >  1:1469702/1‑71 (MQ=255)
                                                    aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                              >  1:1096724/1‑38 (MQ=21)
                                                    aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcg             >  1:1116052/1‑71 (MQ=255)
                                                      aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                         <  1:4251447/41‑1 (MQ=25)
                                                         tcTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg           >  1:2789976/2‑68 (MQ=39)
                                                             aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGtt    >  1:4900800/1‑71 (MQ=255)
                                                             aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGtt    >  1:2553411/1‑71 (MQ=255)
                                                               tGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                        <  1:156811/49‑1 (MQ=38)
                                                               tGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGTTTg  >  1:3273487/1‑71 (MQ=255)
                                                                  |                                                                   
TCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGTTTG  >  AP009048/3825147‑3825280

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: