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)41 (1.130) 31/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)39 (1.070) 32/144 0.2 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

GAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTC  >  REL606/2772781‑2772923
                                                                          |                                                                      
gAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAcccc                                                                      <  2:926575/76‑1 (MQ=255)
  atttttAATTCATCGTAAGGCCCGTTTTCCCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                                                                    <  1:614330/74‑1 (MQ=25)
     tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                 <  1:725897/76‑1 (MQ=255)
         aTTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                                                             <  1:934593/76‑1 (MQ=255)
         aTTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                                                             <  1:556813/76‑1 (MQ=255)
           tCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGAGGGGAAAAAACCCCATAGCGCTaa                                                           >  2:1077719/1‑76 (MQ=25)
             aTCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGGACACCCCGAACGGGGGGGGAAAACCCCATAAGCGCTAACt                                                        >  2:679777/1‑76 (MQ=37)
              tCGTGATGACCGGTTTAGCCCCGCTGGCGCGGGGAACACCCTNAACGGTGGGGAAAAGCCCCGTAAGCGCTAACtt                                                       <  2:1032056/76‑1 (MQ=37)
               cGTGATGACCGGTTTACCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                      <  1:528128/76‑1 (MQ=255)
                gTGATGCCCGGTTTACCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                     <  1:597158/76‑1 (MQ=255)
                 tgatggCCNGTTTATCCCCGCTGGCGCGGGGAACACCCTCAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAg                                                    <  1:1066105/76‑1 (MQ=255)
                   atgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAggg                                                  >  2:94436/1‑76 (MQ=255)
                   atgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAACCCATACGCGCTAACTTAAggg                                                  >  2:38259/1‑76 (MQ=255)
                   atgaCCGGTTTATCCCCCCTGGCGCGGGGGAAACCCCGAACGGGGGGGGAAAAACCCCTCAGCGCTAACTTAAagg                                                  >  2:816919/1‑76 (MQ=255)
                     gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                <  2:366536/76‑1 (MQ=255)
                     gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                >  2:1055296/1‑76 (MQ=255)
                      aCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtg                                               <  2:546552/76‑1 (MQ=255)
                       ccGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgt                                              <  1:794740/76‑1 (MQ=255)
                       ccGGTTTATCCCCGCGGGCGCGTGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgt                                              <  2:41158/76‑1 (MQ=255)
                       ccGGTTTATCCCCCCCGGCGCGGGGGACACCCCGAACGGGGGGGGAAAAACCCCTAAGCGCCAACTTTAAGGGtgt                                               >  2:471753/1‑74 (MQ=16)
                         ggTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                            <  2:460460/76‑1 (MQ=255)
                           tttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAATGGTTGTGGTa                                          <  1:1111082/76‑1 (MQ=255)
                            ttACCCCCGCGGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                         <  1:575870/76‑1 (MQ=255)
                              aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATCAGAGCTAACTTAANGGTTGTGGTaatc                                       >  1:150989/1‑73 (MQ=38)
                              aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                       >  2:815020/1‑76 (MQ=255)
                              aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                       <  1:123940/76‑1 (MQ=255)
                               tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                      <  1:1055655/76‑1 (MQ=255)
                               cccccGCGGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                      <  1:426000/75‑1 (MQ=255)
                                 cccGCTGGCGCGGGGAACACCCTGAACGGTGGGGGAAAACCCCATAAGCGCTAACTTAAGGGGTGTGGGATTACGc                                    >  2:749415/1‑76 (MQ=25)
                                 cccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                    >  2:253688/1‑76 (MQ=255)
                                      tGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                               <  1:750134/76‑1 (MQ=255)
                                           cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                          >  1:624483/1‑76 (MQ=255)
                                           cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                          >  1:943045/1‑76 (MQ=255)
                                            ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACATAAGGGTTGTGGTATTACGCCTGATATGAtt                         >  2:196927/1‑76 (MQ=255)
                                                      cTGAACGGTGGGGAAAAACCCCATATGCGCTCACTTAAGCNNTGTGGTATTANNCCTGATATGATTTACCGTGCCg               >  1:1090732/1‑76 (MQ=17)
                                                       tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa              >  1:378011/1‑76 (MQ=255)
                                                        gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt             >  2:1170128/1‑76 (MQ=255)
                                                           cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa          <  1:848670/76‑1 (MQ=255)
                                                            ggTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt         <  1:658670/76‑1 (MQ=255)
                                                              tGGGGAAAAACCCCATAAGCGCTCACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa       <  1:1039030/76‑1 (MQ=255)
                                                              tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa       <  1:387236/76‑1 (MQ=40)
                                                                 ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc    <  2:782749/76‑1 (MQ=35)
                                                                 ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc    <  1:304231/76‑1 (MQ=35)
                                                                  gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct   >  2:623662/1‑76 (MQ=34)
                                                                    aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc  >  1:1164582/1‑76 (MQ=32)
                                                                          |                                                                      
GAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTC  >  REL606/2772781‑2772923

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.