Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,108 C→T 100% L154L (CTC→CTT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,1080CT98.7% 224.0 / ‑4.3 76L154L (CTC→CTTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/44);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
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.

CGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCA  >  AP009048/3825041‑3825175
                                                                   |                                                                   
cGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGAt                                                                  >  1:4939862/1‑71 (MQ=255)
 gATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                                 <  1:4188034/71‑1 (MQ=255)
   ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                               <  1:1307027/71‑1 (MQ=255)
   ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                               <  1:4841797/71‑1 (MQ=255)
    tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                               <  1:3794901/70‑1 (MQ=255)
    tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                                                               <  1:2974398/70‑1 (MQ=255)
     cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTc                                                              <  1:557728/70‑1 (MQ=255)
       cgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                                                             >  1:2111461/1‑69 (MQ=255)
       cgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                                           <  1:4351201/71‑1 (MQ=255)
        gaatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGa                                                          <  1:2676720/68‑1 (MQ=255)
         aatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGaa                                                         >  1:486598/3‑71 (MQ=255)
          atCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                                                        <  1:3211332/70‑1 (MQ=255)
           tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                                <  1:1901699/62‑1 (MQ=255)
           tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                                                             >  1:49131/1‑65 (MQ=255)
           tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                                                             >  1:4035541/1‑65 (MQ=255)
           tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcg                                                       <  1:4945953/71‑1 (MQ=255)
             gATGGAGATTATTCAGGCTAATGGCGGGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                                                     <  1:4502927/71‑3 (MQ=255)
             gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcggg                                                     >  1:1258178/1‑69 (MQ=255)
             gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                                                      <  1:2912654/70‑2 (MQ=255)
             gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                                                      <  1:4375421/70‑2 (MQ=255)
              aTGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                                                           <  1:2598413/64‑1 (MQ=255)
               tGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGgggg                                                   <  1:4910219/71‑1 (MQ=255)
                ggAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc                                                  >  1:1833793/1‑71 (MQ=255)
                ggAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc                                                  >  1:2954676/1‑71 (MQ=255)
                 gagaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                                                                 >  1:2046396/1‑55 (MQ=255)
                  agaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcgc                                                >  1:5158911/1‑71 (MQ=255)
                     ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc                                             >  1:4599990/1‑71 (MQ=255)
                      tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                                                        <  1:1408462/59‑1 (MQ=255)
                       attCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTc                                                              >  1:2909494/1‑52 (MQ=255)
                        ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                                                             >  1:1824494/1‑52 (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:1146799/3‑71 (MQ=255)
                              ctaaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCg                                    <  1:1901090/69‑1 (MQ=255)
                                aaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTc                                     <  1:1762238/68‑1 (MQ=255)
                                aaTGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCgg                                   >  1:969761/1‑70 (MQ=255)
                                  tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                                <  1:3139225/39‑1 (MQ=37)
                                  tGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcg                                                       <  1:938154/48‑1 (MQ=39)
                                   ggCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                                                                <  1:4367227/38‑1 (MQ=38)
                                     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:4261465/70‑1 (MQ=255)
                                                          gATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTAc         <  1:663771/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)
                                                                   |                                                                   
CGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCA  >  AP009048/3825041‑3825175

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