Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 737,252 C→A 4.2% I185I (ATC→ATA sdhA → succinate dehydrogenase flavoprotein subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606737,2520C→A 4.2% 2.2 0.0 ‑0.0137I185I (ATC→ATAsdhAsuccinate dehydrogenase flavoprotein subunit
Reads supporting (aligned to +/- strand):  new base (3/3):  ref base (56/74):  total (60/77)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 143 total aligned reads displayed.

ATGGCGCGGTGGTGGGTTGTACCGCACTGTGCATCGAAACCGGTGAAGTGGTTTATTTCAAAGCCCGCG  >  REL606/737218‑737286
                                  |                                  
aTGGCGCGGTGGTGGGTTGTACCGCACTGTGCATc                                    <  1:9513906/35‑1 (MQ=255)
aTGGCGCGGTGGTGGGTTGTACCGCACTGTGCATc                                    <  1:18904270/35‑1 (MQ=255)
 tGGCGCGGTGGTGGGTTGTACCGCACTGTGCATCg                                   <  1:8341898/35‑1 (MQ=255)
 tGGCGCGGTGGTGGGTTGTACCGCACTGTGCATCg                                   <  1:17301669/35‑1 (MQ=255)
 tGGCGCGGTGGTGGGTTGTACCGCACTGTGCATCg                                   <  1:13040051/35‑1 (MQ=255)
  ggCGCGGTGGTGGGTTGTACCGCACTGTGCATCGa                                  <  1:7782918/35‑1 (MQ=255)
  ggCGCGGTGGTGGGTTGTACCGCACTGTGCATCGa                                  <  1:6407907/35‑1 (MQ=255)
  ggCGCGGTGGTGGGTTGTACCGCACTGTGCATCGa                                  <  1:1894154/35‑1 (MQ=255)
   gcgcGGTGGTGGGTTGTACCGCACTGTGCATCGaa                                 <  1:15024005/35‑1 (MQ=255)
   gcgcGGTGGTGGGTTGTACCGCACTGTGCATCGaa                                 >  1:6106155/1‑35 (MQ=255)
   gcgcGGTGGTGGGTTGTACCGCACTGTGCATCGaa                                 <  1:6116250/35‑1 (MQ=255)
   gcgcGGTGGTGGGTTGTACCGCACTGTGCATAGaa                                 >  1:10972369/1‑35 (MQ=255)
    cgcgGTGGTGGGTTGTACCGCACTGTGCATCGaaa                                <  1:2034037/35‑1 (MQ=255)
     gcgGTGGTGGGTTGTACCGCACTGTGCATCGAAAc                               <  1:5097875/35‑1 (MQ=255)
     gcgGTGGTGGGTTGTACCGCACTGTGCATCGAAAc                               <  1:15612169/35‑1 (MQ=255)
     gcgGTGGTGGGTTGTACCGCACTGTGCATAGAAAc                               <  1:2585284/35‑1 (MQ=255)
      cgGTGGTGGGTTGTACCGCACTGTGCATCGAAAcc                              <  1:4640142/35‑1 (MQ=255)
      cgGTGGTGGGTTGTACCGCACTGTGCATCGAAAcc                              <  1:2797774/35‑1 (MQ=255)
       ggtggtGGGTTGTACCGCACTGTGCATGGAAACCg                             >  1:19371731/1‑35 (MQ=255)
       ggtggtGGGTTGTACCGCACTGTGCATCGAAACCg                             <  1:5483398/35‑1 (MQ=255)
       ggtggtGGGTTGTACCGCACTGTGCATCGAAACCg                             >  1:5264706/1‑35 (MQ=255)
       ggtggtGGGTTGTACCGCACTGTGCATCGAAACCg                             >  1:11126561/1‑35 (MQ=255)
       ggtggtGGGTTGTACCGCACTGTGCATCGAAACCg                             <  1:13555693/35‑1 (MQ=255)
       ggtggtGGGGTGTACCGCACTGTGCATCGAAACCg                             >  1:8408975/1‑35 (MQ=255)
        gtggtgGGTTGTACCGCACTGTGCATCGAAACCgg                            >  1:18771309/1‑35 (MQ=255)
        gtggtgGGTTGTACCGCACTGTGCATCGAAACCgg                            >  1:13716341/1‑35 (MQ=255)
        gtggtgGGTTGTACCGCACTGTGCATCGAAACCgg                            <  1:10485510/35‑1 (MQ=255)
         tggtggGTTGTACCGCACTGTGCATCGAAACCGGt                           >  1:12507967/1‑35 (MQ=255)
          ggtggGTTGTACCGCACTGTGCATCGAAACCGGTg                          <  1:7437166/35‑1 (MQ=255)
          ggtggGTTGTACCGCACTGTGCATCGAAACCGGTg                          <  1:4329620/35‑1 (MQ=255)
            tggGTTGTACCGCACTGTGCATCGAAACCGGTGaa                        <  1:2320249/35‑1 (MQ=255)
            tggGTTGTACCGCACTGTGCATCGAAACCGGTGaa                        <  1:2367515/35‑1 (MQ=255)
            tggGTTGTACCGCACTGTGCATCGAAACCGGTGaa                        >  1:4539361/1‑35 (MQ=255)
            tggGTTGTACCGCACTGTGCATCGAAACCGGTGaa                        <  1:15857310/35‑1 (MQ=255)
            tggGTTGTACCGCACTGTGCATCGAAACCGGTGaa                        <  1:5844845/35‑1 (MQ=255)
            tggGTTGTACCGCACTGTGCATCGAAACCGGTGaa                        <  1:10579309/35‑1 (MQ=255)
            tggGTTGTACCGCACTGTGCATCGAAACCGGTGaa                        >  1:14549127/1‑35 (MQ=255)
            tggGTTGTACCGAACTGTGCATCGAAACCGGTGaa                        <  1:12365752/35‑1 (MQ=255)
             gggTTGTACCGCACTGTGCATCGAAACCGGTGAAg                       >  1:12320761/1‑35 (MQ=255)
             gggTTGTACCGCACTGTGCATCGAAACCGGTGAAg                       <  1:14335751/35‑1 (MQ=255)
             gggTTGTACCGCACTGTGCATCGAAACCGGTGAAg                       <  1:13239354/35‑1 (MQ=255)
               gTTGTACCGCACTGTGCATCGAAACCGGTGAAGTg                     >  1:447508/1‑35 (MQ=255)
                ttGTACCGCACTGTGCATCGAAACCGGTGAAGTgg                    <  1:3784811/35‑1 (MQ=255)
                ttGTACCGCACTGTGCATCGAAACCGGTGAAGTgg                    >  1:10237397/1‑35 (MQ=255)
                 tGTACCGCACTGTGCATCGAAACCGGTGAAGTGGt                   >  1:1636201/1‑35 (MQ=255)
                 tGTACCGCACTGTGCATCGAAACCGGTGAAGTGGt                   >  1:1700000/1‑35 (MQ=255)
                 tGTACCGCACTGTGCATCGAAACCGGTGAAGTGGt                   >  1:9450082/1‑35 (MQ=255)
                 tGTACCGCACTGTGCATCGAAACCGGTGAAGTGGt                   <  1:2926332/35‑1 (MQ=255)
                 tGTACCGCACTGTGCATCGAAACCGGTGAAGTGGt                   >  1:15588755/1‑35 (MQ=255)
                   tACCGCACTGTGCATCGAAACCGGTGAAGTGGttt                 >  1:18166948/1‑35 (MQ=255)
                     ccGCACTGTGCATCGAAACCGGTGAAGTGGTTTAt               <  1:753200/35‑1 (MQ=255)
                      cGCACTGTGCATCGAAACCGGTGAAGTGGTTTAtt              <  1:2230351/35‑1 (MQ=255)
                      cGCACTGTGCATCGAAACCGGTGAAGTGGTTTAtt              >  1:7279826/1‑35 (MQ=255)
                      cGCACTGTGCATCGAAACCGGTGAAGTGGTTTAtt              <  1:19933569/35‑1 (MQ=255)
                      cGCACTGTGCATCGAAACCGGTGAAGTGGTTTAtt              >  1:14839678/1‑35 (MQ=255)
                       gCACTGTGCATCGAAACCGGTGAAGTGGTTTAttt             >  1:19611560/1‑35 (MQ=255)
                       gCACTGTGCATCGAAACCGGTGAAGTGGTTTAttt             >  1:16313331/1‑35 (MQ=255)
                       gCACTGTGCATCGAAACCGGTGAAGTGGTTTAttt             <  1:20324386/35‑1 (MQ=255)
                       gCACTGTGCATCGAAACCGGTGAAGTGGTTTAttt             >  1:3327971/1‑35 (MQ=255)
                       gCACTGTGCATAGAAACCGGTGAAGTGGTTTAttt             <  1:13499946/35‑1 (MQ=255)
                        cACTGTGCATCGAAACCGGTGAAGTGGTTTATTTc            <  1:16298407/35‑1 (MQ=255)
                        cACTGTGCATCGAAACCGGTGAAGTGGTTTATTTc            <  1:4104603/35‑1 (MQ=255)
                        cACTGTGCATCGAAACCGGTGAAGTGGTTTATTTc            <  1:11661831/35‑1 (MQ=255)
                        cACTGTGCATCGAAACCGGTGAAGTGGTTTATTTc            <  1:18146023/35‑1 (MQ=255)
                        cACTGTGCATCGAAACCGGTGAAGTGGTTTATTTc            <  1:3658635/35‑1 (MQ=255)
                        cACTGTGCATCGAAACCGGTGAAGTGGTTTATTTc            <  1:18535072/35‑1 (MQ=255)
                         aCTGTGCATCGAAACCGGTTAAGTGGTTTATTTCa           >  1:1685032/1‑35 (MQ=255)
                         aCTGTGCATCGAAACCGGTGAAGTGGTTTATTTCa           >  1:3334505/1‑35 (MQ=255)
                         aCTGTGCATCGAAACCGGTGAAGTGGTTTATTTCa           >  1:16320358/1‑35 (MQ=255)
                          cTGTGCATCGAAACCGGTGAAGTGGTTTATTTCaa          >  1:17128832/1‑35 (MQ=255)
                          cTGTGCATCGAAACCGGTGAAGTGGTTTATTTCaa          <  1:6740826/35‑1 (MQ=255)
                          cTGTGCATCGAAACCGGTGAAGTGGTTTATTTCaa          <  1:20179319/35‑1 (MQ=255)
                           tgtgCATCTAAACCGGTGAAGTGGTTTATTTCaaa         >  1:15073823/1‑35 (MQ=255)
                           tgtgCATCGAAACCGGTGAAGTGGTTTATTTCaaa         <  1:18730026/35‑1 (MQ=255)
                           tgtgCATCGAAACCGGTGAAGTGGTTTATTTCaaa         >  1:17736857/1‑35 (MQ=255)
                           tgtgCATCGAAACCGGTGAAGTGGTTTATTTCaaa         <  1:8427827/35‑1 (MQ=255)
                           tgtgCATCGAAACCGGTGAAGTGGTTTATTTCaaa         >  1:1520234/1‑35 (MQ=255)
                            gtgCATCGAAACCGGTGAAGTGGTTTATTTCAAAg        <  1:2084901/35‑1 (MQ=255)
                            gtgCATCGAAACCGGTGAAGTGGTTTATTTCAAAg        >  1:1616804/1‑35 (MQ=255)
                             tgCATCGAAACCGGTGAAGTGGTTTATTTCAAAGc       >  1:11165295/1‑35 (MQ=255)
                             tgCATCGAAACCGGTGAAGTGGTTTATTTCAAAGc       >  1:3803697/1‑35 (MQ=255)
                             tgCATCGAAACCGGTGAAGTGGTTTATTTCAAAGc       >  1:1319704/1‑35 (MQ=255)
                              gCATCGAAACCGGTGAAGTGGTTTATTTCAAATcc      >  1:20047854/1‑35 (MQ=255)
                              gCATCGAAACCGGTGAAGTGGTTTATTTCAAAGcc      <  1:5702505/35‑1 (MQ=255)
                              gCATCGAAACCGGTGAAGTGGTTTATTTCAAAGcc      <  1:13818304/35‑1 (MQ=255)
                              gCATCGAAACCGGTGAAGTGGTTTATTTCAAAGcc      <  1:10332666/35‑1 (MQ=255)
                              gCATCGAAACCGGTGAAGTGGTTTATTTCAAAGcc      <  1:13330605/35‑1 (MQ=255)
                              gCATCGAAACCGGTGAAGTGGTTTATTTCAAAGcc      >  1:13208059/1‑35 (MQ=255)
                               cATCGAAACCGGTGAAGTGGTTTATTTCAAAGccc     >  1:17816531/1‑35 (MQ=255)
                               cATCGAAACCGGTGAAGTGGTTTATTTCAAAGccc     >  1:18147066/1‑35 (MQ=255)
                               cATCGAAACCGGTGAAGTGGTTTATTTCAAAGccc     <  1:17623797/35‑1 (MQ=255)
                               cATCGAAACCGGTGAAGTGGTTTATTTCAAAGcac     <  1:12634063/35‑3 (MQ=255)
                                aTCGAAACCGGTGAAGTGGTTTATTTCAAAGCCcg    >  1:18796495/1‑35 (MQ=255)
                                aTCGAAACCGGTGAAGTGGTTTATTTCAAAGCCcg    <  1:10384684/35‑1 (MQ=255)
                                aTCGAAACCGGTGAAGTGGTTTATTTCAAAGCCcg    <  1:641782/35‑1 (MQ=255)
                                aTCGAAACCGGTGAAGTGGTTTATTTCAAAGCCcg    <  1:3932434/35‑1 (MQ=255)
                                aTCGAAACCGGTGAAGTGGTTTATTTCAAAGCCcg    <  1:4811031/35‑1 (MQ=255)
                                 tCGAAACCGGTTAAGTGGTTTATTTCAAAGCCcgc   >  1:13816712/1‑35 (MQ=255)
                                 tCGAAACCGGTGAAGTGGTTTATTTCAAAGCCcgc   <  1:6684462/35‑1 (MQ=255)
                                  cGAAACCGGTGAAGTGGTTTATTTCAAAGCCcgcg  <  1:11640535/35‑1 (MQ=255)
                                  |                                  
ATGGCGCGGTGGTGGGTTGTACCGCACTGTGCATCGAAACCGGTGAAGTGGTTTATTTCAAAGCCCGCG  >  REL606/737218‑737286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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.