Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,183,348 G→A 32.6% G80S (GGC→AGC)  yqjK → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,183,3480G→A32.6% 70.7 ‑0.6 ‑0.189G80S (GGC→AGC) yqjKhypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/25):  ref base (12/48):  total (16/73)
Fisher's exact test strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.62e-01

Only 100 of 106 total aligned reads displayed.

TAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACGAC  >  REL606/3183314‑3183382
                                    |                                  
tAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTg                                    <  1:10745931/35‑1 (MQ=255)
tAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTg                                    <  1:2015684/35‑1 (MQ=255)
tAATATGCTGGTCCGCTGGGCAAGACGCGGTTTTg                                    <  1:4825314/35‑1 (MQ=255)
 aaTATGCTGGTCCGCTGGGCCAGACGCGGTTTTgg                                   <  1:159491/35‑1 (MQ=255)
  gtatGCTGGTCCGCTGGGCCAGACGCGGTTTTGGc                                  <  1:14270836/34‑1 (MQ=255)
   tatGCTGGTCCGCTGGGCCAGACGCGGGTTTTggg                                  >  1:14621058/1‑34 (MQ=21)
   tatGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCg                                 <  1:8350863/35‑1 (MQ=255)
   tatGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCg                                 <  1:12537868/35‑1 (MQ=255)
   tatGCTGGTCCGCTGGGCCAGACGCGGTTTTAGCg                                 <  1:15507068/35‑1 (MQ=255)
    gtGCTGGTCCGCTGGGCCAGACGCGGGTTTTGCgg                                >  1:1116865/2‑34 (MQ=21)
    atGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGt                                <  1:7445360/35‑1 (MQ=255)
    atGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGt                                <  1:1249704/35‑1 (MQ=255)
    atGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGt                                <  1:14070393/35‑1 (MQ=255)
    atGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGt                                <  1:6965298/35‑1 (MQ=255)
     tGCTGGTCCGCTGGGCCAGACGCGGCTATTGGCGt                                >  1:6909050/1‑35 (MQ=18)
     tGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTa                               <  1:12106068/35‑1 (MQ=255)
     tGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTa                               <  1:3439423/35‑1 (MQ=255)
     tGCTGGTCCGCTGGGCCAGACGCGGGTTTTGCgga                               >  1:11697890/1‑33 (MQ=21)
      gCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTAt                              <  1:13388297/35‑1 (MQ=255)
       cTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATg                             <  1:20240096/35‑1 (MQ=255)
       cTGGTCCGCTGGGCCAGACGCGGTTTTAGCGTATg                             <  1:19770096/35‑1 (MQ=255)
       cTGGTCCGCTGGGCCAGACGCGGGGTTTGGGTATg                             >  1:17714816/1‑35 (MQ=16)
        tGGTCCGCTGGGCCAGACGCGGGTTTTGGGGttgg                             >  1:20127238/1‑32 (MQ=16)
        tGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATgg                            <  1:17830352/35‑1 (MQ=255)
        tGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATgg                            <  1:17782168/35‑1 (MQ=255)
        tGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATgg                            <  1:16496004/35‑1 (MQ=255)
        tGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATgg                            <  1:16003448/35‑1 (MQ=255)
        tGGTCCGCTGGGCCAGACGCGGGTTTTGGGTGTgg                            >  1:854996/1‑35 (MQ=16)
         ggTCAGCTGGGCCAGACGCGGTTTTAGCGTATGGa                           <  1:5635186/35‑1 (MQ=37)
          gTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAg                          <  1:11348/35‑1 (MQ=255)
           tcagCTGGGCCAGACGCGGTTTTGGCGTATGGAgc                         <  1:16900796/32‑1 (MQ=255)
           tCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAgc                         <  1:7773906/35‑1 (MQ=255)
           tCCGCTGGGCCAGACGCGGTTTTAGCGTATGGAgc                         <  1:19110097/35‑1 (MQ=255)
           tCCGCTGGGCCAGACGCGGTTTTAGCGTATGGAgc                         <  1:5996519/35‑1 (MQ=255)
           tCCGCTGGGCCAGACGCGGTTTTAGCGTATGGAgc                         <  1:18072449/35‑1 (MQ=255)
           tCCGCTGGGCCAGACGCGGGTTTTGGGGATggggg                         >  1:16257619/1‑32 (MQ=12)
            ccGCTGGGCCAGACGCGGTTTTGGCGTATGGAgcg                        <  1:11897593/35‑1 (MQ=255)
             cGCTGGGCCAGACGCGGTTTTGGCGTATGGAgcgc                       <  1:1576323/35‑1 (MQ=255)
             cGCTGGGCCAGACGCGGGTTTTGCGTATGGGgcgc                       >  1:1987806/1‑35 (MQ=25)
              gCTGGGCCAGACGCGGTTTTGGCGTATGGAGCGcc                      <  1:5618176/35‑1 (MQ=255)
              gCTGGGCCAGACGCGGTTTTAGCGTATGGAGCGcc                      <  1:16365913/35‑1 (MQ=255)
               cTGGGCCAGACGCGGTTTTGGCGTATGGAGCGCCt                     <  1:18178555/35‑1 (MQ=255)
               cTGGGCCAGACGCGGTTTTAGCGTATGGAGCGCCt                     <  1:19228166/35‑1 (MQ=255)
                tGGGCCAGACGCGGTTTTGGCGTATGGAGCGCCTg                    <  1:4553197/35‑1 (MQ=255)
                tGGGCCAGACGCGGTTTTAGCGTATGGAGCGCCTg                    <  1:6895846/35‑1 (MQ=255)
                 ggtccAGACGCGGTTTTGGCGTATGGAGCGCCTgg                   <  1:11707846/32‑1 (MQ=255)
                 gggCGAGACGCGGTTTTGGCGTATGGAGCGCCTgg                   <  1:3260692/35‑1 (MQ=255)
                 gggCCAGACGCGGTTTTGGCGTATGGAGCGCCTgg                   <  1:20352675/35‑1 (MQ=255)
                 gggCCAGACGCGGTTTTAGCGTATGGAGCGCCTgg                   <  1:10860669/35‑1 (MQ=255)
                  ggCCAGACGCGGTTTTGGCGTATGGAGCGCCTGGc                  <  1:10660709/35‑1 (MQ=255)
                    ccAGACGCGGTTTTGGCGTATGGAGCGCCTGGCGt                <  1:898166/35‑1 (MQ=255)
                    ccAGACGCGGTTTTGGCGTATGGAGCGCCTGGCGt                <  1:20152792/35‑1 (MQ=255)
                    ccAGACGCGGTTTTAGCGTATGGAGCGCCTGGCGt                <  1:20371362/35‑1 (MQ=255)
                    ccAGACGCGGTTTTAGCGTATGGAGCGCCTGGCGt                <  1:2480825/35‑1 (MQ=255)
                    ccAGACGCGGTTTTAGCGTATGGAGCGCCTGGCGt                <  1:3869624/35‑1 (MQ=255)
                     cAGACGCGGTTTTGGCGTATGGAGCGCCTGGCGTc               <  1:8958750/35‑1 (MQ=255)
                      aGACGCGGTTTTAGCGTATGGAGCGCCTGGCGTCt              <  1:10858888/35‑1 (MQ=255)
                       gACGCGGTTTTGGCGTATGGAGCGCCTGGCGTCTg             <  1:18486524/35‑1 (MQ=255)
                       gACGCGGTTTTGGCGTATGGAGCGCCTGGCGTCTg             <  1:15508924/35‑1 (MQ=255)
                       gACGCGGTTTTGGCGTATGGAGCGCCTGGCGTCTg             <  1:20588376/35‑1 (MQ=255)
                       gACGCGGTTTTAGCGTATGGAGCGCCTGGCGTCTg             <  1:19712612/35‑1 (MQ=255)
                       gACGCGGTTTTAGCGTATGGAGCGCCTGGCGTCTg             <  1:7417235/35‑1 (MQ=255)
                        gcgcgGTTTTAGCGTATGGAGCGCCTGGCGTCTgg            <  1:18269675/34‑1 (MQ=255)
                        aCGCGGTTTTGGCGTATGGAGCGCCTGGCGTCTgg            >  1:19870545/1‑35 (MQ=255)
                        aCGCGGTTTTAGCGTATGGAGCGCCTGGCGTCTgg            <  1:8625080/35‑1 (MQ=255)
                        aCGCGGTTTTAGCGTATGGAGCGCCTGGCGTCTgg            <  1:9461796/35‑1 (MQ=255)
                        aCGCGGTTTTAGCGTATGGAGCGCCTGGCGTCTgg            >  1:1416461/1‑35 (MQ=255)
                         cgcgGTTTTGGCGTATGGAGCGCCTGGCGTCTGGt           <  1:3400389/35‑1 (MQ=255)
                         cgcgGTTTTGGCGTATGGAGCGCCTGGCGTCTGGt           >  1:17637056/1‑35 (MQ=255)
                         cgcgGTTTTGGCGTATGGAGCGCCTGGCGTCTGGt           <  1:8685938/35‑1 (MQ=255)
                         cgcgGTTTTGGCGTATGGAGCGCCTGGCGTCTGGt           <  1:12323617/35‑1 (MQ=255)
                         cgcgGTTTTAGCGTATGGAGCGCCTGGCGTCTGGt           <  1:5250198/35‑1 (MQ=255)
                          gcgGTTTTGGCGTATGGAGCGCCTGGCGTCTGGtt          >  1:1762714/1‑35 (MQ=255)
                             tgTTTTAGCGTATGGAGCGCCTGGCGTCTGGTTaa        <  1:1350393/34‑1 (MQ=255)
                             ggTTTTGGCGTATGGAGCGCCTGGCGTCTGGTTaa        >  1:3705861/1‑35 (MQ=255)
                             ggTTTTAGCGTATGGAGCGCCTGGCGTCTGGTTaa        <  1:14551050/35‑1 (MQ=255)
                             ggTTTTAGCGTATGGAGCGCCTGGCGTCTGGTTaa        >  1:19696644/1‑35 (MQ=255)
                             ggTTTTAGCGTATGGAGCGCCTGGCGTCTGGTTaa        <  1:16967565/35‑1 (MQ=255)
                             ggTTTTAGCGTATGGAGCGCCTGGCGTCTGGTTaa        >  1:5645391/1‑35 (MQ=255)
                              gTTTTGGCGTATGGAGCGCCTGGCGTCTGGTTaaa       >  1:3927372/1‑35 (MQ=255)
                                ttttGGCGTATGGAGCGCCTGGCGTCTGGTTaaaa      <  1:3323008/35‑1 (MQ=255)
                                ttttGGCGTATGGAGCGCATGGCGTCTGGTTaaaa      <  1:6480825/35‑1 (MQ=255)
                                 tttGGCGTATGGAGCGCCTGGCGTCTGGTTAAAAc     <  1:11719595/35‑1 (MQ=255)
                                 tttGGCGTATGGAGCGCCTGGCGTCTGGTTAAAAc     <  1:10470224/35‑1 (MQ=255)
                                 tttGGCGTATGGAGCGCCTGGCGTCTGGTTAAAAc     >  1:10026178/1‑35 (MQ=255)
                                 tttAGCGTATGGAGCGCCTGGCGTCTGGTTAAAAc     <  1:16718047/35‑1 (MQ=255)
                                 tttAGCGTATGGAGCGCCTGGCGTCTGGTTAAAAc     >  1:8265638/1‑35 (MQ=255)
                                 tttAGCGTATGGAGCGCCTGGCGTCTGGTTAAAAc     <  1:20250740/35‑1 (MQ=255)
                                  ttGGCGTATGGAGCGCCTGGCGTCTGGTTaaaaaa    >  1:13125987/1‑33 (MQ=255)
                                  ttGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACg    <  1:1400570/35‑1 (MQ=255)
                                  ttGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACg    <  1:7838392/35‑1 (MQ=255)
                                  ttGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACg    >  1:12938509/1‑35 (MQ=255)
                                  ttGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACg    <  1:4895497/35‑1 (MQ=255)
                                  ttGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACg    <  1:3521028/35‑1 (MQ=255)
                                   tGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACGa   >  1:15565964/1‑35 (MQ=255)
                                   tGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACGa   >  1:7814721/1‑35 (MQ=255)
                                   tGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACGa   <  1:5161618/35‑1 (MQ=255)
                                   tGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACGa   >  1:123184/1‑35 (MQ=255)
                                    ggCGTATGGAGCGCCTGGCGTCTGGTTAAAACGAc  <  1:14168020/35‑1 (MQ=255)
                                    ggCGTATGGAGCGCCTGGCGTCTGGTTAAAACGAc  <  1:17456969/35‑1 (MQ=255)
                                    |                                  
TAATATGCTGGTCCGCTGGGCCAGACGCGGTTTTGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACGAC  >  REL606/3183314‑3183382

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.