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)100 (1.910) 65/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)95 (1.820) 68/144 0.0 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

AAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772782‑2772924
                                                                          |                                                                      
aaaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCa                                                                     >  2:1579791/1‑76 (MQ=255)
  aTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATa                                                                   <  1:1141566/76‑1 (MQ=255)
  aTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATa                                                                   <  1:1230933/76‑1 (MQ=255)
        aTTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCt                                                             >  2:1697642/1‑76 (MQ=255)
         ttCATCGTGATGACCGGTTTATCCCCCCCGGCGCGGGGGGCAACCCCAACGGGGGGGGGAAAAACACCAAACgcgc                                                              >  1:1310826/1‑76 (MQ=255)
           cATCGTGATGACCGGTTTATCCCCCCCGGCGCGGGGGAAACCCTGTAC‑‑GGGGGGGGAAAAACCCATAAGCGCTAAc                                                          >  1:598720/1‑76 (MQ=25)
             tCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        <  2:129242/76‑1 (MQ=255)
              cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                       <  1:196146/76‑1 (MQ=255)
                tgatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGGAAAACCCCATAAGCGCTAACTTAAg                                                     >  1:1280824/1‑76 (MQ=255)
                tgatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAg                                                     <  2:882651/76‑1 (MQ=255)
                  atgaCCGGTTTATCCCCGCTGGCGCGGGGAAAACCCTGAAC‑‑GGTGGGGAAAAACCCCCTAAGGGGTAACTTAAggg                                                   >  2:628970/1‑76 (MQ=37)
                     aCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtg                                                <  1:611141/76‑1 (MQ=255)
                        ggTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                             >  2:1373930/1‑76 (MQ=255)
                         gTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                            <  1:1503704/76‑1 (MQ=255)
                          tttATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                           >  1:1066949/1‑76 (MQ=255)
                          tttATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                           >  2:1337303/1‑76 (MQ=255)
                          tttATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                           >  2:509241/1‑76 (MQ=255)
                           ttATCCCCGCTGGCGCGGGGGACACCCTGAAC‑‑GGGGGGGGAAAACCCCATAAGCGCTAACTTAAGGGGTGTGGTAt                                          >  2:610393/1‑76 (MQ=25)
                            tagccccGCGGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:1286422/73‑1 (MQ=255)
                            tagccccGATGGCTCGGGGAACACACTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  2:632017/73‑1 (MQ=25)
                            tATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATATGCGCTAACTTAAGGGTTGTGGTAtt                                         >  2:362147/1‑76 (MQ=255)
                            tATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         >  2:466211/1‑76 (MQ=255)
                            tATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:1286403/76‑1 (MQ=255)
                             aTCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  1:290250/76‑1 (MQ=255)
                             aTCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  2:701189/76‑1 (MQ=255)
                             aTCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  1:1587504/76‑1 (MQ=255)
                             aTCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGCGGGGAAAAACCCCATAAGCGGTAACTTAAGGGTTGTGGTATTa                                        >  2:1538281/1‑76 (MQ=255)
                                     tGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                >  1:1423487/1‑76 (MQ=255)
                                      ggCGCGGGGTACACCCTGAAC‑‑GGTGGGGAAAAACCCGGTAAGTGCTAACTTAAGGGATGTGGCATGACGCCTGatg                               >  2:413363/1‑75 (MQ=255)
                                        cgcgGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             <  1:1566458/76‑1 (MQ=255)
                                        cgcgGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             >  1:695398/1‑76 (MQ=255)
                                          cgGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                           >  1:497965/1‑76 (MQ=255)
                                                acacCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                     >  2:1315544/1‑76 (MQ=255)
                                                acacCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                     >  2:239254/1‑76 (MQ=255)
                                                 cacCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    >  2:1411133/1‑76 (MQ=255)
                                                 cacCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    <  2:1683046/76‑1 (MQ=255)
                                                  acCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  1:824845/76‑1 (MQ=255)
                                                    ccTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                 >  2:1525218/1‑76 (MQ=255)
                                                     cTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                >  1:210646/1‑76 (MQ=255)
                                                      tGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa               <  2:517892/76‑1 (MQ=255)
                                                        aaC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg             >  2:1097727/1‑76 (MQ=255)
                                                         aC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa            <  1:1172954/76‑1 (MQ=255)
                                                              gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         <  2:451318/76‑1 (MQ=255)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        >  2:1115998/1‑76 (MQ=40)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        <  1:31551/76‑1 (MQ=40)
                                                                ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc       <  1:247530/76‑1 (MQ=39)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  2:1274296/76‑1 (MQ=31)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  1:672296/76‑1 (MQ=31)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  >  1:1368627/1‑76 (MQ=31)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  >  1:109451/1‑76 (MQ=31)
                                                                          |                                                                      
AAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAC‑‑GGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772782‑2772924

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.