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)66 (1.300) 53/144 0.0 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)47 (0.930) 37/144 0.4 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

ATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772784‑2772924
                                                                        |                                                                      
aTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATa                                                                   <  2:1083695/76‑1 (MQ=255)
  tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAg                                                                 <  2:567215/76‑1 (MQ=255)
  tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAg                                                                 <  2:567197/76‑1 (MQ=255)
    tGATTCATCGTGATGACCGGTTTATCCCCCCCGGCGCGGGGGACACCCCGAAC‑‑GGGGGGGGAAAAACCCCTAAacg                                                               >  2:1050676/1‑76 (MQ=255)
      aTTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGGAAAACCCCATAAGCGCt                                                             >  2:350531/1‑76 (MQ=255)
       ttCATCGTGATGACCGGTTTATCCCCCCCGGCGCGGGGGGGACCCCCAACGGGGGGGGGAAAACCCCCAAAAgcta                                                              >  1:332561/1‑74 (MQ=255)
         cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAAc                                                          <  1:1398878/76‑1 (MQ=255)
         cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAAc                                                          <  1:1411776/76‑1 (MQ=255)
           tCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        <  2:1596340/76‑1 (MQ=255)
              tgatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAg                                                     <  1:1494873/76‑1 (MQ=255)
               gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    <  1:333920/76‑1 (MQ=255)
                atgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAggg                                                   >  2:1549199/1‑76 (MQ=255)
                  gaCCGGTTTCCCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                 <  1:613210/76‑1 (MQ=255)
                  gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGGGGGGAAAAAACCCATAAGCGCTAACTTAAgggct                                                 >  2:500445/1‑74 (MQ=37)
                       gTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTACGGGTTGTGGt                                            >  2:176690/1‑76 (MQ=255)
                         ttATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  2:177841/1‑76 (MQ=255)
                         ttACCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          <  1:1274957/76‑1 (MQ=255)
                          tATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         >  1:1575237/1‑76 (MQ=255)
                          tATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  2:862273/76‑1 (MQ=255)
                            tCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTCACTTAAGGGTTGTGGGATTaa                                       >  2:253142/1‑75 (MQ=255)
                            tCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  1:878858/76‑1 (MQ=255)
                              cccGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                     <  1:1291290/76‑1 (MQ=255)
                                cGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                   <  2:454042/76‑1 (MQ=255)
                                cGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                   <  1:393199/76‑1 (MQ=255)
                                   tGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                >  1:1593236/1‑76 (MQ=255)
                                      cgcgGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             <  1:732198/76‑1 (MQ=255)
                                      cgcgGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             <  1:438382/76‑1 (MQ=255)
                                        cgGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                           <  1:1521870/76‑1 (MQ=255)
                                        cgGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGGTGTGGTATTACGCCTGATATGAt                           >  2:1331239/1‑76 (MQ=255)
                                         ggggAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          <  1:1104572/76‑1 (MQ=255)
                                          gggAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                         <  2:1575237/76‑1 (MQ=255)
                                           ggAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        <  2:446217/76‑1 (MQ=255)
                                            gAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       <  1:505699/76‑1 (MQ=255)
                                               cacCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    >  2:1374362/1‑76 (MQ=255)
                                               cacCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    <  1:1430471/76‑1 (MQ=255)
                                                   cTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                <  2:16226/76‑1 (MQ=255)
                                                     gAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt              <  1:1246622/76‑1 (MQ=255)
                                                      aaC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg             >  2:859/1‑76 (MQ=255)
                                                       aC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATgn            <  1:1524927/76‑2 (MQ=255)
                                                       aC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa            >  1:919724/1‑76 (MQ=255)
                                                            gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         >  2:207945/1‑76 (MQ=255)
                                                             tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        >  2:866194/1‑76 (MQ=40)
                                                             tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        <  2:1440391/76‑1 (MQ=40)
                                                             tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        >  2:1323898/1‑76 (MQ=40)
                                                              ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc       <  1:988496/76‑1 (MQ=39)
                                                                ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc     <  1:1617086/76‑1 (MQ=35)
                                                                   aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  1:1447557/76‑1 (MQ=31)
                                                                        |                                                                      
ATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772784‑2772924

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