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.100) 35/144 0.1 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)34 (0.870) 29/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)

AAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAC  >  REL606/2772782‑2772919
                                                                          |                                                                 
aaaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCa                                                                <  1:688302/76‑1 (MQ=255)
    tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                            >  1:878480/1‑76 (MQ=255)
      tGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgcg                                                          >  1:995746/1‑76 (MQ=255)
          tCATCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                                                      <  1:755780/76‑1 (MQ=255)
          tCATCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                                                      <  1:1125295/76‑1 (MQ=255)
          tCATCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGGGGGGGAAAACCCCATAAGCGCTaa                                                      >  1:361363/1‑76 (MQ=255)
           cagcGTGATGCCGGGTTTACCCCCGCGGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                     <  1:594379/73‑1 (MQ=37)
            aTCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGGACACCCTGAACGGTGGGGAAAAACCCCCTAAGCGCTAACt                                                    >  2:733013/1‑76 (MQ=255)
            aTCGNGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACANCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACt                                                    <  1:149986/76‑1 (MQ=255)
             tCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                   <  2:450120/76‑1 (MQ=255)
              cGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAANCCCCATAAGCGCTAACTTa                                                  <  1:1219375/76‑1 (MQ=255)
               gTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGGGGGGGAAAACCCAATAAGCGCTAACTTaa                                                 >  2:322684/1‑76 (MQ=37)
                 gatgaCCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAAACCCATAAGCGCTAACTTTAgg                                               >  2:460830/1‑76 (MQ=255)
                   tgaCCGGTTTATCCCCCCCGGCGCGGGGGGAAAACCCAAACGGGGGGGGAAAACCCCCAACAGCGCACCTTAAggg                                              >  2:364295/1‑76 (MQ=255)
                    gggcGGTTTTCCCCCGGGGGGGG‑‑GGGGGACCCCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                            <  1:433555/73‑1 (MQ=18)
                      ccGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGCTtgt                                          <  2:1280290/76‑1 (MQ=255)
                         gTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                       <  1:1200000/76‑1 (MQ=255)
                          tttATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACNGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                      >  2:1106384/1‑76 (MQ=255)
                           ttATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGACAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                     >  2:836482/1‑76 (MQ=255)
                           ttATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                     <  2:996411/76‑1 (MQ=255)
                           ttACCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                     <  1:429561/76‑1 (MQ=255)
                            tATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGGTGTGGTAtt                                    >  2:455716/1‑76 (MQ=255)
                             aTCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                   >  2:652278/1‑76 (MQ=255)
                                   gcgggCGC‑‑GGGGANCACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                             <  1:174394/73‑1 (MQ=255)
                                   gCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                             <  1:826075/76‑1 (MQ=255)
                                   gCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                             <  2:743670/76‑1 (MQ=255)
                                    cTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                            <  1:685658/76‑1 (MQ=255)
                                     tGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                           <  1:742478/76‑1 (MQ=255)
                                     tGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTCAAGGGTTGTGNTATTACGCCTTat                           >  1:167803/1‑76 (MQ=255)
                                      ggCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                          <  1:428847/76‑1 (MQ=255)
                                          c‑‑gGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                      >  1:1017779/1‑76 (MQ=255)
                                          c‑‑gGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                      <  2:758115/76‑1 (MQ=255)
                                             ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                     >  1:1106215/1‑76 (MQ=255)
                                              gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                    >  1:181112/1‑76 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAGACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                  <  1:1115124/76‑1 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                  <  1:400817/76‑1 (MQ=255)
                                                 aaCACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                 <  1:95965/76‑1 (MQ=255)
                                                         gAACGGTGGGGAAANACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGan         <  2:178195/76‑2 (MQ=255)
                                                          aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg        <  2:361363/76‑1 (MQ=255)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTNAGGGTTGTGATATTACGCCTGATATGATTTAACGTGCCGATGAATTa   <  1:111460/76‑1 (MQ=255)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa   <  1:961281/76‑1 (MQ=40)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa   <  1:1100952/76‑1 (MQ=40)
                                                                ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc  >  1:1069710/1‑76 (MQ=39)
                                                                          |                                                                 
AAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGC‑‑GGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAC  >  REL606/2772782‑2772919

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