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)55 (1.180) 45/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)43 (0.920) 35/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)

AAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772782‑2772924
                                                                        |                                                                      
aaaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCa                                                                     <  2:1051418/76‑1 (MQ=255)
   ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                  <  1:112184/76‑1 (MQ=255)
    tttAATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                 <  2:1433385/76‑1 (MQ=255)
        aTTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                                                             <  1:353487/76‑1 (MQ=255)
             tCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        <  2:338123/76‑1 (MQ=255)
               gTGATCACCGATTTAACCCCGCTGGCGCGGGGAACACCCTGATCGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                      <  2:1181414/76‑1 (MQ=38)
                     aCCGGTTTATCCCCGCTGGCCCGGGGAACACCCAGAAAGGGGGGGAAAAACCCCCTAAGCGCTAACTTAAGGGGtg                                                >  2:71695/1‑76 (MQ=25)
                        ggTTTATCCCCGCTGGCGCGGGGTACACCCCGAACGGCGGGGGAAAACCCCCTAAGCGCTGAATTTAGGGTTGTgg                                             >  2:1116504/1‑76 (MQ=21)
                         gTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGCTGTggg                                            >  2:604773/1‑75 (MQ=255)
                           ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  2:1457316/1‑76 (MQ=255)
                            tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:1537993/76‑1 (MQ=255)
                            tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  2:1048881/76‑1 (MQ=255)
                            tATCCCCGCTGGCGCCGGGAACACCCTGAACGGCGGGGAAAAACCCCAAAACCGCCAACTTAACGGTTCTGGTAtt                                         >  1:1364152/1‑76 (MQ=25)
                             aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        >  2:1369294/1‑76 (MQ=255)
                             aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  1:61718/76‑1 (MQ=255)
                              tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  2:729954/76‑1 (MQ=255)
                              cccccGCGGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  1:899099/75‑1 (MQ=255)
                               ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      <  1:144035/76‑1 (MQ=255)
                                 ccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                    >  2:474946/1‑76 (MQ=255)
                                   tcTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  >  1:930941/2‑76 (MQ=255)
                                   gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  <  1:1178377/76‑1 (MQ=255)
                                   gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  <  2:828958/76‑1 (MQ=255)
                                   gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  <  2:144437/76‑1 (MQ=255)
                                         gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                            >  1:396464/1‑76 (MQ=255)
                                           ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          >  1:558466/1‑76 (MQ=255)
                                           ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          >  2:1017153/1‑76 (MQ=255)
                                           ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          <  1:1056826/76‑1 (MQ=255)
                                              gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       <  2:60397/76‑1 (MQ=255)
                                                acacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                     <  1:352191/76‑1 (MQ=255)
                                                 cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    <  1:874815/76‑1 (MQ=255)
                                                 cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    >  2:199416/1‑76 (MQ=255)
                                                  acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  1:174910/76‑1 (MQ=255)
                                                  acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  1:1389097/76‑1 (MQ=255)
                                                       gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGan              <  1:216802/76‑2 (MQ=255)
                                                       gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt              <  1:1387563/76‑1 (MQ=255)
                                                       gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt              <  1:162140/76‑1 (MQ=255)
                                                         aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGTTATGATTTAACGTGCCGATGa            <  2:470738/76‑1 (MQ=255)
                                                          cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa           <  1:1072602/76‑1 (MQ=255)
                                                          cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa           >  1:447613/1‑76 (MQ=255)
                                                           ggTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt          >  2:611501/1‑76 (MQ=255)
                                                             tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        >  2:105036/1‑76 (MQ=40)
                                                               gggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct      >  2:1339674/1‑76 (MQ=38)
                                                                ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc     <  2:534189/76‑1 (MQ=35)
                                                                 gtaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  2:611280/74‑1 (MQ=14)
                                                                 gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  2:708363/76‑1 (MQ=34)
                                                                  aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc   <  2:503548/76‑1 (MQ=32)
                                                                   aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  2:84852/76‑1 (MQ=31)
                                                                        |                                                                      
AAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772782‑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.