Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp 94.5% coding (494‑499/549 nt) fimA → major type 1 subunit fimbrin (pilin)

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2772854 =NA (NA)84 (0.930) 36/42 NT 95.2% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 3 (0.050)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 2774196NA (NA)70 (0.780) 36/42 NT 89.3% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 6 (0.090)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

Only 100 of 162 total aligned reads displayed.

GGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTA  >  REL606/2772820‑2772883
                                  |                             
ggCGCGGGGAACACCCTGAACGGTGGGGAAAAAcc                               >  1:11391514/1‑35 (MQ=255)
 gcgcGGGGAACACCCTGAACGGTGGGGAAAAAccc                              >  1:2721636/1‑35 (MQ=255)
 gcgcGGGGAACACCCTGAACGGTGGGGAAAAAccc                              >  1:12745510/1‑35 (MQ=255)
 gcgcGGGGAACACCCTGAACGGTGGGGAAAAAccc                              >  1:1181460/1‑35 (MQ=255)
 gcgcGGGGAACACCCTGAACGGTGGGGAAAAAccc                              <  1:7856019/35‑1 (MQ=255)
  cgcgGGGAACACCCTGAACGGTGGGGAAAAAcccc                             >  1:19209884/1‑35 (MQ=255)
  cgcgGGGAACACCCTGAACGGTGGGGAAAAAcccc                             >  1:3054139/1‑35 (MQ=255)
  cgcgGGGAACACCCTGAACGGTGGGGAAAAAcccc                             <  1:17702286/35‑1 (MQ=255)
  cgcgGGGAACACCCTGAACGGTGGGGAAAAAcccc                             <  1:13421053/35‑1 (MQ=255)
   gcgGGGACCAACCTGAACGGTGGGGAAAAACCCCa                            <  1:2396188/35‑1 (MQ=37)
   gcgGGGAACACCCTGAACGGTGGGGAAAAAcccag                            >  1:196088/1‑33 (MQ=255)
   gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCa                            <  1:18059795/35‑1 (MQ=255)
   gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCa                            >  1:3507069/1‑35 (MQ=255)
   gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCa                            >  1:11615539/1‑35 (MQ=255)
   gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCa                            >  1:10476094/1‑35 (MQ=255)
    cgGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                           >  1:4746353/1‑35 (MQ=255)
    cgGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                           >  1:17159487/1‑35 (MQ=255)
    cgGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                           <  1:7193024/35‑1 (MQ=255)
     ggggAACACCCTGAACGGTGGGGAAAAACCCCATa                          <  1:12590956/35‑1 (MQ=255)
     ggggAACACCCTGAACGGTGGGGAAAAACCCCATa                          >  1:360224/1‑35 (MQ=255)
     ggggAACACCCTGAACGGTGGGGAAAAACCCCATa                          >  1:16351504/1‑35 (MQ=255)
       ggAACACCCTGAACGGTGGGGAAAAACCCCATAAg                        <  1:16183188/35‑1 (MQ=255)
        gAACACCCTGAACGGTGGGGAAAAACCCCATAAgc                       <  1:19952754/35‑1 (MQ=255)
         aaCACCCTGAACGGTGGGGAAAACCCCATAAgcgc                     <  1:3150767/35‑1 (MQ=37)
         aaCACACTGAACGGTGGGGAAAAACCCCATAAgcg                      >  1:7562416/1‑35 (MQ=255)
          acacCCTGAACGGTGGGGAAAAACCCCATAAgcgc                     <  1:5854637/35‑1 (MQ=255)
          acacCCTGAACGGTGGGGAAAAACCCCATAAgcgc                     <  1:6240313/35‑1 (MQ=255)
          acacCCTGAACGGTGGGGAAAAACCCCATAAgcgc                     >  1:12927828/1‑35 (MQ=255)
          acacCCTGAACGGTGGGGAAAAACCCCATAAgcgc                     >  1:1357057/1‑35 (MQ=255)
           cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                    <  1:20196423/35‑1 (MQ=255)
           cacCCTGAACGGTGGGGAAAAACCCCATAAACGCt                    >  1:13866599/1‑35 (MQ=255)
            acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                   <  1:19897573/35‑1 (MQ=255)
             cccTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                  <  1:3454486/35‑1 (MQ=255)
             cccTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                  <  1:10193548/35‑1 (MQ=255)
             cccTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                  <  1:595185/35‑1 (MQ=255)
             cccTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                  <  1:911404/35‑1 (MQ=255)
             cccTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                  <  1:6127792/35‑1 (MQ=255)
             cccTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                  >  1:15232953/1‑35 (MQ=255)
              ncTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                 <  1:6558160/34‑1 (MQ=255)
              ccTGAACTGTGGGGAAAAACCCCATAAGCGATAAc                 >  1:8811814/1‑35 (MQ=25)
              ccTGAACGGTTGGGAAAAACCCCATAAGCGCTAAc                 >  1:934726/1‑35 (MQ=255)
              ccTGAACGGTGGGGAAAAACCCCATAAGCTCTAAc                 >  1:13749086/1‑35 (MQ=255)
              ccTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                 >  1:373127/1‑35 (MQ=255)
              ccTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                 >  1:2857713/1‑35 (MQ=255)
              ccTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                 <  1:39226/35‑1 (MQ=255)
              ccTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                 >  1:17658671/1‑35 (MQ=255)
               cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACt                <  1:8139710/35‑1 (MQ=255)
               cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACt                <  1:237482/35‑1 (MQ=255)
               cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACt                <  1:1364224/35‑1 (MQ=255)
               cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACt                <  1:13183601/35‑1 (MQ=255)
                tGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt               >  1:2522966/1‑35 (MQ=255)
                tGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt               <  1:1311822/35‑1 (MQ=255)
                tGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt               >  1:7254433/1‑35 (MQ=255)
                tGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt               >  1:6926761/1‑35 (MQ=255)
                 gAACGGTGGGTAAAAACCCCATAAGCGCTAACTTa              >  1:3795501/1‑35 (MQ=255)
                 gAACGGTGGGGGAAAACCCCATAAGCGCTAACTTa              <  1:4337295/35‑1 (MQ=255)
                 gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa              >  1:2158681/1‑35 (MQ=255)
                 gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa              >  1:1575261/1‑35 (MQ=255)
                 gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa              <  1:19565655/35‑1 (MQ=255)
                  aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTaa             >  1:15585599/1‑35 (MQ=255)
                  aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTaa             <  1:19432573/35‑1 (MQ=255)
                  aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTaa             <  1:9583602/35‑1 (MQ=255)
                  aaCGGTGGGGAAAAACCCCATAAGAGCTAACTTaa             <  1:15926668/35‑1 (MQ=255)
                   aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAg            <  1:2035563/35‑1 (MQ=255)
                   aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAg            <  1:12374056/35‑1 (MQ=255)
                    cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg           >  1:5106204/1‑35 (MQ=255)
                    cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg           >  1:968714/1‑35 (MQ=255)
                    cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg           >  1:11973866/1‑35 (MQ=255)
                    cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg           <  1:1047088/35‑1 (MQ=255)
                     ggTGGGGAAAAACCCCATAAGCGCTAACTTAAggg          <  1:14126372/35‑1 (MQ=255)
                     ggTGGGGAAAAACCCCATAAGCGCTAACTTAAggg          >  1:14485136/1‑35 (MQ=255)
                     ggTGGGGAAAAACCCCATAAGCGCTAACTTAAggg          <  1:15547320/35‑1 (MQ=255)
                     ggTGGGGAAAAACCCCATAAGCGCTAACTTAAggg          >  1:10666634/1‑35 (MQ=255)
                      gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt         <  1:11429923/35‑1 (MQ=255)
                       tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt        >  1:19196254/1‑35 (MQ=255)
                       tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt        <  1:9472935/35‑1 (MQ=255)
                        ggggAAAAACCCCATAAGCGCTAACTTAAGGGTtg       <  1:8956749/35‑1 (MQ=255)
                        ggggAAAAACCCCATAAGCGCTAACTTAAGGGTtg       >  1:17312594/1‑35 (MQ=255)
                        ggggAAAAACCCCATAAGCGCTAACTTAAGGGTtg       <  1:11493048/35‑1 (MQ=255)
                        ggggAAAAACCCCATAAGCGCTAACTTAAGGGTtg       >  1:4020473/1‑35 (MQ=255)
                         gggAAAAACCCCATAAGCGCTAACTTAAGGGTtgt      >  1:20059/1‑35 (MQ=255)
                         gggAAAAACCCCATAAGCGCTAACTTAAGGGTtgt      <  1:12463219/35‑1 (MQ=255)
                         gggAAAAACCCCATAAGCGCTAACTTAAGGGTtgt      <  1:20463792/35‑1 (MQ=255)
                         gggAAAAACCCCATAAGCGCTAACTTAAGGGTtgt      >  1:11195525/1‑35 (MQ=255)
                         gggAAAAACCCCATAAGCGCTAACTTAAGGGTtgt      >  1:3736820/1‑35 (MQ=255)
                         gggAAAAACCCCATAAGCGCTAACTTAAGGGTtgt      <  1:13158807/35‑1 (MQ=255)
                          ggAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg     >  1:4116664/1‑35 (MQ=40)
                          ggAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg     <  1:6929894/35‑1 (MQ=40)
                          ggAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg     >  1:2269212/1‑35 (MQ=40)
                           gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg    <  1:15205727/35‑1 (MQ=38)
                           gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg    >  1:9178577/1‑35 (MQ=38)
                           gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg    <  1:18910403/35‑1 (MQ=38)
                            aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGt   >  1:10170114/1‑35 (MQ=35)
                            aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGt   >  1:17264292/1‑35 (MQ=35)
                            aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGt   <  1:5348286/35‑1 (MQ=35)
                            aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGt   <  1:18277986/35‑1 (MQ=35)
                             aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTa  <  1:7053081/35‑1 (MQ=32)
                             aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTa  >  1:14343118/1‑35 (MQ=32)
                             aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTa  >  1:5896189/1‑35 (MQ=32)
                             aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTa  >  1:2896081/1‑35 (MQ=32)
                                  |                             
GGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTA  >  REL606/2772820‑2772883

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