Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,176,730 T→C 100% Y349C (TAC→TGC)  uxaC ← glucuronate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,176,7300TC99.1% 425.4 / ‑4.8 112Y349C (TAC→TGC) uxaCglucuronate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (52/59);  total (52/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-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.

Only 100 of 120 total aligned reads displayed.

TTACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGTAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAGC  >  REL606/3176682‑3176777
                                                |                                               
ttACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGcagcag                                               <  1:2513683/51‑1 (MQ=255)
ttACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGcagcag                                               >  1:5597554/1‑51 (MQ=255)
ttACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGcagcag                                               <  1:1619420/51‑1 (MQ=255)
 tACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCaga                                              <  1:5095983/51‑1 (MQ=255)
 tACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCaga                                              <  1:10627767/51‑1 (MQ=255)
 tACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCaga                                              <  1:2637289/51‑1 (MQ=255)
  aCCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCagag                                             <  1:6428291/51‑1 (MQ=255)
   ccGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAgaga                                            <  1:2914545/51‑1 (MQ=255)
    cGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGaa                                           <  1:10747740/51‑1 (MQ=255)
    cGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGaa                                           <  1:3169263/51‑1 (MQ=255)
    cGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGaa                                           <  1:7447507/51‑1 (MQ=255)
     gATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAAt                                          <  1:1765262/51‑1 (MQ=255)
     gATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAAt                                          >  1:3980043/1‑51 (MQ=255)
     gATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAAt                                          <  1:5166842/51‑1 (MQ=255)
     gATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAAt                                          <  1:6452193/51‑1 (MQ=255)
      aTCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATg                                         >  1:10366350/1‑51 (MQ=255)
      aTCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATg                                         <  1:9325173/51‑1 (MQ=255)
       tCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATgg                                        >  1:2222926/1‑51 (MQ=255)
       tCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATgg                                        <  1:3389183/51‑1 (MQ=255)
        cATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGt                                       >  1:1171104/1‑51 (MQ=255)
         aTGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtt                                      >  1:3892888/1‑51 (MQ=255)
         aTGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtt                                      <  1:6815832/51‑1 (MQ=255)
         aTGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtt                                      <  1:5993719/51‑1 (MQ=255)
          tGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGttt                                     >  1:10190399/1‑51 (MQ=255)
          tGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGttt                                     <  1:4500008/51‑1 (MQ=255)
           ggTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtttt                                    >  1:2102316/1‑51 (MQ=255)
           ggTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGtttt                                    <  1:66350/51‑1 (MQ=255)
            gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   <  1:3229702/51‑1 (MQ=255)
            gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   >  1:3483954/1‑51 (MQ=255)
            gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   <  1:5284360/51‑1 (MQ=255)
            gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   >  1:1412901/1‑51 (MQ=255)
            gTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTc                                   <  1:6722022/51‑1 (MQ=255)
              gCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCgg                                 <  1:8438161/51‑1 (MQ=255)
               cccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgc                                >  1:6653531/1‑51 (MQ=255)
               cccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgc                                >  1:2361622/1‑51 (MQ=255)
               cccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgc                                >  1:7986450/1‑51 (MQ=255)
               cccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgc                                >  1:1309288/1‑51 (MQ=255)
                ccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCTGAGAATGGTTTTCGgca                               <  1:5545965/51‑1 (MQ=38)
                ccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgca                               >  1:6668000/1‑51 (MQ=255)
                ccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgca                               >  1:4232747/1‑51 (MQ=255)
                ccAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgca                               >  1:10473553/1‑51 (MQ=255)
                 cAGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgcag                              >  1:10377504/1‑51 (MQ=255)
                  aGCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgcagc                             >  1:7550984/1‑51 (MQ=255)
                   gCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgcagca                            <  1:4745076/51‑1 (MQ=255)
                   gCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgcagca                            <  1:9172133/51‑1 (MQ=255)
                   gCACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGgca                               >  1:5007329/1‑48 (MQ=255)
                    cACTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGcagcag                           >  1:3188497/1‑51 (MQ=255)
                     aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          >  1:9467692/1‑51 (MQ=255)
                     aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          <  1:6474601/51‑1 (MQ=255)
                     aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          <  1:5448173/51‑1 (MQ=255)
                     aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          >  1:4408224/1‑51 (MQ=255)
                     aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          <  1:2764111/51‑1 (MQ=255)
                     aCTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          >  1:1042065/1‑51 (MQ=255)
                      cTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGt                         <  1:8762335/51‑1 (MQ=255)
                      cTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGt                         >  1:7329677/1‑51 (MQ=255)
                      cTTCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGt                         >  1:2244285/1‑51 (MQ=255)
                       ttCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          >  1:5339435/1‑49 (MQ=255)
                       ttCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGtt                        <  1:2257041/51‑1 (MQ=255)
                        tCGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGcagcag                           <  1:582316/47‑1 (MQ=255)
                         cGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAgg                          >  1:2566394/1‑47 (MQ=255)
                         cGTTATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGtttt                      <  1:9347112/51‑1 (MQ=255)
                           ttATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTttct                    >  1:4255353/1‑51 (MQ=255)
                            tATCGCGCGGGTTCAGGCAGCGGAGAATGGTTTTCGGCAGCAGGTTttctt                   >  1:8274689/1‑51 (MQ=255)
                            tATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTttctt                   >  1:10223250/1‑51 (MQ=255)
                            tATCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTttctt                   >  1:6323834/1‑51 (MQ=255)
                              tCGCGCGGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCa                 <  1:10844456/51‑1 (MQ=255)
                               cgcgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTttct                    <  1:5197235/47‑1 (MQ=255)
                               cgcgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCAt                >  1:8083999/1‑51 (MQ=255)
                               cgcgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCAt                <  1:7873785/51‑1 (MQ=255)
                               cgcgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCAt                <  1:3219023/51‑1 (MQ=255)
                                gcgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCAtt               >  1:3705320/1‑51 (MQ=255)
                                 cgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTc              >  1:3286005/1‑51 (MQ=255)
                                 cgcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTc              <  1:2920526/51‑1 (MQ=255)
                                  gcgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCt             <  1:1164899/51‑1 (MQ=255)
                                   cgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtg            >  1:3892246/1‑51 (MQ=255)
                                   cgGGTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtg            >  1:10914213/1‑51 (MQ=255)
                                    gggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgc           <  1:3035602/51‑1 (MQ=255)
                                    gggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgc           >  1:391559/1‑51 (MQ=255)
                                    gggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgc           >  1:4550342/1‑51 (MQ=255)
                                    gggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgc           <  1:3592645/51‑1 (MQ=255)
                                     ggTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgct          >  1:1664651/1‑51 (MQ=255)
                                      gTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgctt         <  1:9364744/51‑1 (MQ=255)
                                      gTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgctt         <  1:101007/51‑1 (MQ=255)
                                      gTTCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgctt         <  1:5712666/51‑1 (MQ=255)
                                       ttCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttg        <  1:153430/51‑1 (MQ=255)
                                       ttCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttg        <  1:8027009/51‑1 (MQ=255)
                                       ttCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttg        >  1:10765616/1‑51 (MQ=255)
                                        tCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgc       >  1:3324381/1‑51 (MQ=255)
                                        tCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgc       <  1:6388064/51‑1 (MQ=255)
                                        tCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgc       >  1:404963/1‑51 (MQ=255)
                                        tCAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgc       <  1:2038108/51‑1 (MQ=255)
                                         cAGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCtgcttgct      >  1:7281051/1‑51 (MQ=255)
                                          aGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTc     <  1:234703/51‑1 (MQ=255)
                                          aGGCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTc     <  1:7758845/51‑1 (MQ=255)
                                           ggCAGCAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCa    <  1:8307856/51‑1 (MQ=255)
                                            gcagcaGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAg   >  1:2727651/1‑51 (MQ=255)
                                            gcagcaGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAg   <  1:7403572/51‑1 (MQ=255)
                                            gcagcaGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAg   <  1:8663418/51‑1 (MQ=255)
                                             cagcagAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAGc  >  1:4273599/1‑51 (MQ=255)
                                             cagcagAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAGc  >  1:4274854/1‑51 (MQ=255)
                                                |                                               
TTACCGATCATGGTGCCCAGCACTTCGTTATCGCGCGGGTTCAGGCAGTAGAGAATGGTTTTCGGCAGCAGGTTTTCTTCATTCTGCTTGCTCAGC  >  REL606/3176682‑3176777

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: