Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,342,867 T→G 13.4% Q299P (CAG→CCG)  nuoG ← NADH dehydrogenase subunit G

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,342,8670T→G13.4% 16.3 ‑0.3 ‑0.459896nuoGNADH dehydrogenase subunit G
Reads supporting (aligned to +/- strand):  new base (5/3):  ref base (27/24):  total (32/27)
Fisher's exact test strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.28e-01

ACTGACGCAGAATATCTGCTGCGCCCTGCATTGCCTGTTCGGCGTTGAGGGTAATGAAATCATCGCCAC  >  REL606/2342832‑2342900
                                   |                                 
aCTGACGCAGAATATCTGCTGCGCCCTGCATTGCCt                                   <  1:5650239/36‑1 (MQ=255)
aCTGACGCAGAATATCTGCTGCGCCCTGCATTGCCt                                   >  1:3587047/1‑36 (MQ=255)
aCTGACGCAGAATATCTGCTGAGCCCTGCATTGCCt                                   >  1:5407036/1‑36 (MQ=255)
  tcaCGCATAATATCTGCTGCGCCCTGCATTGCCTGt                                 <  1:984580/34‑1 (MQ=37)
  tGACGCAGAATATCTGCTGCGCCCTGCATTGCCTGt                                 <  1:2304841/36‑1 (MQ=255)
  tGACGCAGAATATCTGCTGCGCCCTGCATTGCCTGt                                 >  1:7851303/1‑36 (MQ=255)
  tGACGCAGAATATATGCTGCGCACTGCATTGCCTGt                                 >  1:7615681/1‑36 (MQ=38)
    aCGCAGAATATCTGCTGCGCCCTGCATTGCCTGTTc                               <  1:5680777/36‑1 (MQ=255)
    aCGCAGAATATCTGCTGCGCCCTGCATTGCCTGTTc                               <  1:2288413/36‑1 (MQ=255)
    aCGCAGAATATCTGCTGCGCCCTGCATTGCCTGTTc                               <  1:4961156/36‑1 (MQ=255)
     cGCAGAATATCTGCTGCGCCCTGCATTGCCGGTTCg                              <  1:2829070/36‑1 (MQ=255)
      gCAGAATATCTGCTGCGCCCTGCATTGCCTGTTCgg                             <  1:5063841/36‑1 (MQ=255)
      gCAGAATATCTGCTGCGCCCTGCATTGCCTGTTCgg                             >  1:4787574/1‑36 (MQ=255)
      gCAGAATATCTGCTGCGCCCTGCATTGCCTGTTCgg                             <  1:6364471/36‑1 (MQ=255)
       cAGAATATCTGCTGCGCCCTGCATTGCCTGTTCGGc                            >  1:2188250/1‑36 (MQ=255)
       cAGAATATCTGCTGCGCCCTGCATTGCCTGTTCGGc                            >  1:494825/1‑36 (MQ=255)
       cAGAATATCTGCTGCGCCCTGCATTGCCTGTTCGGc                            >  1:4637708/1‑36 (MQ=255)
       cAGAATATCTGCTGCGCCCTGCATTGCCGGTTCGGc                            <  1:6168314/36‑1 (MQ=255)
        aGAATATCTGCTGCGCCCTGCATTGCCTGTTCGGCg                           >  1:4081480/1‑36 (MQ=255)
         gAATATCTGCTGCGTCCTGCATTGCCGGTTAGGCGt                          >  1:6816648/1‑36 (MQ=21)
          aaTATCTGCTGCGCCCTGCATTGCCTGTTCGGCGtt                         <  1:7337638/36‑1 (MQ=255)
          aaTATCTGCTGCGCCCTGCATTGCCTGTTCGGCGtt                         <  1:668420/36‑1 (MQ=255)
           atatCTGCTGCGCCCTGCATTGCCTGTTCGGCGTTg                        <  1:3002175/36‑1 (MQ=255)
            tatCTGCTGCGCCCTGCATTGCCTGTTCGGCGTTGa                       >  1:127210/1‑36 (MQ=255)
            tatCTGCTGCGCCCTGCATTGCCTGTTCGGCGTTGa                       >  1:1870640/1‑36 (MQ=255)
               ctgctgCGCCCTGCATTGCCTGTTCGGCGTTGAggg                    <  1:7688660/36‑1 (MQ=255)
               ctgctgCGCCCTGCATTGCCTGTTCGGCGTTGAggg                    >  1:7146164/1‑36 (MQ=255)
               ctgctgCGCCCTGCATTGCCTGTTCGGCGTTGAggg                    >  1:4295586/1‑36 (MQ=255)
               ctgctgCGCCCTGCATTGCCGGTTCGGCGTTGAggg                    >  1:429652/1‑36 (MQ=255)
                tgctgcGCCCTGCATTGCCTGTTCGGCGTTGAGGGt                   <  1:4323469/36‑1 (MQ=255)
                 gctgcGCCCTGCATTGCCGGTTCGGCGTTGAGGGTa                  >  1:3598692/1‑36 (MQ=255)
                  ctgcGCCCTGCATTGCCTGTTCGGCGTTGAGGGTaa                 <  1:417971/36‑1 (MQ=255)
                  ctgcGCCCTGCATTGCCTGTTCGGCGTTGAGGGTaa                 >  1:5214756/1‑36 (MQ=255)
                  ctgcGCCCTGCATTGCCTGTTCGGCGTTGAGGGTaa                 >  1:6034920/1‑36 (MQ=255)
                   tgcGCCCTGCATTGCCTGTTCGGCGTTGAGGGTAAt                <  1:3601669/36‑1 (MQ=255)
                   tgcGCCCTGCATTGCCTGTTCGGCGTTGAGGGTAAt                >  1:5190153/1‑36 (MQ=255)
                   tgcGCCCTGCATTGCCTGTTCAGCGTTGAGGGTAAt                >  1:7692373/1‑36 (MQ=255)
                   tgcGCCCTGCATTGCCGGTTCGGCGTTGAGGGTAAt                >  1:6284339/1‑36 (MQ=255)
                     cgcCCTGCATTGCCTGTTCGGCGTTGAGGGTAATGa              <  1:4529135/36‑1 (MQ=255)
                     cgcCCTGCATTGCCTGTTCGGCGTTGAGGGTAATGa              <  1:7161972/36‑1 (MQ=255)
                     cgcCCTGCATTGCCTGTTCGGCGTTGAGGGTAATGa              >  1:3151873/1‑36 (MQ=255)
                     cgcCCTGCATTGCCTGTTCGGCGTTGAGGGTAATGa              >  1:3073100/1‑36 (MQ=255)
                     cgcCCTGCATTGCCTGTTCGGCGTTGAGGGTAATGa              <  1:2348057/36‑1 (MQ=255)
                      gcCCTGCATTGCCTGTTCGGCGTTGAGGGTAATGaa             <  1:7678140/36‑1 (MQ=255)
                        ccTGCATTGCCTGTTCGGCGTTGAGGGTAATGAAat           >  1:5959333/1‑36 (MQ=255)
                        ccTGCATTGCCGGTTCGGCGTTGAGGGTAATGAAat           >  1:5796841/1‑36 (MQ=255)
                         cTGCATTGCCTGTTCGGCGTTGAGGGTAATGAAatc          >  1:863230/1‑36 (MQ=255)
                          tGCATTGCCTGTTCGGCGTTGAGGGTAATGAAatca         >  1:5944060/1‑36 (MQ=255)
                          tGCATTGCCTGTTCGGCGTTGAGGGTAATGAAatca         >  1:4131629/1‑36 (MQ=255)
                           gCATTGCCTGTTCGGCGTTGAGGGTAATGAAatcat        >  1:5426274/1‑36 (MQ=255)
                           gCATTGCCTGTTCGGCGTTGAGGGTAATGAAatcat        <  1:437238/36‑1 (MQ=255)
                           gCATTGCCGGTTCGGCGTTGAGGGTAATGAAatcat        <  1:1428773/36‑1 (MQ=255)
                             aTTGCCTGTTCGGCGTTGAGGGTAATGAAATCATCg      >  1:3935708/1‑36 (MQ=255)
                             aTTGCCTGTTCGGCGTTGAGGGTAATGAAATCATCg      >  1:1672554/1‑36 (MQ=255)
                             aTTGCCTGTTCGGCGTTGAGGGTAATGAAATCATCg      <  1:3238757/36‑1 (MQ=255)
                             aTTGCCTGTTCGGCGTTGAGGGTAATGAAATCATCg      >  1:1905350/1‑36 (MQ=255)
                               tGCCTGTTCGGCGTTGAGGGTAATGAAATCATCGcc    >  1:2760955/1‑36 (MQ=255)
                               tGCCTGTTCGGCGTTGAGGGTAATGAAATCATCGcc    <  1:101419/36‑1 (MQ=255)
                                gCCTGTTCGGCGTTGAGGGTAATGAAATCATCGCCa   <  1:4619181/36‑1 (MQ=255)
                                gCCTGTTCGGCGTTGAGGGTAATGAAATCATCGCCa   <  1:1523477/36‑1 (MQ=255)
                                 ccTGTTCGGCGTTGAGGGTAATGAAATCATCGCCac  >  1:2111785/1‑36 (MQ=255)
                                 ccTGTTCGGCGTTGAGGGTAATGAAATCATCGCCac  <  1:2089493/36‑1 (MQ=255)
                                   |                                 
ACTGACGCAGAATATCTGCTGCGCCCTGCATTGCCTGTTCGGCGTTGAGGGTAATGAAATCATCGCCAC  >  REL606/2342832‑2342900

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG