Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,805,203 C→A 14.0% V108L (GTA→TTA)  astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,805,2030C→A14.0% 8.2 ‑5.5 0.070V108L (GTA→TTA) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  new base (8/2):  ref base (18/42):  total (26/44)
Fisher's exact test strand distribution p-value = 4.09e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGC  >  REL606/1805169‑1805237
                                  |                                  
tCGGCAGATGGCGCGATCGTTGCCGCATTGGCTAc                                    <  1:1858187/35‑1 (MQ=255)
tCGGCAGATGGCGCGATCGTTGCCGCATTGGCTAc                                    <  1:2381087/35‑1 (MQ=255)
tCGGCAGATGGCGCGATCGTTGCCGAATTGGCTAc                                    <  1:15579314/35‑1 (MQ=255)
 cGGCGGATGGCGCGATCGTTGCCGCATTGGCTAcc                                   <  1:15272569/35‑1 (MQ=255)
 cGGCAGATGGCGCGATCGTTGCCGCATTGGCTAcc                                   <  1:19075228/35‑1 (MQ=255)
  ggCAGATTGCGCGATCGTTGCCGCATTGGCTAacc                                  >  1:3277291/1‑35 (MQ=37)
  ggCAGATGGCGCGATCGTTGCCGCATTGGCTAccc                                  <  1:3566454/35‑1 (MQ=255)
  ggCAGATGGCGCGATCGTTGCCGCATTGGCTAccc                                  <  1:1309637/35‑1 (MQ=255)
  ggCAGATGGCGCGATCGTTGCCGCATTGGCTAcac                                  >  1:7233774/1‑33 (MQ=255)
  ggCAGATGGCGCGATCGTTGCCGAATTGGCTAcca                                  >  1:2770488/1‑34 (MQ=255)
     aGATGGCGCGATCGTTGCCGCATTGGCTACCccca                               >  1:14850811/1‑35 (MQ=255)
     aGATGGCGCGATCGTTGCCGCATTGGCTAACcaca                               >  1:9958504/1‑35 (MQ=255)
     aGATGGCGCGATCGTTGCCGCATTAGCTACCcaca                               <  1:18014753/35‑1 (MQ=255)
      gATGGCGCGATCGTTGCCGCATTGGCTACCCACAt                              <  1:523776/35‑1 (MQ=255)
      gATGGCGCGATCGTTGCCGCATTGGCTACCCACAt                              <  1:18225284/35‑1 (MQ=255)
      gATGGCGCGATCGTTGCCGCATTGGCTACCCACAt                              <  1:5076669/35‑1 (MQ=255)
      gATGGCGCGATCGTTGCAGCATTGGCTACCCACAt                              <  1:8602817/35‑1 (MQ=255)
      gATAGCGCGATCGTTGCCGCATTGGCTACCCACAt                              <  1:13796190/35‑1 (MQ=255)
       aTGGCGCGATCGTTGCCGCATTGGCTACCCACAtt                             <  1:5149817/35‑1 (MQ=255)
       aTGGCGCGATCGTTGCCGCATTGGCTAACCCCAtt                             >  1:4282269/1‑35 (MQ=37)
        tGGCGCGATCGTTGCCGCATTGGCTAACCACATTg                            >  1:4995712/1‑35 (MQ=255)
        tGGCGCGATCGTTGCCGCATTGGCTAACCACATTg                            >  1:3575820/1‑35 (MQ=255)
        gggCGCGATCTTTGCCGCATTGGCTACCCACATTg                            <  1:15104561/34‑1 (MQ=255)
        gggCGCGATCGTTGCCGCATTGGCTACCCACATTg                            <  1:3195045/34‑1 (MQ=255)
         ggCGCGATCGTTGCCGCATTGGCTACCCACATTgg                           >  1:12083031/1‑35 (MQ=255)
         ggCGCGATCGTTGCCGCATTGGCTACACACATTgg                           <  1:12549584/35‑1 (MQ=255)
          gcgcGTTCGTTGCCGCATTGGCTACCCACATTGGc                          <  1:5211650/35‑1 (MQ=255)
          gcgcGCTCGTTGCCCCATTGGCTACCCACATTGGc                          <  1:4314880/35‑1 (MQ=38)
          gcgcGATCTTTGCCGCATTGGCTACCCACATTGGc                          <  1:16095777/35‑1 (MQ=255)
          gcgcGATCGTTGCCGCATTGGCTACCCACATTGGc                          >  1:19951079/1‑35 (MQ=255)
          gcgcGATCGTTGCCGCATTGCCTAACCACATTGGc                          >  1:4883144/1‑35 (MQ=37)
           cgcgATCGTTGCCGCATTTGCTACCCACATTGGCg                         >  1:11268019/1‑35 (MQ=255)
           cgcgATCGTTGCCGCATTGGCTACCCCCATTGGCg                         >  1:9874174/1‑35 (MQ=255)
           cgcgATCGTTGCCGCATTGGCTACCCACATTGGCg                         >  1:15294797/1‑35 (MQ=255)
           cgcgATCGTTGCCGCATTGGCTACCCACAATGGCg                         >  1:4668477/1‑35 (MQ=255)
           cgcgATCGTTCCCGCATTGGCTACCCACATTGGCg                         <  1:18335028/35‑1 (MQ=255)
            gcgATCGTTGCCGCATTGGCTACCCACATTGGCGa                        <  1:3482936/35‑1 (MQ=255)
            gcgATCGTTGCCGCATTGGCTACCCACATTGGCGa                        >  1:9401598/1‑35 (MQ=255)
            gcgATCGTTGCCGCATTGGCTACCCACATTGGCGa                        >  1:12852741/1‑35 (MQ=255)
            gcgATCGTTGCCGCATTGGCTAACCACATTggggg                        >  1:19152966/1‑32 (MQ=37)
            gcgATCGTTGCCGCATTAGCTACCCACATTGGCGa                        <  1:62140/35‑1 (MQ=255)
             ctaTCGTTGCCGCATTGGCTACCCACATTGGCGaa                       <  1:1704970/33‑1 (MQ=255)
             cgATCTTTGCCGCATTGGCTACCCACATTGGCGag                       <  1:19007968/35‑2 (MQ=255)
             cgATCGTTGCCGCATTGGCTACCCACATTGGCGaa                       <  1:6957291/35‑1 (MQ=255)
             cgATCGTTGCCGCATTGGCTACCCACATTGGCGaa                       <  1:10772896/35‑1 (MQ=255)
              gATCGTTGCCGCATTGGCTACCCACATTGGCGAAg                      >  1:6990297/1‑35 (MQ=255)
              gATCGTTGCCGCATTGGCTACCCACATTGGCGAAg                      >  1:12078833/1‑35 (MQ=255)
               aTCTTTGCCGCATTGGCTACCCACATTGGCGAAGc                     <  1:18128544/35‑1 (MQ=255)
               aTCGTTGCCGCATTGGCTACCCCCATTGGCGAAGc                     >  1:19465782/1‑35 (MQ=255)
               aTCGTTGCCGCATTGGCTACCCACATTGGCGAAGc                     <  1:5463139/35‑1 (MQ=255)
               aTCGTTGCCGCATTGGCTAACCACATTGGCGAAGc                     >  1:3100606/1‑35 (MQ=255)
                tctttGCCGCATTGGCTACCCACATTGGCGAAGCg                    <  1:13834720/32‑1 (MQ=255)
                tctttGCCGCATTGGCTACCCACATTGGCGAAGCg                    <  1:18117822/32‑1 (MQ=255)
                tCGTTGCCGCATTGGCTACCCACATTGGCGAAGCg                    <  1:13456260/35‑1 (MQ=255)
                tCGTTGCCCCATTGGCTACCCACATTGGCGAAGCg                    <  1:11317601/35‑1 (MQ=255)
                   ttGCCGCATTGGCTACCCACATTGGCGAAGCGGaa                 <  1:8227839/35‑1 (MQ=255)
                    tGCCGCATTGGCTACCCACATTGGCGAAGCGGAAc                <  1:19467861/35‑1 (MQ=255)
                     gCCGCATTGGCTACCCACATTGGCGAAGCGGAACt               <  1:20021806/35‑1 (MQ=255)
                     gCCGCATTGGCTACCCACATTGGCGAAGCGGAACt               >  1:635948/1‑35 (MQ=255)
                      ccccATTGGCTACCCACATTGGCGAAGCGGAACTg              <  1:10751026/32‑1 (MQ=255)
                      ccccATTGGCTACCAACATTGGCGAAGCGGAACTg              <  1:3483997/32‑1 (MQ=37)
                      ccGCATTGGCTACCCACATTGGCGAAGCGGAACTg              <  1:15026994/35‑1 (MQ=255)
                      ccGCATTGGCTACCCACATTGGCGAAGCGGAACTg              <  1:11720176/35‑1 (MQ=255)
                      ccGCATTGGCTACCCACATTGGCGAAGCGGAACTg              <  1:17279122/35‑1 (MQ=255)
                      ccGCATTGGCTAACCACATTGGCGAAGCGGAACTg              >  1:5438817/1‑35 (MQ=255)
                       cccATTGGCTACCCACATTGGCGAAGCGGAACTGa             <  1:450436/33‑1 (MQ=255)
                       cccATTGGCTAACCACATTGGCGAAGCGGAACTGa             <  1:1816252/33‑1 (MQ=38)
                       cGCGTTGGCTAACCACATTGGCGAAGCGGAACTGa             <  1:18380520/35‑1 (MQ=37)
                       cGCATTGGCTACCCACATTGGCGAAGCGGAACTGa             >  1:6902510/1‑35 (MQ=255)
                       cGCATTGGCTAACCACATTGGCGAAGCGGAAATGa             >  1:10755898/1‑35 (MQ=37)
                        gCATTGGCTACCCACATTGGCGAAGCGGAACTGAc            <  1:13488004/35‑1 (MQ=255)
                        gCATTGGCTACCCACATTGGCGAAGCGGAACTGAc            <  1:13782024/35‑1 (MQ=255)
                        gCATTGGCTACCCACATTGGCGAAGCGGAACTGAc            <  1:3962879/35‑1 (MQ=255)
                        ccATTGGCTACCCACATTGGCGAAGCGGAACTGAc            <  1:7659671/34‑1 (MQ=255)
                        ccATTGGCTACCCACATTGGCGAAGCGGAACTGAc            <  1:19257884/34‑1 (MQ=255)
                         ctttGGCTCCCCACATTGGCGAAGCGGAACTGACg           <  1:19134809/33‑1 (MQ=38)
                         cATTGGCTACCCAAATTGGCGAAGCGGAACTGACg           <  1:10759285/35‑1 (MQ=255)
                         cATTGGCTAACCACATTGGCGGAGCGGGAATGACg           >  1:16900110/1‑35 (MQ=16)
                         cATTGGCTAACCACATTGGCGAAGCGGAACTGACg           >  1:5662526/1‑35 (MQ=255)
                         cATTGGCTAACCACATTGGCGAAGCGGAAATGACg           >  1:6304006/1‑35 (MQ=37)
                          aTTGGCTCCCCACATTGGCGAAGCGGAACTGACGc          <  1:15161544/35‑1 (MQ=255)
                          aTTGGCTACCCACATTGGCGAAGCGGAACTGAAGc          <  1:7991166/35‑1 (MQ=255)
                            tGGCTACCCACATTGGCGAAGCGGAACTGACGCTg        >  1:14867226/1‑35 (MQ=255)
                            tGGCTACCCACATTGGCGAAGCGGAACTGACGCTg        >  1:17409237/1‑35 (MQ=255)
                              gCTACCCACATTGGCGAAGCGGAACTGACGCTGGa      >  1:13239047/1‑35 (MQ=255)
                              gCTACCCACATTGGCGAAGCGGAACTGACGCTGGa      >  1:17418491/1‑35 (MQ=255)
                              gCTACCCACATTGGCGAAGCGGAACTGACGCTGGa      >  1:2006283/1‑35 (MQ=255)
                               cTACCCACATTGGCGAAGCGGAACTGACGCTGGaa     >  1:20300638/1‑35 (MQ=255)
                                tACCCACATTGGCGAAGCGGAACTGACGCTGGaaa    >  1:8538497/1‑35 (MQ=255)
                                tACCCACATTGGCGAAGCGGAACTGACGCTGGaaa    >  1:9764871/1‑35 (MQ=255)
                                  cccACATTGGCGAAGCGGAACTGACGCTGGAAAGc  >  1:1420631/1‑35 (MQ=255)
                                  cccACATTGGCGAAGCGGAACTGACGCTGGAAAGc  >  1:13214333/1‑35 (MQ=255)
                                  |                                  
TCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACATTGGCGAAGCGGAACTGACGCTGGAAAGC  >  REL606/1805169‑1805237

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.