Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 9,972 T→G 40.1% N174T (AAC→ACC)  yaaH ← conserved inner membrane protein associated with acetate transport

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6069,9720T→G40.1% 230.1 ‑1.9 ‑0.0165N174T (AAC→ACC) yaaHconserved inner membrane protein associated with acetate transport
Reads supporting (aligned to +/- strand):  new base (31/35):  ref base (58/41):  total (89/76)
Fisher's exact test strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.51e-01

Only 100 of 179 total aligned reads displayed.

ACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACT  >  REL606/9938‑10006
                                  |                                  
aCCAATCGGCATAACGGTGCGACCAAACTGCTCGt                                    >  1:7502530/1‑35 (MQ=255)
aCCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                    <  1:16609582/35‑1 (MQ=255)
aCCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                    <  1:3975010/35‑1 (MQ=255)
aCCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                    >  1:11934010/1‑35 (MQ=255)
aCCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                    <  1:7924306/35‑1 (MQ=255)
aCCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                    <  1:8491930/35‑1 (MQ=255)
 ccAATCGGCAGAACGGTGCGACCAAACTGCTCGtt                                   <  1:5205696/35‑1 (MQ=255)
 ccAATCGGCAGAACGGTGCGACCAAACTGCTCGtt                                   <  1:6768819/35‑1 (MQ=255)
  cAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                  <  1:7083706/35‑1 (MQ=255)
  cAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                  <  1:9186253/35‑1 (MQ=255)
  cAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                  >  1:12074713/1‑35 (MQ=255)
  cAATCGGCAGAACGGTGCGACCAAACTGCTCGGTc                                  <  1:1252286/35‑1 (MQ=255)
   aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                                 >  1:10948095/1‑35 (MQ=255)
   aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                                 >  1:13727092/1‑35 (MQ=255)
   aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                                 >  1:15237997/1‑35 (MQ=255)
   aaTCGGCAGAACGGTGCGACCAAACTGCTCGGTCa                                 <  1:10076875/35‑1 (MQ=255)
    aTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                <  1:2927952/35‑1 (MQ=255)
    aTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                <  1:12746351/35‑1 (MQ=255)
    aTCGGCAGAACGGTGCGACCAAACTGCTCGGTCAg                                <  1:14739259/35‑1 (MQ=255)
     tCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                               >  1:18768618/1‑35 (MQ=255)
     tCGGCAGAACGGTGCGACCAAACTGCTCGGTCAGt                               <  1:10335872/35‑1 (MQ=255)
      cGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                              <  1:2876113/35‑1 (MQ=255)
      cGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                              <  1:13516901/35‑1 (MQ=255)
      cGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                              >  1:13199134/1‑35 (MQ=255)
      cGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                              >  1:17918566/1‑35 (MQ=255)
       ggCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                             <  1:9417143/35‑1 (MQ=255)
       ggCAGAACGGTGCGACCAAACTGCTCGGTCAGTAc                             >  1:1113033/1‑35 (MQ=255)
        gCAGAACGGTGCGACCAAACTGCTCGGTCAGTACt                            >  1:6525198/1‑35 (MQ=255)
        gCAGAACGGTGCGACCAAACTGCTCGGTCAGTACt                            <  1:14564899/35‑1 (MQ=255)
        gCAGAACGGTGCGACCAAACTGCTCGGTCAGTACt                            <  1:8384632/35‑1 (MQ=255)
        gCAGAACCGTGCGACCAAACTGCTCGTTCAGTACt                            >  1:979829/1‑35 (MQ=255)
         cAGAACGGTGCGACCAAACTGCTCGTTCAGTACtt                           >  1:4285625/1‑35 (MQ=255)
         cAGAACGGTGCGACCAAACTGCTCGTTCAGTACtt                           <  1:11562181/35‑1 (MQ=255)
         cAGAACGGTGCGACCAAACTGCTCGTTCAGTACtt                           <  1:8747597/35‑1 (MQ=255)
         cAGAACGGTGCGACCAAACTGCTCGGTCAGTACtt                           >  1:14818196/1‑35 (MQ=255)
         cAGAACGGTGCGACCAAACTGCTCGGTCAGTACtt                           >  1:16710861/1‑35 (MQ=255)
          aGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                          >  1:18889197/1‑35 (MQ=255)
          aGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                          >  1:18426722/1‑35 (MQ=255)
          aGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                          >  1:16948897/1‑35 (MQ=255)
           gAACGGTGCGACCAAACTGCTCGTTCAGGACTTCa                         <  1:6595246/35‑1 (MQ=255)
           gAACGGTGCGACCAAACTGCTCGGTCAGTACTTCa                         <  1:3285465/35‑1 (MQ=255)
           gAACGGTGCGACCAAACTGCTCGGTCAGTACTTCa                         >  1:3400671/1‑35 (MQ=255)
           gAACGGTGCGACCAAACTGCTCGGTCAGTACTTCa                         <  1:12346924/35‑1 (MQ=255)
            aaCGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                        <  1:19638790/35‑1 (MQ=255)
            aaCGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                        >  1:11713919/1‑35 (MQ=255)
            aaCGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                        >  1:4255913/1‑35 (MQ=255)
            aaCGGTGCGACCAAACTGCTCGGTCAGTACTTCAc                        <  1:15128092/35‑1 (MQ=255)
            aaCGGTGCGACCAAACTGCTCGGTCAGTACTTCAc                        <  1:17947564/35‑1 (MQ=255)
             aCGGTGCGACAAAACTGCTCGTTCAGTACTTCAcc                       <  1:2287872/35‑1 (MQ=255)
              cGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                      >  1:606555/1‑35 (MQ=255)
               ggTGCGACCAAACTGCTCGTTCAGTACTTCACCCa                     >  1:8995436/1‑35 (MQ=255)
               ggTGCGACCAAACTGCTCGGTCAGTACTTAACCCa                     <  1:10516658/35‑1 (MQ=25)
                gTGCGACCAAACTGCTCGGTCAGTACTTCACCCAt                    <  1:5471665/35‑1 (MQ=255)
                 tGCGACCAAACTGCTCGGTCAGTACTTCACCCATc                   >  1:16017513/1‑35 (MQ=255)
                 tGCGACCAAACTGCTCGGTCAGTACTTCACCCATc                   >  1:20302485/1‑35 (MQ=255)
                 tGCGACAAAACTGCTCGTTCAGTACTTCACCCATc                   <  1:18356781/35‑1 (MQ=255)
                  gCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                  >  1:7698771/1‑35 (MQ=255)
                  gCGACCAAACTGCTCGTTCAGTACTTCACCCATCg                  >  1:1128250/1‑35 (MQ=255)
                  gCGACCAAACTGCTCGGTCAGTACTTCACCCATCg                  >  1:3907893/1‑35 (MQ=255)
                   ngACCAAACTGCTCGTTCAGTACTTCACCCATCGc                 >  1:12024584/2‑35 (MQ=255)
                   cGACCAAACTGCTCGTTCAGTACTTCACCCATCGc                 >  1:17445731/1‑35 (MQ=255)
                   cGACCAAACTGCTCGTTCAGTACTTCACCCATCGc                 >  1:16141884/1‑35 (MQ=255)
                   cGACCAAACTGCTCGGTCAGTACTTCACCCATCGc                 >  1:6677627/1‑35 (MQ=255)
                    gACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                >  1:20203112/1‑35 (MQ=255)
                    gACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                <  1:11314962/35‑1 (MQ=255)
                    gACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                >  1:13995770/1‑35 (MQ=255)
                     aCCAAACTGCTCGTTCAGTACTTCACCCATCGCCa               <  1:1366581/35‑1 (MQ=255)
                     aCCAAACTGCTCGGTCAGTACTTCACCCATCGCCa               >  1:14325220/1‑35 (MQ=255)
                     aCCAAACTGCTCGGTCAGTACTTCACCCATCGCCa               <  1:11193559/35‑1 (MQ=255)
                     aCCAAACTGCTAGGTCAGTACTTCACCCATCGCCa               <  1:19957053/35‑1 (MQ=37)
                      ccAAACTGCTCGTTCAGTACTTCACCCATCGCCag              <  1:15408761/35‑1 (MQ=255)
                      ccAAACTGCTCGGTCAGTACTTCACACATCGCCag              <  1:8319223/35‑1 (MQ=37)
                       cAAACTGCTCGTTCAGTACTTCACCCATCGCCaga             >  1:20306749/1‑35 (MQ=255)
                       cAAACTGCTCGGTCAGTACTTCCCCCATCGCCaga             <  1:16066579/35‑1 (MQ=37)
                       cAAACTGCTCGGTCAGTACTTCACCCATCGCCaga             <  1:19193766/35‑1 (MQ=255)
                       cAAACTGCTCGGTCAGTACTTCACACATCGCCaga             <  1:8951380/35‑1 (MQ=37)
                        aaaCTGCTCGTTCAGTACTTCACCCATCGCCagat            >  1:10651505/1‑35 (MQ=255)
                        aaaCTGCTCGTTCAGTACTTCACCCATCGCCagat            <  1:18135006/35‑1 (MQ=255)
                        aaaCTGCTCGGTCAGTACTTCACCCATCGCCagat            <  1:10032715/35‑1 (MQ=255)
                         aaCTGCTCGTTCAGTACTTCACCCATCGCCagata           >  1:10270588/1‑35 (MQ=255)
                           cTGCTCGTTCAGTACTTCACCCATCGCCagataga         >  1:3258690/1‑35 (MQ=255)
                           cTGCTCGTTCAGTACTTCACCCATCGCCagataga         >  1:7568738/1‑35 (MQ=255)
                           cTGCTCGGTCAGTACTTCACCCATCGCCagataga         >  1:5589467/1‑35 (MQ=255)
                           cTGCTCGGTCAGTACTTCACCCATCGCCagataga         >  1:12859707/1‑35 (MQ=255)
                           cTGCTCGGTCAGTACTTCACCCATCGCCagataga         <  1:6859610/35‑1 (MQ=255)
                            tGCTCGGTCAGTACTTCACCCATCGCCagatagat        <  1:1142913/35‑1 (MQ=255)
                             gCTCGTTCAGTACTTCACCCATCGCCAGATAGAtt       >  1:17024870/1‑35 (MQ=255)
                             gCTCGTTCAGTACTTCACCCATCGCCAGATAGAtt       <  1:17617832/35‑1 (MQ=255)
                             gCTCGTTCAGTACTTCACCCATCGCCAGATAGAtt       <  1:4469249/35‑1 (MQ=255)
                             gCTCGTTCAGTACTTCACCCATCGCCAGATAGAtt       <  1:7299896/35‑1 (MQ=255)
                              cTCGGTCAGTACTTCACCCATCGCCAGATAGATTg      >  1:19498524/1‑35 (MQ=255)
                              cTCGGTCAGTACTTCACCCATCGCCAGATAGATTg      <  1:2620136/35‑1 (MQ=255)
                               tCGTTCAGTACTTCACCCATCGCCAGATAGATTGc     <  1:3841129/35‑1 (MQ=255)
                               tCGTTCAGTACTTCACCCATCGCCAGATAGATTGc     >  1:5073203/1‑35 (MQ=255)
                                cGTTCAGTACTTCACCCATCGCCAGATAGATTGCa    <  1:8335424/35‑1 (MQ=255)
                                cGTTCAGTACTTCACCCATCGCCAGATAGATTGCa    <  1:3581123/35‑1 (MQ=255)
                                cGTTCAGTACTTCACCCATCGCCAGATAGATTGCa    <  1:18636572/35‑1 (MQ=255)
                                 gTTCAGTACTTCACCCATCGCCAGATAGATTGCAc   >  1:3519382/1‑35 (MQ=255)
                                  ttCAGTACTTCACCCATCGCCAGATAGATTGCACt  <  1:19410937/35‑1 (MQ=255)
                                  ttCAGTACTTCACCCATCGCCAGATAGATTGCACt  >  1:18602946/1‑35 (MQ=255)
                                  |                                  
ACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACT  >  REL606/9938‑10006

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