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,1230CG99.0% 249.0 / ‑6.0 72R159R (CGC→CGGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (27/44);  minor base T (1/0);  total (28/44)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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.

TGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTAT  >  AP009048/3825056‑3825188
                                                                   |                                                                 
tgcagaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGgggg                                                                <  1:6495493/68‑1 (MQ=255)
  gagaTTATTCAGGCTAGTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcg                                                              <  1:7508201/71‑1 (MQ=255)
   agaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcgc                                                             <  1:3476906/71‑1 (MQ=255)
      ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                                          >  1:357762/1‑71 (MQ=255)
        attCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggcg                                                         <  1:2487428/70‑1 (MQ=255)
         ttgaGGCTATTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCga                                                        <  1:6790716/67‑1 (MQ=255)
         ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCga                                                        <  1:6430320/70‑1 (MQ=255)
          tCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                                          >  1:5225201/1‑67 (MQ=255)
           cAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAt                                                     >  1:4139146/1‑71 (MQ=255)
            aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAtt                                                    <  1:3833071/71‑1 (MQ=255)
             ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt                                                   <  1:7350767/71‑1 (MQ=255)
              gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                                  <  1:4689352/68‑1 (MQ=255)
              gctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                                  <  1:2315815/68‑1 (MQ=255)
                 aaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                                               >  1:3837055/1‑71 (MQ=255)
                  aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCgg                                                <  1:2229671/69‑1 (MQ=255)
                  aTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCgg                                                <  1:2024813/69‑1 (MQ=255)
                   tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGa                                             <  1:800517/71‑1 (MQ=255)
                     gcgcGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                                           <  1:7416865/71‑1 (MQ=255)
                     gcgcGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                                           <  1:7566626/71‑1 (MQ=255)
                      cgcgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTc                                          <  1:6785786/71‑1 (MQ=255)
                      cgcgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGCCGATTc                                          <  1:4525007/71‑1 (MQ=255)
                       gcgACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCa                                         <  1:7031708/71‑1 (MQ=255)
                          aCGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGgggg                                                                <  1:3725865/45‑1 (MQ=38)
                            gCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGTGGgcg                                                              >  1:918467/1‑45 (MQ=38)
                            gCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAg                                    >  1:848342/1‑71 (MQ=255)
                             cTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                                               <  1:5283355/59‑1 (MQ=255)
                             cTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                                               <  1:3750757/59‑1 (MQ=255)
                             cTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGc                                               <  1:5265666/59‑1 (MQ=255)
                                gCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt                                  <  1:7671517/69‑1 (MQ=255)
                                 cTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa                                <  1:2018133/70‑1 (MQ=255)
                                   ggCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCg                             >  1:7003781/1‑71 (MQ=255)
                                   ggCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCg                             >  1:506524/1‑71 (MQ=255)
                                    gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGtt                                  <  1:1610088/65‑1 (MQ=255)
                                    gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGc                              <  1:3917853/69‑1 (MQ=255)
                                    gCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGt                            <  1:6366694/71‑1 (MQ=255)
                                      gtattGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGa                          >  1:466300/4‑71 (MQ=255)
                                       tattGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                                  <  1:3346246/44‑1 (MQ=38)
                                        attGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGAt                         >  1:3564878/2‑70 (MQ=255)
                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGAtt                        <  1:4341040/70‑1 (MQ=255)
                                         ttGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTa                       <  1:6153803/71‑1 (MQ=255)
                                          tGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTa                       >  1:3099649/1‑70 (MQ=255)
                                            aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACaa                    <  1:1214513/71‑1 (MQ=255)
                                            aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACaa                    >  1:442962/1‑71 (MQ=255)
                                            aTTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACa                     <  1:6057039/70‑1 (MQ=255)
                                             tttCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAgg                                       <  1:7902882/51‑1 (MQ=255)
                                               tCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGAtt                                           >  1:4306547/1‑45 (MQ=255)
                                                cGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTg                 >  1:2452764/1‑70 (MQ=255)
                                                   ttgATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCa               >  1:7226851/3‑69 (MQ=255)
                                                    tgATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGaa                                     >  1:7514079/2‑46 (MQ=255)
                                                    tgATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCaa              <  1:1831210/68‑1 (MQ=255)
                                                    tgATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAg            >  1:669132/2‑71 (MQ=255)
                                                     gATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAg                               >  1:7634300/1‑51 (MQ=255)
                                                     gATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGa                          >  1:4041512/1‑56 (MQ=255)
                                                     gATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGt           >  1:6755176/1‑71 (MQ=255)
                                                     gATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGt           >  1:1618125/1‑71 (MQ=255)
                                                     gATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGt           >  1:1627321/1‑71 (MQ=255)
                                                      aTCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTg          <  1:435955/71‑1 (MQ=255)
                                                        cGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGAt        >  1:3531612/1‑71 (MQ=255)
                                                          tCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa                                <  1:3495287/45‑1 (MQ=255)
                                                           cAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGa                                >  1:769484/1‑44 (MQ=255)
                                                           cAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAg                               >  1:5448236/1‑45 (MQ=255)
                                                           cAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAg                               >  1:2260655/1‑45 (MQ=255)
                                                            aGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACt                  >  1:3576937/1‑57 (MQ=255)
                                                            aGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACt                  >  1:2306142/1‑57 (MQ=255)
                                                            aGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATctct    <  1:7169939/71‑1 (MQ=255)
                                                             ggAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTa   <  1:6830450/71‑1 (MQ=255)
                                                              gAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGa                                      <  1:2212033/35‑1 (MQ=255)
                                                              gAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGa                                      <  1:3225603/35‑1 (MQ=255)
                                                              gAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGa                                      <  1:5129365/35‑1 (MQ=255)
                                                              gAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGAtc       <  1:7916849/66‑1 (MQ=255)
                                                              gAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTat  <  1:1462365/71‑1 (MQ=255)
                                                              gAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTat  <  1:1079423/71‑1 (MQ=255)
                                                               aaCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGt                            <  1:7561623/44‑1 (MQ=255)
                                                                aCGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTa   <  1:5970889/68‑1 (MQ=255)
                                                                   |                                                                 
TGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTAT  >  AP009048/3825056‑3825188

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