Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,204 C→T 100% D186D (GAC→GAT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,2040CT100.0% 201.0 / NA 67D186D (GAC→GATpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/35);  total (32/35)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCC  >  AP009048/3825136‑3825264
                                                                    |                                                            
aTTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCt                                                            >  1:2291189/1‑71 (MQ=255)
 tttCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTg                                                           <  1:4262951/71‑1 (MQ=255)
  ttCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGa                                                          >  1:4415318/1‑71 (MQ=255)
   tCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGa                                                          <  1:606085/70‑1 (MQ=255)
     ggCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAtt                                                        <  1:2975268/70‑1 (MQ=255)
        gATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGcc                                                     >  1:32415/1‑69 (MQ=255)
        gATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGcc                                                     >  1:3930705/1‑69 (MQ=255)
        gATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTATATCATCACCCTGAAAGATCTGATTGcc                                                     >  1:4047861/1‑69 (MQ=255)
         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:3682894/71‑1 (MQ=255)
                   ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                         <  1:437638/71‑1 (MQ=255)
                    tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                         >  1:793342/1‑54 (MQ=255)
                    tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                         >  1:4369589/1‑70 (MQ=255)
                    tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                         <  1:645461/70‑1 (MQ=255)
                     gAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                         >  1:2769709/1‑53 (MQ=255)
                     gAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                         >  1:4090439/1‑69 (MQ=255)
                      aGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                      <  1:1848380/71‑1 (MQ=255)
                       gCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                         >  1:2921963/1‑51 (MQ=255)
                       gCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                     <  1:3587294/71‑1 (MQ=255)
                           gATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAag                                           <  1:2919795/61‑1 (MQ=255)
                            aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                 <  1:3154490/70‑1 (MQ=255)
                            aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                 <  1:1243729/70‑1 (MQ=255)
                              tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                         >  1:3774777/1‑44 (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)
                                   cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTg                                                           <  1:1942916/37‑1 (MQ=255)
                                   cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTg                                                           <  1:1327946/37‑1 (MQ=255)
                                   cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                         <  1:892133/71‑1 (MQ=255)
                                   cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                          >  1:3141773/1‑70 (MQ=255)
                                   cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                          >  1:4268880/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:3254212/1‑63 (MQ=255)
                                                 ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATc                   >  1:3235700/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)
                                                                    |                                                            
ATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCC  >  AP009048/3825136‑3825264

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: