Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,336,064 T→G 15.4% I168L (ATC→CTC)  nuoM ← NADH dehydrogenase subunit M

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,336,0640T→G15.4% 19.8 0.0 ‑0.065502nuoMNADH dehydrogenase subunit M
Reads supporting (aligned to +/- strand):  new base (5/5):  ref base (28/26):  total (34/31)
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.77e-01

CTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACA  >  REL606/2336029‑2336099
                                   |                                   
cTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGAt                                     <  1:6951050/36‑1 (MQ=255)
cTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGAt                                     >  1:209966/1‑36 (MQ=255)
 tGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATa                                    >  1:7212305/1‑36 (MQ=255)
 tGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATa                                    >  1:6343624/1‑36 (MQ=255)
   ggTGTAAATGAAGAACTTGGTTGCCGCCGTGATAcg                                  >  1:6243000/1‑36 (MQ=255)
   ggTGTAAATGAAGAACTTGGTTGCCGCCGTGAGAcg                                  <  1:6192150/36‑1 (MQ=255)
   ggTGTAAATGAAGAACTTGGTTGCCGCCGTGAGAcg                                  >  1:765977/1‑36 (MQ=255)
    gtgtAAATGAAGAACTTGGTTGCCGACGTGAGAcgc                                 >  1:2955505/1‑36 (MQ=37)
     tgtAAATGAAGAACTTTGTTGCCGCCGTGATAcgcg                                >  1:493147/1‑36 (MQ=255)
     tgtAAATGAAGAACTTGGTTGCCGCCGTGATAcgcg                                >  1:3081062/1‑36 (MQ=255)
      gtAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                               >  1:1839912/1‑36 (MQ=255)
      gtAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                               >  1:1681536/1‑36 (MQ=255)
      gtAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                               <  1:7181525/36‑1 (MQ=255)
      gtAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                               <  1:7187337/36‑1 (MQ=255)
      gtAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                               >  1:926984/1‑36 (MQ=255)
      gcaaaTGAAGAACTGCGTTGCCGCCGTGACACGCGt                               >  1:624553/3‑36 (MQ=21)
       tAAATGAAGAACTTGGTTGCCGCCGTGATACGCGtt                              >  1:200903/1‑36 (MQ=255)
       tAAATGAAGAACTTGGTTGCCGCCGTGATACGCGtt                              <  1:584936/36‑1 (MQ=255)
        aaaTGAAGAACTTGGTTGCCGCCGTGATACGCGttt                             >  1:4911258/1‑36 (MQ=255)
        aaaTGAAGAACTTGGTTGCCGCCGTGAGACGCGttt                             <  1:4332089/36‑1 (MQ=255)
        aaaAGAAGAACTTGGTTGCCGCCGTGATCCGCGttt                             <  1:3503189/36‑1 (MQ=37)
         aaTGAAGAACTTGGTTGCCGCCGTGATTCGCGtttt                            >  1:6871687/1‑36 (MQ=255)
         aaTGAAGAACTTGGTTGCCGCCGTGATACGCGtttt                            <  1:2389255/36‑1 (MQ=255)
         aaTGAAGAACTTGGTTGCCGCCGTGATACGCGtttt                            <  1:7406522/36‑1 (MQ=255)
         aaTGAAGAACTTGGTTGCCGCCGTGATACGCGtttt                            >  1:376551/1‑36 (MQ=255)
         aaTGAAGAACTTGGTTGCCGCCGAGATACGCGtttt                            >  1:7796589/1‑36 (MQ=255)
          ataaagaaCTTGGTTGCCGCCGTGATACGCGTTTTa                           <  1:1025680/33‑1 (MQ=255)
          aTGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTa                           <  1:7369057/36‑1 (MQ=255)
          aTGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTa                           <  1:7223548/36‑1 (MQ=255)
            gaagaaCTTGGTTGCCGCCGTGATACGCGTTTTAcc                         <  1:7040336/36‑1 (MQ=255)
             aagaaCTTGGTTGCCGCCGTGATACGCGTTTTACCg                        >  1:1289776/1‑36 (MQ=255)
              agaaCTTGGTTGCCGCCGTGATACGCGTTTTACCGt                       >  1:1157696/1‑36 (MQ=255)
                aaCTTGGTTGCCGCCGTGATACGCGTTTTACCGTCa                     <  1:3605764/36‑1 (MQ=255)
                 aCTTGGTTGCCGCCGTGATACGCGTTTTACCGTCag                    >  1:5480435/1‑36 (MQ=255)
                 aCTTGGTTGCCGCCGTGAGACGCGTTTTACCGTCag                    >  1:7735432/1‑36 (MQ=255)
                  cTTGGTTGCCGCCGTGATACGCGTTTTACCGTCaga                   <  1:4359434/36‑1 (MQ=255)
                   ttGGTTGCCGCCGTGATACGCGTTTTACCGTCagag                  >  1:4005700/1‑36 (MQ=255)
                   ttGGTTGCCGCCGTGAGACGCGTTTTACCGTCagag                  >  1:7305937/1‑36 (MQ=255)
                    tGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAgg                 >  1:5367908/1‑36 (MQ=255)
                      gTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCt               >  1:5414703/1‑36 (MQ=255)
                      gTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCt               <  1:1425870/36‑1 (MQ=255)
                      gTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCt               >  1:4438729/1‑36 (MQ=255)
                      gTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCt               <  1:1913083/36‑1 (MQ=255)
                      gTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCt               <  1:6387512/36‑1 (MQ=255)
                       ttGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt              <  1:1206725/36‑1 (MQ=255)
                       ttGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt              <  1:4658298/36‑1 (MQ=255)
                        tGCCGCCGTGATACGCGTTTTACCGTCAGAGGCttt             >  1:7606990/1‑36 (MQ=255)
                           cgccgTGATACGCGTTTTACCGTCAGAGGCTTTATg          <  1:137616/36‑1 (MQ=255)
                           cgccgTGATACGCGTTTTACCGTCAGAGGCTTTATg          <  1:6023592/36‑1 (MQ=255)
                           cgccgTGATACGCGTTTTACCGTCAGAGGCTTTATg          <  1:6515906/36‑1 (MQ=255)
                           cgccgTGATACGCGTTTTACCGTCAGAGGCTTTATg          <  1:5462871/36‑1 (MQ=255)
                           cgccgTGAGACGCGTTTTACCGTCAGAGGCTTTATg          <  1:4639382/36‑1 (MQ=255)
                            gccgTGATACGCGTTTTACCGTCAGAGGCTTTATGc         >  1:7035556/1‑36 (MQ=255)
                            gccgTGATACGCGTTTTACCGTCAGAGGCTTTATGc         <  1:6893443/36‑1 (MQ=255)
                            gccgTGATACGCGTTTTACCGTCAGAGGCTTTATGc         <  1:3279675/36‑1 (MQ=255)
                            gccgTGATACGCGTTTTACCGTCAGAGGCTTTATGc         >  1:5340869/1‑36 (MQ=255)
                            gccgTGAGACGCGTTTTACCGTCAGAGGCTTTATGc         <  1:1558973/36‑1 (MQ=255)
                            gccgTGAGACGCGTTTTACCGTCAGAGGCTTTATGc         >  1:39429/1‑36 (MQ=255)
                             ccgTGATACGCGTTTTACCGTCAGAGGCTTTATGcc        >  1:3551874/1‑36 (MQ=255)
                              cgTGATACGCGTTTTACCGTCAGAGGCTTTATGccc       >  1:729279/1‑36 (MQ=255)
                              cgTGAGACGCGTTTTACCGTCAGAGGCTTTATGccc       <  1:3896485/36‑1 (MQ=255)
                               gTGATACGCGTTTTACCGTCAGAGGCTTTATGcccc      <  1:4862537/36‑1 (MQ=255)
                               gTGATACGCGTTTTACCGTCAGAGGCTTTATGcccc      <  1:5264206/36‑1 (MQ=255)
                                tGATACGCGTTTTACCGTCAGAGGCTTTATGccccc     >  1:3393655/1‑36 (MQ=255)
                                 gATACGCGTTTTACCGTCAGAGGCTTTATGCCCCca    <  1:1485470/36‑1 (MQ=255)
                                 gATACGCGTTTTACCGTCAGAGGCTTTATGCCCCca    >  1:5092667/1‑36 (MQ=255)
                                 gATACGCGTTTTACCGTCAGAGGCTTTATGCCCCca    >  1:548725/1‑36 (MQ=255)
                                   tACGCGTTTTACCGTCAGAGGCTTTATGCCCCcaca  <  1:2023370/36‑1 (MQ=255)
                                   |                                   
CTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACA  >  REL606/2336029‑2336099

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG