Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp 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)61 (0.890) 43/144 0.5 100% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 0 (0.000)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 2774196NA (NA)55 (0.810) 41/144 0.7 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

CGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTC  >  REL606/2772780‑2772923
                                                                           |                                                                     
cGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGGACACCCCGAACGGTGGGGGAAAAccc                                                                      >  1:349182/1‑76 (MQ=37)
cGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAccc                                                                      <  2:633004/76‑1 (MQ=255)
cGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAccc                                                                      <  2:1554391/76‑1 (MQ=255)
   aaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                                                                   <  2:1600338/76‑1 (MQ=255)
    aTTTTGATTCATCGTGATGACCGGTTTATCCCCCCCGGCGCGGGGGAAAACCCGCACAGGGGGGGAAAAACCCCaa                                                                   >  1:1025164/1‑75 (MQ=255)
     ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                 <  1:517633/76‑1 (MQ=255)
     ttttGATTCATCGTGATGACCGGTTTATCCCCCCCGGCGCGGGGGAAAACCCCGACGGGGGGGGGAAAAcccccaa                                                                  >  1:1782835/1‑73 (MQ=255)
      tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                <  1:1370936/76‑1 (MQ=255)
        tGATTCATCGTGATGACCGGTTTATCCCCGCCGGCGCGGGGAACACCCTGAACGGGGGGGAAAAACCCCATAAgcg                                                              >  1:2187630/1‑76 (MQ=255)
         gATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgcgc                                                             <  1:1143875/76‑1 (MQ=255)
          aTTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                                                            >  2:177866/1‑76 (MQ=255)
           ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                           <  2:1115925/76‑1 (MQ=255)
            tCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                                                          <  2:1240262/76‑1 (MQ=255)
             cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                         <  1:1806856/76‑1 (MQ=255)
               tcctgatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                       <  1:1599689/73‑1 (MQ=255)
                cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGGGGGGAAAAACCCCATAAGCGCTAACTTa                                                      >  2:2106511/1‑76 (MQ=255)
                 gTGATGCCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                     <  1:954898/76‑1 (MQ=255)
                  tgatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGGGGGGAAAAACCCCATAAGCGCTAACTTAAg                                                    <  1:834540/76‑1 (MQ=255)
                    atgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAggg                                                  >  2:201228/1‑76 (MQ=255)
                      gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                <  1:883111/76‑1 (MQ=255)
                      gaCCGGTTTATCCCCGCCGGCGCGGGGAACACCCTGAACGGAGGGGAAAAATCCCATAAGCGCTAACTTAAGGGtt                                                >  2:2180583/1‑76 (MQ=37)
                       aCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGGtg                                               >  2:1900505/1‑76 (MQ=255)
                        ccGGTTTATCCCCGCCTGCGCGGGGGACACCCTGAACCGCGGGGAAAAAACCCATAAGCGCTAACCTAAAGGgtgt                                              >  2:1976783/1‑76 (MQ=21)
                         cGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg                                             >  1:686525/1‑76 (MQ=255)
                         cGGTTTATCCCCCCCGGCGCGGGGGACACCCCGAACGGGGGGGGAAAAACCCCCAAGCGGTAACTTAAGGGGtggg                                             >  2:1271837/1‑74 (MQ=16)
                          ggTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                            <  2:1151817/76‑1 (MQ=255)
                           gTTTACCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                           <  1:1211964/76‑1 (MQ=255)
                            tttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                          <  1:977002/76‑1 (MQ=255)
                            tttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                          <  1:1595895/76‑1 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGNGCTAACTTAAGGGTTGTGGTAt                                         <  2:300457/76‑1 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                         <  1:315322/76‑1 (MQ=255)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                        >  2:1278177/1‑76 (MQ=255)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                        >  2:2174595/1‑76 (MQ=255)
                               ccccccGCGGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                       <  1:1051060/74‑1 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                       <  2:682552/76‑1 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                       <  2:1304863/76‑1 (MQ=255)
                                tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                      <  2:415473/76‑1 (MQ=255)
                                 ccccGTTGGGGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                     <  1:657923/76‑1 (MQ=255)
                                 ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                     >  2:1925724/1‑76 (MQ=255)
                                  cccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                    >  2:2142924/1‑76 (MQ=255)
                                   ccGCGGGCGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                   <  1:1262901/76‑1 (MQ=255)
                                    cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                  >  1:830972/1‑76 (MQ=255)
                                     gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                 <  2:1658396/76‑1 (MQ=255)
                                      cTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                <  1:1394521/76‑1 (MQ=255)
                                        ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                              >  1:785005/1‑76 (MQ=255)
                                         gcgcGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                             <  2:1938076/76‑1 (MQ=255)
                                          cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                            <  1:264349/76‑1 (MQ=255)
                                              ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                         <  1:1858321/76‑1 (MQ=255)
                                               gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                        <  2:904698/76‑1 (MQ=255)
                                               gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                        <  2:540920/76‑1 (MQ=255)
                                                ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                       <  2:1778290/76‑1 (MQ=255)
                                                ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                       <  1:1175456/76‑1 (MQ=255)
                                                 gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                      >  2:421157/1‑76 (MQ=255)
                                                        cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg               <  2:460795/76‑1 (MQ=255)
                                                         tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa              <  1:229875/76‑1 (MQ=255)
                                                         tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa              >  1:2019627/1‑76 (MQ=255)
                                                         tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa              >  2:1578134/1‑76 (MQ=255)
                                                          gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt             <  1:1704312/76‑1 (MQ=255)
                                                          gAACAGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt             >  1:1389250/1‑76 (MQ=255)
                                                             cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa          <  1:2131977/76‑1 (MQ=255)
                                                             cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa          <  2:693430/76‑1 (MQ=255)
                                                             cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATAATTTAACGTGCCGATGaa          <  1:474659/76‑1 (MQ=255)
                                                               gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt        <  1:386396/76‑1 (MQ=255)
                                                                 ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc      >  2:1383973/1‑76 (MQ=39)
                                                                 ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc      >  1:343251/1‑76 (MQ=39)
                                                                   ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc    <  2:1542317/76‑1 (MQ=35)
                                                                     aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc  <  1:1068153/76‑1 (MQ=32)
                                                                     aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc  <  1:609134/76‑1 (MQ=32)
                                                                           |                                                                     
CGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTC  >  REL606/2772780‑2772923

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