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)48 (1.120) 39/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)46 (1.080) 39/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)

ATTTTGATTCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772784‑2772924
                                                                          |                                                                      
aTTTTGATTCATCGTGATGACCGGTTTAC‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATa                                                                   <  1:584413/76‑1 (MQ=255)
  tttGATTCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                 <  1:1201633/76‑1 (MQ=255)
     gATTCATCGTGATGACCGGTTTATCCCCCCCCGGGC‑‑GGGGGGGACCCCCCAACGGGGGGGGAAAAACCCCaaaacc                                                                >  1:358592/1‑74 (MQ=255)
         cATCGTGATGACCGGTTTAT‑‑CCCCCCCGGCGGGGGGGGGAAACCCAAACGGGGGGGGAAAAACCACAAAAgcggac                                                            >  1:1087616/1‑73 (MQ=255)
          aTCGTGATGACCGGTTTAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAca                                                         <  2:1417439/76‑2 (MQ=255)
          aTCGTGATGACCGGTTTAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGGGGGGGAAAAACCCATAAGCGCTAACt                                                         >  2:417941/1‑76 (MQ=37)
            cGTGATGACCGGTTTAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGGAAAACCCCATAACCGCTAACTTa                                                       >  2:1099358/1‑76 (MQ=255)
            cGTGATGACCGGTTTAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                       <  1:839089/76‑1 (MQ=255)
             gTGATGACCGGTTTAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGNGGAAAAACCCCATAAGCGCTAACTTaa                                                      <  1:1422539/76‑1 (MQ=255)
              tgatgaCCGGTTTAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAg                                                     <  2:1020921/76‑1 (MQ=255)
                       gTTTAC‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                            <  1:335523/76‑1 (MQ=255)
                        tttAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                           >  2:1315289/1‑76 (MQ=255)
                         ttAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          <  1:272924/76‑1 (MQ=255)
                         ttAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  1:540383/1‑76 (MQ=255)
                          tAT‑‑CCCCACTGGC‑‑GCGGGGAACACCCTGATCGGTGGGGAAAAACCCCATACGTGCTAACGTAAGGGTTGTGGTAtt                                         <  1:71556/76‑1 (MQ=37)
                           aT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  1:1007872/76‑1 (MQ=255)
                            t‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  1:863290/76‑1 (MQ=255)
                                  cGCTGGC‑‑GCGGGGAACACCCTGAACCGTNGGGAAAAACCACCTAAGCTCTAACTTAAGGGTTGTGGTATTACGCCt                                   <  1:119167/76‑1 (MQ=25)
                                   gCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  <  2:1051897/76‑1 (MQ=255)
                                    cTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                 <  1:1088417/76‑1 (MQ=255)
                                    cTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                 >  2:168088/1‑76 (MQ=255)
                                        c‑‑gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             <  1:1379900/76‑1 (MQ=255)
                                        c‑‑gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             <  1:307674/76‑1 (MQ=255)
                                             ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          <  2:860336/76‑1 (MQ=255)
                                             ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          >  2:291482/1‑76 (MQ=255)
                                              gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                         <  1:393462/76‑1 (MQ=255)
                                               ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        <  1:1192714/76‑1 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       <  1:640587/76‑1 (MQ=255)
                                                  acacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                     <  1:1169864/76‑1 (MQ=255)
                                                    acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  2:1109668/76‑1 (MQ=255)
                                                    acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  2:900707/76‑1 (MQ=255)
                                                    acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  2:1114188/76‑1 (MQ=255)
                                                          aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg             >  1:600315/1‑76 (MQ=255)
                                                           aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa            <  2:895466/76‑1 (MQ=255)
                                                            cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa           <  1:1144201/76‑1 (MQ=255)
                                                                 gggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct      <  2:1015296/76‑1 (MQ=38)
                                                                 gggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct      <  2:1015279/76‑1 (MQ=38)
                                                                  ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc     <  2:495907/76‑1 (MQ=35)
                                                                   gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  2:936764/76‑1 (MQ=34)
                                                                    aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc   >  2:712081/1‑76 (MQ=32)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGGATTACTCTCa  <  1:1319291/76‑1 (MQ=14)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGANTTACTCTCa  <  1:30008/76‑1 (MQ=14)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCNGATATGATTTAACGTGCCGATGAATTACTCTCa  <  2:1200137/76‑1 (MQ=14)
                                                                          |                                                                      
ATTTTGATTCATCGTGATGACCGGTTTAT‑‑CCCCGCTGGC‑‑GCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772784‑2772924

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