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)45 (0.830) 38/144 0.3 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)70 (1.290) 52/144 0.1 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

ATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772784‑2772924
                                                                      |                                                                      
aTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGGACACCCTGAACGGGGGGGGAAAACCCCATa                                                                   >  1:543783/1‑76 (MQ=38)
 ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAGCCCATaa                                                                  <  1:1178649/76‑1 (MQ=255)
 ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                  <  1:1494382/76‑1 (MQ=255)
 ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACACCATaa                                                                  <  1:1481949/76‑1 (MQ=255)
  tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                 <  1:490840/76‑1 (MQ=255)
      aTTCTTCTTGAAGTCCGCTTAACCCCGGCTGGAGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                                                             <  1:831916/76‑1 (MQ=255)
       ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                            <  1:1515557/76‑1 (MQ=255)
           tCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        <  1:1657687/76‑1 (MQ=255)
             gTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                      >  2:1701876/1‑76 (MQ=255)
                       gTTTATCCCCGCTGGCGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                            <  1:1114540/76‑1 (MQ=255)
                         ttATCCCCGCTGGCGCGGGGAAAACCCTGAACGGTGGGGTAAAACCCCATAAGCGCTAACTTAAGGGGTGTGGTAt                                          >  2:575160/1‑76 (MQ=37)
                          tcgccccGCGGGGGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTTTGGTAtt                                         <  1:660904/73‑1 (MQ=25)
                          tacccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:669066/73‑1 (MQ=255)
                          tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:468184/76‑1 (MQ=255)
                          tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:178651/76‑1 (MQ=255)
                          tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:1521369/76‑1 (MQ=255)
                          tATCCCCGCGGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:1221417/76‑1 (MQ=255)
                           aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        >  2:130364/1‑76 (MQ=255)
                           aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  1:29334/76‑1 (MQ=255)
                            tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGGAAAACCCCATAAGCGCTAACTTAAGGGGTGTGGTATTAc                                       >  2:628582/1‑76 (MQ=255)
                            tCCCCGCCCGCGCGGGGTACAACCTTCACCGAGGGGGAAAAACCCACAAGCGCTAACCTAAAGGTTGATGGATTAc                                       >  2:1321145/1‑76 (MQ=255)
                             ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      <  2:1386379/76‑1 (MQ=255)
                                cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                   <  2:391951/76‑1 (MQ=255)
                                  cTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                 >  2:467513/1‑76 (MQ=255)
                                   tGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                <  2:496090/76‑1 (MQ=255)
                                   cggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                <  1:1153026/75‑1 (MQ=255)
                                   aagcggGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTGGTGGTATTACGCCTGat                                <  2:1082393/74‑1 (MQ=21)
                                       gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                            <  1:239673/76‑1 (MQ=255)
                                        cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                           >  1:1710234/1‑76 (MQ=255)
                                         ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          <  1:668807/76‑1 (MQ=255)
                                         ggggAACACCCCGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTNCGCCTGATATGAtt                          <  1:1001164/76‑1 (MQ=255)
                                           ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGNCTGATATGATTTa                        <  1:1231613/76‑1 (MQ=255)
                                           ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        <  2:128792/76‑1 (MQ=255)
                                            gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       >  1:687963/1‑76 (MQ=255)
                                             gacacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                      <  1:1454263/75‑1 (MQ=255)
                                             aaCACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCNGATATGATTTAAc                      <  1:1233766/76‑1 (MQ=255)
                                                acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  1:946146/76‑1 (MQ=255)
                                                acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  1:391070/76‑1 (MQ=255)
                                                acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  1:1195471/76‑1 (MQ=255)
                                                   cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                >  1:1528296/1‑76 (MQ=255)
                                                    tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa               <  1:1050379/76‑1 (MQ=255)
                                                      gaCGGTGGGGAAAATCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg             <  1:34697/75‑1 (MQ=255)
                                                      aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg             >  2:299601/1‑76 (MQ=255)
                                                       aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa            <  1:315149/76‑1 (MQ=255)
                                                       aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa            <  1:1160514/76‑1 (MQ=255)
                                                         ggTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt          >  1:790198/1‑76 (MQ=255)
                                                          gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAATGTGCCGATGAAtt         <  1:1221966/76‑1 (MQ=255)
                                                          gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         <  2:1492009/76‑1 (MQ=255)
                                                           tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAttn        <  1:1449337/76‑2 (MQ=255)
                                                           tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        <  1:1442919/76‑1 (MQ=40)
                                                           tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        <  1:1315601/76‑1 (MQ=40)
                                                           tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        <  2:1756702/76‑1 (MQ=40)
                                                           tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        >  2:219751/1‑76 (MQ=40)
                                                              ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc     >  1:1037658/1‑76 (MQ=35)
                                                               gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  2:480500/76‑1 (MQ=34)
                                                               gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  1:1273407/76‑1 (MQ=34)
                                                               gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  2:1245624/76‑1 (MQ=34)
                                                                 aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  2:1356979/76‑1 (MQ=31)
                                                                      |                                                                      
ATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772784‑2772924

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