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% 177.2 / NA 60Y190Y (TAC→TATpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/33);  total (27/33)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAG  >  AP009048/3825149‑3825276
                                                                   |                                                            
aGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTg                                                           <  1:156869/71‑1 (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)
                                                                   |                                                            
AGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCGAAGAG  >  AP009048/3825149‑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: