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 (0.830) 36/144 0.4 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)56 (1.080) 46/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)

AAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAC  >  REL606/2772782‑2772925
                                                                          |                                                                       
aatttttGTTTCATCGTGATGCCCGTTTTACCCCCGCGGGCGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCa                                                                      <  1:1145281/73‑1 (MQ=255)
 aaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                                                                     <  1:111500/76‑1 (MQ=255)
    tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGGGGGGAAAAACCCCATAAg                                                                  >  2:1609249/1‑76 (MQ=255)
     ttGATTCATCGTGATGACCGGTTTACCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgc                                                                 <  1:336138/76‑1 (MQ=255)
        tttCGTCGTGATGCCCGTTTTCCCCCCGGGGGCGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                                                              <  1:681408/75‑1 (MQ=255)
         ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGGGGGGGAAAAACCCCAAAGCGCTa                                                             >  2:1428863/1‑76 (MQ=25)
         ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                             <  2:1327004/76‑1 (MQ=255)
           cgtCGTGAGGACCGTTTTACCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                           <  1:467869/74‑1 (MQ=38)
           cATCGTGATGACCGGTTTCTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                           <  1:1215694/76‑1 (MQ=255)
           cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                           <  2:153597/76‑1 (MQ=255)
               gTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                       <  1:1543444/76‑1 (MQ=255)
               gTGATGACCGGTTTATCCCCGCCGGCGCGGGGAACACCCTGAACGGGGGGGAAAAACCCCATAAGCGCTAACTTaa                                                       >  2:454073/1‑76 (MQ=38)
                 gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                     <  1:1255794/76‑1 (MQ=255)
                 gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                     <  2:530939/76‑1 (MQ=255)
                    gaCCGGTTTATCCCCGCTGGGGCGGGGGACACCCCCCACCGGGGGGGCGAAACCCCCAAGCGCTAACTTAAggggt                                                  >  2:253263/1‑74 (MQ=9)
                    gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                  <  1:133332/76‑1 (MQ=255)
                     aCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAAAGGGGGGGAAAAACCCCCAAAGCGCCAACTTAAGGGGtg                                                 >  2:210928/1‑76 (MQ=21)
                      ccGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGGGGGGAAAAACCCCACAAGCGCTAACTTAAAGGgtgt                                                >  1:1458450/1‑76 (MQ=25)
                         gTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGGGGGGAAAAACCCCATAAGCGCTAACTTAAGGGGTGTGGt                                             >  2:1303154/1‑76 (MQ=255)
                          tttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                            >  2:1504091/1‑76 (MQ=255)
                           ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                           <  2:1198393/76‑1 (MQ=255)
                           ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                           >  2:312592/1‑76 (MQ=255)
                           ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                           >  2:790080/1‑76 (MQ=255)
                           ttATCCCCGCCGGCGGGGGGAACACCCTGAACGGTGGGGAAAAAACCCCTAAGCGCTAAATTAAGGATTGTGGGAt                                           >  2:1587736/1‑76 (MQ=255)
                            tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAGAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                          <  2:1595155/76‑1 (MQ=255)
                            tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                          <  1:1525083/76‑1 (MQ=255)
                            tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                          <  2:458307/76‑1 (MQ=255)
                              cccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                        <  1:735035/75‑1 (MQ=255)
                               ccccGCGGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                       <  2:1092273/76‑1 (MQ=255)
                                 ccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                     >  2:1248681/1‑76 (MQ=255)
                                  cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                    >  2:124633/1‑76 (MQ=255)
                                    cTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                  >  2:1481275/1‑76 (MQ=255)
                                         gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                             <  1:496664/76‑1 (MQ=255)
                                          cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                            >  1:1212023/1‑76 (MQ=255)
                                          cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                            >  1:626454/1‑76 (MQ=255)
                                           ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTACGGGATGTGGGATTCAGCCAGATATGAtt                           >  2:1138658/1‑76 (MQ=21)
                                            gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                          <  1:71382/76‑1 (MQ=255)
                                            gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                          <  2:594518/76‑1 (MQ=255)
                                            gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                          <  2:638367/76‑1 (MQ=255)
                                             ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTCTTACGCCTGATATGCTTTa                         >  1:196430/1‑76 (MQ=255)
                                             ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                         <  1:321002/76‑1 (MQ=255)
                                               aaCACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                       <  2:1656579/76‑1 (MQ=255)
                                                acacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                      >  2:1458711/1‑76 (MQ=255)
                                                 cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                     >  2:1244951/1‑76 (MQ=255)
                                                         aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg              <  1:578902/76‑1 (MQ=255)
                                                         aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg              <  1:484427/76‑1 (MQ=255)
                                                          aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa             <  2:1458450/76‑1 (MQ=255)
                                                           cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa            <  1:417475/76‑1 (MQ=255)
                                                           cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa            >  1:1379334/1‑76 (MQ=255)
                                                              gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt          <  2:160287/76‑1 (MQ=255)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa         <  1:1149002/76‑1 (MQ=40)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa         <  1:1014333/76‑1 (MQ=40)
                                                                   gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct     <  1:1212645/76‑1 (MQ=34)
                                                                    aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCNATGAATTACtctc    <  1:1182043/76‑1 (MQ=17)
                                                                    aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc    <  1:460323/76‑1 (MQ=32)
                                                                    aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc    <  1:709074/76‑1 (MQ=32)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa   <  1:1719078/76‑1 (MQ=31)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa   <  1:707538/76‑1 (MQ=31)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTAACGTGCCGATGAATTACTCTCAc  <  1:323201/76‑1 (MQ=9)
                                                                          |                                                                       
AAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAC  >  REL606/2772782‑2772925

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.