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% 215.3 / NA 63A189A (GCT→GCCpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/35);  total (28/35)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAG  >  AP009048/3825145‑3825276
                                                                    |                                                               
aTTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTa                                                               >  1:1231590/1‑71 (MQ=255)
    aGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTg                                                           <  1:156869/71‑1 (MQ=255)
     ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTatc                                                             <  1:1919417/68‑3 (MQ=255)
     ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTat                                                              <  1:3574004/67‑2 (MQ=255)
     ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTg                                                           >  1:1148170/1‑70 (MQ=255)
      gAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                          >  1:4062535/1‑70 (MQ=255)
        aGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAag                                                       <  1:2579895/71‑1 (MQ=255)
         gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAATATCTGATTGCCTATCTGGAaga                                                      <  1:3929/71‑1 (MQ=255)
          ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                     <  1:437638/71‑1 (MQ=255)
          ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                     <  1:3682894/71‑1 (MQ=255)
           tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                     <  1:645461/70‑1 (MQ=255)
           tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                     >  1:4369589/1‑70 (MQ=255)
            gAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                     >  1:4090439/1‑69 (MQ=255)
             aGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                                  <  1:1848380/71‑1 (MQ=255)
              gCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                                 <  1:3587294/71‑1 (MQ=255)
                  gATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAag                                                       <  1:2919795/61‑1 (MQ=255)
                   aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                             <  1:1243729/70‑1 (MQ=255)
                   aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                             <  1:3154490/70‑1 (MQ=255)
                     tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                           >  1:4256143/1‑70 (MQ=255)
                     tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                          >  1:3122057/1‑71 (MQ=255)
                       cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                        <  1:1157502/71‑1 (MQ=255)
                        aaCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGc                                       <  1:1574414/71‑1 (MQ=255)
                          cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                     <  1:892133/71‑1 (MQ=255)
                          cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                      >  1:4268880/1‑70 (MQ=255)
                          cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                      >  1:3141773/1‑70 (MQ=255)
                             cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGaa                                   <  1:1534286/70‑1 (MQ=255)
                             cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGa                                    <  1:3117652/69‑1 (MQ=255)
                               aaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                      >  1:4512178/1‑49 (MQ=37)
                                 gTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCt                              <  1:509692/71‑1 (MQ=255)
                                  tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                        <  1:2658754/60‑1 (MQ=255)
                                    aTCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGAGGAACATCTgg                            >  1:2263009/1‑70 (MQ=255)
                                      ctctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcgg                         >  1:3141809/1‑71 (MQ=255)
                                       tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGc                                       <  1:3654010/56‑1 (MQ=38)
                                       tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                     <  1:2098638/58‑1 (MQ=38)
                                        ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                      <  1:3547468/40‑1 (MQ=25)
                                        ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATc                               >  1:3235700/1‑63 (MQ=255)
                                        ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATc                               >  1:3254212/1‑63 (MQ=255)
                                        ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCt                              >  1:3900517/1‑64 (MQ=255)
                                        ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTg                             >  1:408241/1‑65 (MQ=255)
                                           tcatcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                    <  1:2691774/71‑1 (MQ=255)
                                           tcatcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                    <  1:2049608/71‑1 (MQ=255)
                                            catcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                    >  1:2346770/1‑70 (MQ=255)
                                            catcaACCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                    >  1:2658198/1‑70 (MQ=255)
                                             atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGAGGCGGTTa                   <  1:1229661/70‑1 (MQ=255)
                                              tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                           <  1:808548/45‑1 (MQ=37)
                                              tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                     <  1:1457950/51‑1 (MQ=38)
                                              tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAg                 >  1:115142/1‑71 (MQ=255)
                                                aCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcg                          >  1:1805099/1‑60 (MQ=255)
                                                 cccTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGcc              >  1:1544886/1‑71 (MQ=255)
                                                   cTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAg                 <  1:831298/66‑1 (MQ=39)
                                                     gAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                             <  1:3767437/36‑1 (MQ=21)
                                                      aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                           >  1:4346466/1‑37 (MQ=21)
                                                        aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa                  >  1:3013971/1‑60 (MQ=38)
                                                         gatcTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg       <  1:1406325/67‑1 (MQ=255)
                                                          atcTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                    <  1:825021/54‑1 (MQ=38)
                                                            cTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCt                              <  1:3054520/44‑1 (MQ=38)
                                                             tGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACAt                                >  1:3984845/1‑41 (MQ=38)
                                                             tGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt                    <  1:2935577/53‑1 (MQ=38)
                                                              gATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGaa                                   >  1:2602652/1‑37 (MQ=37)
                                                              gATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGaa                                   >  1:488464/1‑37 (MQ=37)
                                                              gATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGaa                                   >  1:842866/1‑37 (MQ=37)
                                                               aTTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGAGGCGGTTAAGGCCTATCTCGAagag  <  1:22680/69‑1 (MQ=255)
                                                                 tGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTgg                            >  1:2449540/1‑41 (MQ=37)
                                                                    |                                                               
ATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAG  >  AP009048/3825145‑3825276

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