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)77 (1.350) 54/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)69 (1.210) 48/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)

CGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772780‑2772924
                                                                          |                                                                      
cGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAccc                                                                       >  1:1771928/1‑76 (MQ=255)
 gAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAcccc                                                                      <  2:1851356/76‑1 (MQ=255)
      tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGGACACCCTGAACGGGGGGGGAAAAACCCATAAg                                                                 >  2:660057/1‑76 (MQ=37)
      tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                 <  2:531078/76‑1 (MQ=255)
      tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                 <  1:245544/76‑1 (MQ=255)
        tGATTCATCGTGCTGACCTGTTTATCCCCCCCGTCGCGGGGCCCACCCTGAACAGCGGGGACAACACCCCTAAtcg                                                               >  2:985404/1‑76 (MQ=255)
        tGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgcg                                                               <  1:89504/76‑1 (MQ=255)
        tGATTCATCGTGATGACCGGTTTATCCCCGCCGGGCCGGGGAAAACCCTGAACGGGGGGGAAAAACCCCATAAgcc                                                               >  2:985423/1‑75 (MQ=255)
         gATTCATCGTGATGACCGGTTTATCCCCGCCGGCGCGGGGGACAGCCTGGACGGAGGGGGAAAACCCCATAAacgc                                                              >  2:1638891/1‑76 (MQ=25)
             cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                          <  2:866333/76‑1 (MQ=255)
             cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                          <  1:1775228/76‑1 (MQ=255)
              aTCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAACCCATAAGCGCCAACt                                                         >  2:552677/1‑76 (MQ=255)
                 gTGATGCCCGGTTTACCCCCGCGGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                      <  1:1246407/76‑1 (MQ=38)
                 gTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                      <  1:375544/76‑1 (MQ=255)
                   gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    <  1:1774968/76‑1 (MQ=255)
                      gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGGAAAAACCCATAAGCGCTAACTTAAggggt                                                 >  2:1020994/1‑74 (MQ=37)
                       aCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGGGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtg                                                >  2:1601813/1‑76 (MQ=255)
                        ccGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgt                                               <  1:1522186/76‑1 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  2:787921/1‑76 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  2:1608160/1‑76 (MQ=255)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         >  2:234672/1‑76 (MQ=255)
                               ccccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  1:1032532/74‑1 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  2:217689/76‑1 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        >  1:1600877/1‑76 (MQ=255)
                                   ccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                    <  1:1443624/76‑1 (MQ=255)
                                    cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGccn                                   <  1:1233730/76‑2 (MQ=255)
                                    cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                   >  2:1093014/1‑76 (MQ=255)
                                    cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                   <  1:951864/76‑1 (MQ=255)
                                    cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                   <  1:554802/76‑1 (MQ=255)
                                      cTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                 >  1:365609/1‑76 (MQ=255)
                                      cTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                 >  1:1212224/1‑76 (MQ=255)
                                       tGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                <  1:1088619/76‑1 (MQ=255)
                                       tGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                >  1:1134874/1‑76 (MQ=255)
                                         gcgcGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                              <  1:809382/76‑1 (MQ=255)
                                          cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             <  1:772608/76‑1 (MQ=255)
                                          cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             >  1:702863/1‑76 (MQ=255)
                                          cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             <  1:403831/76‑1 (MQ=255)
                                             ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          >  2:1063850/1‑76 (MQ=255)
                                              gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                         <  1:606888/76‑1 (MQ=255)
                                              gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                         >  2:1632040/1‑76 (MQ=255)
                                               ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        >  1:496672/1‑76 (MQ=255)
                                                 aaCCCCCTGTACGGTGGGGAAAGCCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                      <  1:1688462/76‑1 (MQ=17)
                                                   cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    >  2:1467812/1‑76 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                <  1:437884/76‑1 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                <  2:1715271/76‑1 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                <  2:1683749/76‑1 (MQ=255)
                                                        tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa               >  2:783923/1‑76 (MQ=255)
                                                        tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa               <  1:777936/76‑1 (MQ=255)
                                                          aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg             >  2:1788522/1‑76 (MQ=255)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        >  1:1688591/1‑76 (MQ=40)
                                                                ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc       >  1:307417/1‑76 (MQ=39)
                                                                ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc       <  2:1351674/76‑1 (MQ=39)
                                                                   gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    >  2:737590/1‑76 (MQ=34)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  2:871696/76‑1 (MQ=31)
                                                                          |                                                                      
CGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772780‑2772924

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