Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,216 C→T 100% Y190Y (TAC→TAT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,2160CT100.0% 249.5 / NA 84Y190Y (TAC→TATpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/51);  total (33/51)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGTTTG  >  AP009048/3825148‑3825280
                                                                    |                                                                
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)
    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:2727574/1‑68 (MQ=255)
     aGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                              >  1:3746999/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:2769300/70‑1 (MQ=255)
          aGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGaa                                                       <  1:3992676/70‑1 (MQ=255)
          aGCGTCATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTTGAagag                                                         <  1:2199811/68‑1 (MQ=255)
            cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                    <  1:3579771/71‑1 (MQ=255)
            cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                    <  1:1583298/71‑1 (MQ=255)
            cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                    >  1:3836769/1‑71 (MQ=255)
            cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                                    >  1:2159176/1‑71 (MQ=255)
              tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAATGATCTGATTGCCTATCTGGAaga                                                          <  1:4077854/63‑1 (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:1628120/70‑1 (MQ=255)
                aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                                 <  1:5199474/70‑1 (MQ=255)
                aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                                 >  1:4544600/1‑70 (MQ=255)
                 ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGaa                                                            <  1:253477/58‑1 (MQ=38)
                 ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                               >  1:3635434/1‑71 (MQ=255)
                 ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                               <  1:2217108/71‑1 (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:64746/1‑71 (MQ=255)
                                     ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGgcggcg                           >  1:70968/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)
                                                                    |                                                                
CAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAGTTTG  >  AP009048/3825148‑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: