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)43 (1.000) 34/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)38 (0.880) 31/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)

GAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTC  >  REL606/2772781‑2772923
                                                                          |                                                                     
gAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAcccc                                                                     <  1:871613/76‑1 (MQ=255)
  aaTTTTGACTCATCGTGATGACCGGTTTATCCCCGCTGGGGCGGGGAACACCCTGAACGGTGGGGATAAACCCCAt                                                                   >  1:42063/1‑76 (MQ=38)
   aTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAATACCCCATa                                                                  <  1:44064/76‑1 (MQ=255)
   aTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATa                                                                  <  2:1140772/76‑1 (MQ=255)
   aTTTTGATTCATCGTGATGACCGGTTTATCCCCGCCGGCGCGGGGGAAACCCCGAACGGGGGGGGAAAAAccccta                                                                   >  1:1073321/1‑74 (MQ=25)
     tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                <  2:59761/76‑1 (MQ=255)
       tGATTCATCGTGATGACCGGTTTATCCCCGCTGNCGCGGGGAACACCCTGAACGGGNGGGGAAAAGCNCATAAgcg                                                              >  2:1198468/1‑76 (MQ=37)
          ttCATCGTGATGACCGGTTTTTNNCCGCTGGCGCNNCGAACACCGTGAACAGTGGGGAAAAACCCCATAAGCGCTa                                                           <  1:1276020/76‑1 (MQ=255)
            cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                         >  2:154443/1‑76 (MQ=255)
             aTCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACt                                                        <  1:515510/76‑1 (MQ=255)
              tCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                       <  1:889149/76‑1 (MQ=255)
               cGTGATGACCGGTTTATCNCCGCTGGCGCGGGGAACACCTTGAACGCTGGGGAAAAACCCCATAAGCGCTAACTTa                                                      >  1:1366874/1‑76 (MQ=38)
               cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                      >  1:1397180/1‑76 (MQ=255)
               cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                      <  1:444109/76‑1 (MQ=255)
               cGTGATGACCGGTTTAGCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                      <  1:444104/76‑1 (MQ=255)
                 tgatgaCCGGTTTATCCCCCCTGGCGCGGGGGACACCCTGAACGGGGGGGGAAAAACCCCAAAACGCGAACTTAAg                                                    >  2:493983/1‑76 (MQ=21)
                 tgatgaCCGGTGTAGCCCCGCTGGCGCTGGGTACACCGTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAg                                                    <  1:1291290/76‑1 (MQ=37)
                   atgaCCGGNTTATCCCCGCTGGCGCGGGGAACACCCTGAGCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAggg                                                  <  1:1250607/76‑1 (MQ=255)
                        cGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg                                             <  2:920973/76‑1 (MQ=255)
                          gTTTACCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                           <  1:1033027/76‑1 (MQ=255)
                           tttATCCCCGCTGNCGCGGGGAACACCCTGCACGGTGGGGAAAAATCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                          <  2:1297591/76‑1 (MQ=37)
                            ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGCAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                         >  1:234944/1‑76 (MQ=255)
                             tacccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                        <  1:1035291/73‑1 (MQ=255)
                             tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAATCGTTAACTTAAGGGTTGTGGTAtt                                        >  1:1345362/1‑76 (MQ=255)
                               tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGGTGTGGTATTAc                                      >  2:627437/1‑76 (MQ=255)
                               cccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                      <  1:804798/75‑1 (MQ=255)
                                ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGNTGTGGTATTACg                                     <  1:21795/76‑1 (MQ=255)
                                    gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                 <  2:510118/76‑1 (MQ=255)
                                       ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                              <  2:639333/76‑1 (MQ=255)
                                        gcgcGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAGGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                             <  1:955026/76‑1 (MQ=255)
                                         cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                            <  1:777954/76‑1 (MQ=255)
                                          gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                           <  2:1234339/76‑1 (MQ=255)
                                             gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                        <  1:1266619/76‑1 (MQ=255)
                                             gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                        <  1:1266638/76‑1 (MQ=255)
                                               gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGGCTGATATGATTTaa                      <  2:798920/76‑1 (MQ=255)
                                               gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                      <  1:824895/76‑1 (MQ=255)
                                               gAACACCCNGAACGGTGGGGAAAAACCCCATAAGCGCTATCTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                      <  1:1247498/76‑1 (MQ=255)
                                                     ccTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                >  1:739906/1‑76 (MQ=255)
                                                     ccTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                >  1:925142/1‑76 (MQ=255)
                                                      cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg               >  1:780531/1‑76 (MQ=255)
                                                      cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg               <  1:63639/76‑1 (MQ=255)
                                                       tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa              >  2:629826/1‑76 (MQ=255)
                                                           cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa          <  1:790581/76‑1 (MQ=255)
                                                              gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt        <  2:460824/76‑1 (MQ=255)
                                                               tGGGGAAAACCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa       >  2:1226668/1‑76 (MQ=255)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa       <  1:149818/76‑1 (MQ=40)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa       >  1:1120581/1‑76 (MQ=40)
                                                                  ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc    >  2:1049445/1‑76 (MQ=35)
                                                                    aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAANGTGCCGATGAATTACtctc  <  1:1326982/76‑1 (MQ=17)
                                                                    aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc  >  1:376921/1‑76 (MQ=32)
                                                                          |                                                                     
GAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTC  >  REL606/2772781‑2772923

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