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)46 (0.860) 34/144 0.6 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)41 (0.770) 35/144 0.5 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

ATTTTGATTCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772784‑2772924
                                                                        |                                                                      
aTTTTGATTCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATa                                                                   <  1:914104/76‑1 (MQ=255)
aTTTTGATTCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATa                                                                   <  2:1227737/76‑1 (MQ=255)
  tttGTTTCATCGTGATGCCCGGTTTAC‑‑CCCCGCGGGCGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                 <  2:732307/76‑1 (MQ=255)
  tttGATTCATCGTGATGACCGGTTTATTCCCCCCCGGGCGGGGGGCAAACCCAAAGGGGGGGGGAAAAACCCCaaa                                                                   >  1:1051495/1‑74 (MQ=255)
  tttGATTCATCGTGATGACCGGTTTAT‑‑CCCCGCCGGCGCGGGGGACACCCTGAACGGGGGGGGAAAAACCCCTAAg                                                                 >  2:1084658/1‑76 (MQ=25)
    tGATTCAGCGTGAGGACCGGTTTAC‑‑CCCCGCGGGCGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgcg                                                               <  1:666800/76‑1 (MQ=255)
      aTTCATCGTGATGCCCGTTTTAC‑‑CCCCGCTGGCGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                                                             <  1:704256/76‑1 (MQ=255)
      aTTCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                                                             <  2:369952/76‑1 (MQ=255)
       ttCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                            <  2:1312544/76‑1 (MQ=255)
       ttCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                            <  1:819441/76‑1 (MQ=255)
       tccATCGTGATGCCCGGTTTAC‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                            <  1:357686/74‑1 (MQ=38)
          aTCGTGATGACCGGTTTAT‑‑CCCCCCTGGCGGGGGGGACACCCCGAACGGGGGGGGAAAACCCCATACACGCTAACt                                                         >  2:1440735/1‑76 (MQ=255)
           tCGTGATGCCCGTTTTCC‑‑CCCCGCTGGCGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        <  1:820194/76‑1 (MQ=38)
           tCGTGATGACCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        <  1:1592137/76‑1 (MQ=255)
               gatgaCCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    <  1:319969/76‑1 (MQ=255)
                 tgaCCGGTTTAT‑‑CCCCGCCGGCGGGGGGGAGACCCTGAACGGCGGGGGGAAAACCCCTAAGCGCGTACTTTAgggg                                                  >  2:1308395/1‑75 (MQ=255)
                  gaCCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                 <  1:1483188/76‑1 (MQ=255)
                  gaCCGGTTTAC‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                 <  1:771436/76‑1 (MQ=255)
                        tttAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATCAGCGCTAACTTAAGGGTTGTGGTa                                           >  2:1169449/1‑76 (MQ=255)
                        tttAC‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                           <  1:407902/76‑1 (MQ=255)
                         ttAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  2:248703/1‑76 (MQ=255)
                          tAT‑‑CCCCGCTGGCGCGGGGGACACCCTGAACGGCGGGGAAAAACCCCCACAGCGCCAACTTAAGGGGTATGGGAtt                                         >  2:439931/1‑76 (MQ=21)
                          tAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:641640/76‑1 (MQ=255)
                          tAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:1754898/76‑1 (MQ=255)
                           aT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        >  2:976023/1‑76 (MQ=255)
                            t‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  2:800571/76‑1 (MQ=255)
                                   gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  <  1:1585093/76‑1 (MQ=255)
                                        cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             <  1:467435/76‑1 (MQ=255)
                                          cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                           <  1:226400/76‑1 (MQ=255)
                                           ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          >  2:316340/1‑76 (MQ=255)
                                            gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                         <  1:1430646/76‑1 (MQ=255)
                                             ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        <  1:503087/76‑1 (MQ=255)
                                             ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        <  1:826342/76‑1 (MQ=255)
                                              gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       <  1:710666/76‑1 (MQ=255)
                                              gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       <  1:286787/76‑1 (MQ=255)
                                                 cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    <  1:535674/76‑1 (MQ=255)
                                                 cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    <  1:402390/76‑1 (MQ=255)
                                                 cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    >  1:1016308/1‑76 (MQ=255)
                                                  acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  2:720/76‑1 (MQ=255)
                                                  acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   >  1:622579/1‑76 (MQ=255)
                                                    ccTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                 >  1:1270394/1‑76 (MQ=255)
                                                      tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa               >  1:1540619/1‑76 (MQ=255)
                                                        aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg             <  2:418226/76‑1 (MQ=255)
                                                         aCGGTCGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa            <  1:328115/76‑1 (MQ=255)
                                                            gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         >  1:939324/1‑76 (MQ=255)
                                                            gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         <  1:1308395/76‑1 (MQ=255)
                                                            gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         >  1:1239153/1‑76 (MQ=255)
                                                            gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         <  1:404643/76‑1 (MQ=255)
                                                                ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc     >  1:608902/1‑76 (MQ=35)
                                                                 gAAAAACCCCATAAGCGCTAACTTATGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  1:182999/76‑1 (MQ=21)
                                                                   aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  >  2:450445/1‑76 (MQ=31)
                                                                        |                                                                      
ATTTTGATTCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772784‑2772924

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