Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,123 C→G 100% R159R (CGC→CGG pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,1230CG100.0% 201.2 / NA 58R159R (CGC→CGGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (29/29);  total (29/29)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

TGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCA  >  AP009048/3825056‑3825190
                                                                   |                                                                   
tGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGgggg                                                                  <  1:4910219/71‑1 (MQ=255)
 ggAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc                                                                 >  1:1833793/1‑71 (MQ=255)
 ggAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc                                                                 >  1:2954676/1‑71 (MQ=255)
   agaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcgc                                                               >  1:5158911/1‑71 (MQ=255)
      ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                                            >  1:4599990/1‑71 (MQ=255)
           cAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAt                                                       >  1:4052875/1‑71 (MQ=255)
            aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAtt                                                      <  1:592841/71‑1 (MQ=255)
             ggCTAATGGCCCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                                     <  1:2476450/71‑1 (MQ=255)
              gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                                     <  1:2021689/67‑1 (MQ=255)
               ctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCg                                                   <  1:1901090/69‑1 (MQ=255)
               ctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCg                                                   >  1:1146799/3‑71 (MQ=255)
                 aaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                                    <  1:1762238/68‑1 (MQ=255)
                 aaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCgg                                                  >  1:969761/1‑70 (MQ=255)
                      cgcgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                                     <  1:2355370/62‑1 (MQ=255)
                       gcgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                                             >  1:2399176/1‑69 (MQ=255)
                       gcgACGCTTGCTGGCGTATTGATTTCGCTTGATCCTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                                             >  1:80286/1‑69 (MQ=255)
                        cgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAg                                          <  1:3163753/71‑1 (MQ=255)
                         gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg                                         <  1:4849352/71‑1 (MQ=255)
                         gACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAg                                          >  1:353631/1‑70 (MQ=255)
                           cGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg                                         >  1:4973471/1‑69 (MQ=255)
                           cGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGa                                        <  1:2896662/70‑1 (MQ=255)
                                gCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt                                    >  1:1020674/1‑69 (MQ=255)
                                gCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTg                                   <  1:5123487/70‑1 (MQ=255)
                                 cTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa                                  <  1:1316577/70‑1 (MQ=255)
                                  tGGCGTATTGATTTCGCTTGATCGTCAGGAACGgggg                                                                  <  1:164915/37‑1 (MQ=25)
                                    gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCTAGATTTCGGCGATTCAGGAAGTTGAGCGt                              <  1:4206122/71‑1 (MQ=255)
                                    gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCg                               <  1:1903097/70‑1 (MQ=255)
                                      gtattGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGa                            >  1:1583564/4‑71 (MQ=255)
                                      gtattGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGa                            >  1:3820412/4‑71 (MQ=255)
                                       tattGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGt                              <  1:3064294/66‑1 (MQ=255)
                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCgaga                                                        >  1:2607080/1‑40 (MQ=38)
                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa                                  >  1:328850/1‑62 (MQ=255)
                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGAtt                          <  1:4008152/70‑1 (MQ=255)
                                           gATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAg                                 >  1:3957880/1‑61 (MQ=255)
                                           gATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTAc                        <  1:663771/70‑1 (MQ=255)
                                           gATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTAc                        <  1:4261465/70‑1 (MQ=255)
                                           gATTTCGCTCGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTAc                        <  1:3501724/70‑1 (MQ=255)
                                            aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACaa                      <  1:3892138/71‑1 (MQ=255)
                                            aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACaa                      <  1:2118813/71‑1 (MQ=255)
                                             tttCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGt                              <  1:994333/62‑1 (MQ=255)
                                               tCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACt                    <  1:1840356/70‑1 (MQ=255)
                                               tCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTg                   >  1:4396959/1‑71 (MQ=255)
                                                cGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACt                    <  1:3234063/69‑1 (MQ=255)
                                                 gCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                                    >  1:444682/1‑36 (MQ=37)
                                                 gCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                                    >  1:3871366/1‑36 (MQ=37)
                                                 gCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGc                  >  1:4302177/1‑70 (MQ=255)
                                                 gCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGc                  >  1:1466578/1‑70 (MQ=255)
                                                 gCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCa                 <  1:1327931/71‑1 (MQ=255)
                                                  cttgATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCaa                <  1:1414900/68‑1 (MQ=255)
                                                     gATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAAc                     >  1:1720376/1‑63 (MQ=255)
                                                     gATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAg              <  1:4107919/70‑1 (MQ=255)
                                                       tCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCa                 <  1:3089145/65‑1 (MQ=255)
                                                          tCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTg            <  1:706693/67‑1 (MQ=255)
                                                          tCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGAtct        >  1:4259625/1‑71 (MQ=255)
                                                           cAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGAt          >  1:5028077/1‑68 (MQ=255)
                                                           cAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGAt          >  1:1924738/1‑68 (MQ=255)
                                                             ggAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTa     >  1:1688660/1‑71 (MQ=255)
                                                                aCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTatca  <  1:574886/71‑1 (MQ=255)
                                                                aCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTatca  >  1:1479902/1‑71 (MQ=255)
                                                                aCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTatca  >  1:2469118/1‑71 (MQ=255)
                                                                   |                                                                   
TGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCA  >  AP009048/3825056‑3825190

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