Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,335,105 C→A 9.5% A9A (GCG→GCT yjdL ← predicted transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,335,1050C→A 9.5% 8.2 ‑1.6 ‑0.072A9A (GCG→GCTyjdLpredicted transporter
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (20/45):  total (24/48)
Fisher's exact test strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.28e-01

Only 100 of 101 total aligned reads displayed.

ACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTC  >  REL606/4335071‑4335139
                                  |                                  
aCTCCCAGATTTGGATCGCCACGATATAGTATATc                                    >  1:13087189/1‑35 (MQ=255)
 cccccATATTTGGATCGCCACGATATAGTATATcg                                   <  1:17168847/33‑1 (MQ=38)
 cTCCCAGATTTGGATCTCCACGATATAGTATATcg                                   <  1:3035103/35‑1 (MQ=255)
 cTCCCAGATTTGGATCGCCACGATATAGTATATcg                                   >  1:2894663/1‑35 (MQ=255)
 cTCCCAGATTTGGATCGCCACGATATAGTATATcg                                   >  1:13540305/1‑35 (MQ=255)
  tCCCAGATTTGGATCGCCACGATATAGTATATcgc                                  >  1:15191871/1‑35 (MQ=255)
  tCCCAGATTTGGATCGCCACGATATAGTATATcgc                                  >  1:18515183/1‑35 (MQ=255)
  tCCCAGATTTGGATCGCCACGATATAGTATATcgc                                  >  1:15159527/1‑35 (MQ=255)
  ccccAGTTTTGGTTCGCCACGATATAGTATATcgc                                  <  1:6270379/34‑1 (MQ=25)
  ccccAGTTTTGGTTCGCCACGATATAGTATATcgc                                  <  1:2954432/34‑1 (MQ=25)
  ccccAGATTGGGTTCGCCACGATATAGTATATcgc                                  <  1:5972860/34‑1 (MQ=37)
   cccAGTTTTGGTTCGCCACGATATAGTATATcgcg                                 <  1:7417039/35‑1 (MQ=37)
   cccAGATTTGGATCGCCACGATATAGTATATcgcg                                 >  1:7636886/1‑35 (MQ=255)
   cccAGATTTGGATCGCCACGATATAGTATATcgcg                                 >  1:8114024/1‑35 (MQ=255)
    ccAGATTTGGATCGCCACGATATAGTATATcgcgc                                >  1:19619080/1‑35 (MQ=255)
     cAGATTTGGATCGCCACGATATAGTATATcgcgcg                               >  1:10125740/1‑35 (MQ=255)
     cAGATTTGGATCGCCACGATATAGTATATcgcgcg                               >  1:18486745/1‑35 (MQ=255)
     cAGATTTGGATCGCCACGATATAGTATATcgcgcg                               >  1:13941693/1‑35 (MQ=255)
      ggttttGCATCGCCACGATATAGTATATcgcgcgc                              <  1:8639225/32‑1 (MQ=37)
      aGATTTGTTTCGCCACGATATAGTATATcgcgcgc                              <  1:743275/35‑1 (MQ=38)
      aGATTTGGATCGCCACGATATAGTATATcgcgcgc                              <  1:1512599/35‑1 (MQ=255)
      aGATTTGGATCGCCACGATATAGTATATcgcgcgc                              >  1:4409927/1‑35 (MQ=255)
      aGATTTGGATCGCCACGATATAGTATATcgcgcgc                              >  1:330005/1‑35 (MQ=255)
      aGATTGGGTTCGCCACGATATAGTATATcgcgcgc                              <  1:1164781/35‑1 (MQ=38)
      aGATGTGGTTCGCCACGATATAGTATATcgcgcgc                              <  1:13665912/35‑1 (MQ=38)
        ttttGGTTCGCCACGATATAGTATATCGCGCGCgg                            <  1:2081462/34‑1 (MQ=255)
        ttttGAATCGCCACGATATAGTATATCGCGCGCgg                            <  1:16632103/34‑1 (MQ=255)
        aTTTGGATCGCCACTATATAGTATATCGCGCGCgg                            >  1:2807872/1‑35 (MQ=255)
        aTTTGGATCGCCACGATATAGTATATcgcgcgctg                            >  1:12628123/1‑33 (MQ=255)
         tttGGATCGCCACGATATAGTATATCGCGCGCGGc                           >  1:8272315/1‑35 (MQ=255)
           tGGATCGCCACGATATAGTATATCGCGGGCGGGtg                         >  1:4218961/1‑35 (MQ=38)
            gtttCGCCACGATATAGTATATCGCGCGCGGCtgt                        <  1:13255371/32‑1 (MQ=255)
            ggttCGCCACGATATAGTATATAGCGCGCGGCtgt                        <  1:18479056/32‑1 (MQ=18)
            ggttCGCCACGATATAGTATATAGCGCGCGGCtgt                        <  1:8988435/32‑1 (MQ=18)
            ggATCGCCACGATATAGTATATCGCGCGCGGCtgg                        >  1:15143298/1‑34 (MQ=255)
             ggtCGCCACGATATAGTATATCGCGCGCGGCtgtg                       <  1:6003297/33‑1 (MQ=255)
             gATCGCCACTATATAGTATATCGCGCGCGGCtggg                       >  1:2454551/1‑33 (MQ=37)
             gATCGCCACGATATAGTATATCGCGCGCGGCtggg                       >  1:12839058/1‑33 (MQ=255)
             gATCGCCACGATATAGTATATCGCGCGCGGCtggg                       >  1:15480109/1‑33 (MQ=255)
             gATCGCCACGATATAGTATACCGCGCGCGGCtggg                       >  1:12882389/1‑33 (MQ=37)
              aTCGCCACGATATAGTATATCGCGCGCGGCtgggg                      >  1:10103130/1‑32 (MQ=255)
               tCGCCACGATATAGTATATCGCGCGCGGCTGTGAg                     <  1:20483496/35‑1 (MQ=255)
               tCGCCACGATATAGTATATCGCGCGCGGCTGTGAg                     <  1:14844983/35‑1 (MQ=255)
               tCGCCACGATATAGTATATCGCGCGCGGCTGTGAg                     <  1:8516578/35‑1 (MQ=255)
                 gcaaCGACATAGTAAATAGCGCGCGGCTGGGGggg                   >  1:14033400/4‑35 (MQ=25)
                 gCCACGATATAGTATATCGCGCGCGGCTGGGGggg                   >  1:16490361/1‑35 (MQ=25)
                 gCCACGATATAGTATATCGCGCGCGGCTGGGGggg                   >  1:13151774/1‑35 (MQ=25)
                  ccACGATATAGTATATCGCGCGCGGCTGTGAGGgt                  <  1:20581199/35‑1 (MQ=255)
                  ccACGATATAGTATATCGCGCGCGGCTGTGAGGgt                  <  1:20184121/35‑1 (MQ=255)
                  ccACGATATAGTATATCGCGCGCGGCTGTGAGGgt                  <  1:19073845/35‑1 (MQ=255)
                  ccACGATATAGTATATCGCGCGCGGCTGTGAGGgt                  <  1:10084842/35‑1 (MQ=255)
                  ccACGATATAGTATATCGCGCGCGGCTGTGAGGgt                  <  1:1607535/35‑1 (MQ=255)
                  ccACGATATAGTATATCGCGCGCGGCTGTGAGGgt                  <  1:12141949/35‑1 (MQ=255)
                  ccACGATATAGTATATCGCGCGCGGCTGTGAGGgt                  <  1:14070094/35‑1 (MQ=255)
                  ccACGATATAGTATATCGCGCGCGGCTGGGGgggg                  >  1:20528246/1‑34 (MQ=25)
                  caaCGACATAGTAAATAGCGCGCGGCTGGGGgggg                  >  1:5983508/3‑34 (MQ=25)
                   cACGATATAGTATATCGCGCGCGGCTGTGAGGgtg                 <  1:15739141/35‑1 (MQ=255)
                   cACGATATAGTATATCGCGCGCGGCTGGGGggggg                 >  1:15442780/1‑33 (MQ=21)
                   cACGATATAGTATATCGCGCGCGGCTGGGGggggg                 >  1:363440/1‑33 (MQ=21)
                   cACGATATAGTATATCGCGCGCGGCTGGGGggggg                 >  1:18360838/1‑33 (MQ=21)
                   cACGATATAGTATATCGCGCGCGGCTGGGGggggg                 >  1:13844262/1‑33 (MQ=21)
                   aaCGACATAGTAAATAGCGCGCGGCTGGGGggggg                 >  1:9353494/2‑33 (MQ=21)
                    aCGATATAGTATATCGCGCTCGGCTGGGGgggggg                >  1:8075026/1‑32 (MQ=19)
                    aCGATATAGTATATCGCGCGCTGCTGTGAcggggg                >  1:8417031/1‑32 (MQ=25)
                    aCGATATAGTATATCGCGCGCGGCTGTGAGGgtgt                <  1:16725933/35‑1 (MQ=255)
                    aCGATATAGTATATCGCGCGCGGCTGTGAGGgtgt                <  1:14891298/35‑1 (MQ=255)
                    aCGATATAGTATATCGCGCGCGGCTGTGAGGgtgt                <  1:11995007/35‑1 (MQ=255)
                    aCGACATAGTAAATAGCGCGCGGCTGGGGgggggg                >  1:573930/1‑32 (MQ=21)
                     ctatatAGTATATCGCGCGCGGCTGTGAGGGTGtt               <  1:18397024/33‑1 (MQ=255)
                     cGATATAGTATATCGCGCGCGGCTGTGAGGGTGtt               <  1:1406013/35‑1 (MQ=255)
                     cGATATAGTATATAGCGCGCGGCTGTGAGGGTGtt               <  1:7611194/35‑1 (MQ=255)
                      gATATAGTATATCGCGCGCGGCTGTGAGGGTGttt              <  1:7392552/35‑1 (MQ=255)
                       atatAGTATATCGCGCGCGGCTGTGAGGGTGtttt             <  1:18426393/35‑1 (MQ=255)
                        tataGTATATCGCGCGCGGCTGTGAGGGTGTTTTc            <  1:519341/35‑1 (MQ=255)
                        tataGTATATCGCGCGCGGCTGTGAGGGTGTTTTc            <  1:8275320/35‑1 (MQ=255)
                        tataGTATATCGCGCGCGGCTGGGGGGGGGgtttt            >  1:17833259/1‑34 (MQ=16)
                         ataGTATATCGCGCGCGGCTGTGAGGGTGTTTTCa           <  1:2766655/35‑1 (MQ=255)
                         ataGTATATCGCGCGCGGCTGTGAGGGTGTTTTCa           <  1:16126552/35‑1 (MQ=255)
                          taGTATATCGCGTGCGGCTGTGAGGGTGTTTTCAt          <  1:20148786/35‑1 (MQ=255)
                          taGTATATCGCGCGCGGCTGTGAGGGTGTTTTCAt          <  1:5027022/35‑1 (MQ=255)
                          taGTATATCGCGCGCGGCTGTGAGGGTGTTTTCAt          <  1:4265348/35‑1 (MQ=255)
                          taGTATATCGCGCGCGGCTGTGAGGGTGTTTTCAt          <  1:17936594/35‑1 (MQ=255)
                          taGTATATCGCGCGCGGCTGTGAGGGTGTTTTCAt          <  1:7553111/35‑1 (MQ=255)
                           aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCAtg         <  1:13726626/35‑1 (MQ=255)
                            gTATATCGCGCGCGGCTGTGAGGGTGTTTTCAtgt        <  1:8325799/35‑1 (MQ=255)
                            gTATATCGCGCGCGGCTGTGAGGGTGTTTTCAtgt        <  1:7242424/35‑1 (MQ=255)
                            gTATATCGCGCGCGGCTGTGAGGGTGTTTTCAtgt        <  1:16262335/35‑1 (MQ=255)
                            gTATATCGCGCGCGGCTGTGAGGGTGTTTTCAtgt        <  1:15556587/35‑1 (MQ=255)
                            gTATATCGCGCGCGGCTGTGAGGGTGTTTTCAtgt        <  1:3076675/35‑1 (MQ=255)
                             tataTCGCGCGCGGCTGTGAGGGTGTTTTCAtgtg       <  1:1832694/35‑1 (MQ=255)
                             tataTCGCGCGCGGCTGGGGGGGGGGTTTTTtgtg       >  1:10142363/1‑35 (MQ=9)
                              atatCGCGCGCGGCTGTGAGGGTGTTTTCATgtgt      <  1:7529917/35‑1 (MQ=255)
                              atatCGCGCGCGGCTGTGAGGGTGTTTTCATgtgt      <  1:765620/35‑1 (MQ=255)
                              atatCGCGCGCGGCTGTGAGGGTGTTTTCATgtgt      <  1:16279799/35‑1 (MQ=255)
                               tatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGtt     <  1:8383024/35‑1 (MQ=255)
                                atCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTtc    <  1:11873928/35‑1 (MQ=255)
                                 tCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTtct   <  1:13620606/35‑1 (MQ=255)
                                 tCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTtct   <  1:19003125/35‑1 (MQ=255)
                                 tCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTtct   <  1:10487556/35‑1 (MQ=255)
                                  cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTtctc  <  1:11708368/35‑1 (MQ=255)
                                  |                                  
ACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTC  >  REL606/4335071‑4335139

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