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)54 (0.900) 46/144 0.2 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)50 (0.840) 42/144 0.3 97.9% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 1 (0.020)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

CGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772780‑2772924
                                                                          |                                                                      
cGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCCCCGGCGCGGGGGACACCCCGAACGGGGGGGGAAAAacc                                                                       >  2:1069420/1‑76 (MQ=37)
    aTTTTTATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCAtn                                                                   <  1:355250/76‑2 (MQ=255)
        tGATTCATCGTGATGACCGGTTTATCCCCCCTGGCGCGGGGGACACCCTGAACGGGGGGGGAAAAACCCATAAgcg                                                               >  2:858797/1‑76 (MQ=255)
           ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACAGCCTGAACGGTGGGGAAAAACCTCATAAGCGATa                                                            <  1:1870390/76‑1 (MQ=37)
           ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                            <  1:1804710/76‑1 (MQ=255)
           ttCATCGTGATGACCGGTTTATCCCCCCCGGCGCGGGGGACAACCCGAACGGGGGGGGAAAACCCCCTAAGCGCTa                                                            >  1:1561562/1‑76 (MQ=255)
               tCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        <  1:728341/76‑1 (MQ=255)
               tCGTGATGACCGGTTTATCCCCGCGGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        <  1:1641185/76‑1 (MQ=255)
                cGTGATGACCGTTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                       <  1:523768/76‑1 (MQ=255)
                cGTGATGACCGGTTTATCCCCGTTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                       <  1:1485502/76‑1 (MQ=255)
                  tgatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAg                                                     <  1:1227638/76‑1 (MQ=255)
                   gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGAGGGGAAAAACCCCCTAAGCGCTAACTTAAgg                                                    >  2:1375284/1‑76 (MQ=255)
                    atgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTCAAggg                                                   <  1:1916603/76‑1 (MQ=255)
                      gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                 <  1:338856/76‑1 (MQ=255)
                      gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGGGGGGGAAAACCCCATAAGCGCTAACTTAGGGGtg                                                 >  2:217560/1‑75 (MQ=37)
                         cGGTTTATCCCCGCTGGCGCGGGGGACACCCTGAACGGGGGGGAAAAACCCCATAAGCGCTAACTTAAGGGGtgtg                                              >  2:1735405/1‑76 (MQ=38)
                          ggTTTATCTCCGCTGGCGCGCGGAACGATCTACACGGTGGGGCAAAACCCCATAAGCGCTGACGTAAGGGTTGTgg                                             <  1:1968825/76‑1 (MQ=16)
                          ggTTTATCCCCGCTGGCGCGGGGAACACCCTGATCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                             <  1:891783/76‑1 (MQ=255)
                          ggTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGTGAAAAACCCCATAAACGCTAACTTAAGGGTTGTgg                                             >  2:1623448/1‑76 (MQ=255)
                          ggTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGGGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                             >  2:407613/1‑76 (MQ=255)
                            tttCTCCCCGGCTGCGCGGGGAAACCCCTTAACGGTTGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                           <  2:1366461/76‑1 (MQ=18)
                            tttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                           <  2:1366451/76‑1 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          <  2:396095/76‑1 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          <  1:1368820/76‑1 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  2:1850408/1‑76 (MQ=255)
                              tATCCCCGCTGGCGCGGGGGACACCCTGAACGGGGGGGGAAAAACCCATAAGCGCTAACTTAAGGGGTGTGGGAtt                                         >  1:765961/1‑76 (MQ=25)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:373472/76‑1 (MQ=255)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         >  2:1912686/1‑76 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        >  2:1790518/1‑76 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  2:957154/76‑1 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        >  2:1837765/1‑76 (MQ=255)
                                tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTACGGGTTGTGGTATTAc                                       >  2:902831/1‑76 (MQ=255)
                                tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       >  2:902851/1‑76 (MQ=255)
                                  cccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                     <  1:1317477/76‑1 (MQ=255)
                                     gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAATGGTTGTGGTATTACGCCTg                                  <  1:1735339/76‑1 (MQ=255)
                                        ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                               >  2:659142/1‑76 (MQ=255)
                                         gcgcGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGntat                              <  1:1823863/76‑1 (MQ=255)
                                         gcgcGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                              <  1:1518687/76‑1 (MQ=255)
                                             ggggAACACCCTGAGCGATAGGCAAAAACCCCATAAGCGCTAACTTAAGGGTTGTCGTATTACGCCTGATATGAtt                          <  1:98161/76‑1 (MQ=17)
                                             ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          >  2:554655/1‑76 (MQ=255)
                                              gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                         <  2:1517776/76‑1 (MQ=255)
                                              gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                         <  2:1134789/76‑1 (MQ=255)
                                               ggAACACCCTGAACGGTGGGGAGAANCTCCATAAGCGCTAACTTAAGNGTTGTGGTATTACGCCTGATATGATTTa                        <  1:295988/76‑1 (MQ=18)
                                               ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        <  2:1868915/76‑1 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       >  1:1477964/1‑76 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       <  2:1032039/76‑1 (MQ=255)
                                                gAACACCCTGAACGGTGGCGANAANCCCCATACGCGCTAACTTAAGNGTTGTGGNANTACGCCTGATATGATTTaa                       >  2:293173/1‑76 (MQ=18)
                                                 aaCACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                      >  2:1548533/1‑76 (MQ=255)
                                                 aaCACCCTGAACGGTAGGGCAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATANGATTTAAc                      <  1:1964078/76‑1 (MQ=255)
                                                   cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    <  2:1038596/76‑1 (MQ=255)
                                                   cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    <  2:1038591/76‑1 (MQ=255)
                                                    acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  2:1412712/76‑1 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                >  1:1478373/1‑76 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACACCTGATATGATTTAACGTGCCg                >  1:1478386/1‑76 (MQ=255)
                                                         gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt              <  1:977092/76‑1 (MQ=255)
                                                           aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa            <  2:502071/76‑1 (MQ=255)
                                                            cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa           <  1:1213805/76‑1 (MQ=255)
                                                            cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATGTAACGTGCCGATGaa           <  1:178682/76‑1 (MQ=255)
                                                              gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         <  1:862741/76‑1 (MQ=255)
                                                              gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         <  1:1847073/76‑1 (MQ=255)
                                                              gTGGGGAAAAACCCCATAAGCGCTAACGTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         <  1:207730/76‑1 (MQ=255)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        <  1:1002526/76‑1 (MQ=40)
                                                               tGGGGAAAAACACCATCAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        >  1:119957/1‑76 (MQ=255)
                                                                   gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  2:112369/76‑1 (MQ=34)
                                                                   gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  1:367985/76‑1 (MQ=34)
                                                                    aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc   >  1:958390/1‑76 (MQ=32)
                                                                    aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc   <  1:330729/76‑1 (MQ=32)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  1:1533845/76‑1 (MQ=31)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  1:1034156/76‑1 (MQ=31)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  1:1034151/76‑1 (MQ=31)
                                                                          |                                                                      
CGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772780‑2772924

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