Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 1,979,486 IS5 (+) +4 bp 100% intergenic (‑271/‑264) insA ← / → uspC IS1 repressor TnpA/universal stress protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1979486 =0 (0.000)75 (0.800) 55/178 NT 100% intergenic (‑271/‑267) insA/uspC IS1 repressor TnpA/universal stress protein
?NC_000913 2101749 = NA (NA)noncoding (1195/1195 nt) IS5 repeat region
* ? NC_000913 = 19794890 (0.000)71 (0.760) 49/178 NT 100% intergenic (‑274/‑264) insA/uspC IS1 repressor TnpA/universal stress protein
?NC_000913 = 2290113 NA (NA)noncoding (1/1195 nt) IS5 repeat region

TCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGA  >  NC_000913/2101658‑2101843
                                                                                           |                                                                                              
tCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGc                                                                                        >  2:1764965/1‑100 (MQ=255)
tCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGc                                                                                        >  2:1813171/1‑100 (MQ=255)
 cTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCg                                                                                       >  1:1157852/1‑100 (MQ=255)
 cTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCg                                                                                       <  2:1297718/100‑1 (MQ=255)
 cTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCg                                                                                       <  2:2081798/100‑1 (MQ=255)
  tGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGa                                                                                      >  2:659924/1‑100 (MQ=255)
   gAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGaa                                                                                     >  2:2078164/1‑100 (MQ=255)
     aaaaCGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATa                                                                                   >  1:28175/1‑100 (MQ=255)
     aaaaCGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATa                                                                                   <  1:2138948/100‑1 (MQ=255)
      aaaCGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATaa                                                                                  >  1:1298802/1‑100 (MQ=255)
        aCGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGc                                                                                >  1:1675221/1‑100 (MQ=255)
           tACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCggg                                                                             <  1:1813180/100‑1 (MQ=255)
           tACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCggg                                                                             <  2:140865/100‑1 (MQ=255)
            aCCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCgggg                                                                            <  2:206949/100‑1 (MQ=255)
              cAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGaa                                                                          >  2:1695227/1‑100 (MQ=255)
                gCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAt                                                                        >  2:1658548/1‑100 (MQ=255)
                 cATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAtt                                                                       >  1:928584/1‑100 (MQ=255)
                  aTAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACGTTCGTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAtt                                                                       >  2:1515249/1‑99 (MQ=255)
                   tAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAATCGGGGCAAttgt                                                                     >  2:2575101/1‑98 (MQ=255)
                   tAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAttct                                                                     >  1:1881224/1‑100 (MQ=255)
                   tAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAttct                                                                     <  1:2321547/100‑1 (MQ=255)
                   tAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAttct                                                                     <  1:2078175/100‑1 (MQ=255)
                   tAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAttct                                                                     <  2:1328843/100‑1 (MQ=255)
                     aaTGATCCTGCGTACCAATGAATACCCGTATAAGAATAGATACTTTCTTTAGGAATTTTTGTTTTTTTAAGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                   <  2:310953/100‑1 (MQ=37)
                     aaTGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                   <  2:2189714/100‑1 (MQ=255)
                      aTGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATtcttct                                                                  <  1:1680351/100‑1 (MQ=255)
                       tGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTtctc                                                                 >  2:2427674/1‑100 (MQ=255)
                          tCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGc                                                              <  2:2304558/100‑1 (MQ=255)
                           ccTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                                                             >  1:461170/1‑100 (MQ=255)
                            cTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAACTTCTTCTCGGCTg                                                            >  1:2203757/1‑100 (MQ=255)
                            cTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTg                                                            <  2:519479/100‑1 (MQ=255)
                             tGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGa                                                           >  2:730825/1‑100 (MQ=255)
                              gcgcACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGAc                                                          <  1:819581/100‑1 (MQ=255)
                               cgcACCAATCAACAACCGTATCAGAATAGATACATTCTGTAGGAATTTTTGTTGTATTAAGGCAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGTCt                                                         >  2:2260951/1‑100 (MQ=37)
                                   ccAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                     >  1:1450163/1‑100 (MQ=255)
                                     aaTCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtca                                                   >  1:1797998/1‑100 (MQ=255)
                                     aaTCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtca                                                   >  2:36089/1‑100 (MQ=255)
                                         aacaacCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttc                                               >  2:1374427/1‑100 (MQ=255)
                                           caacCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcat                                             >  1:930238/1‑100 (MQ=255)
                                             acCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTcatttcattt                                           >  2:246845/1‑100 (MQ=255)
                                                gTATCAGAATAGATACTTTCTTTAGGAATTTTTTTTTTTTTTTTGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt                                        >  2:1235497/1‑100 (MQ=18)
                                                 tATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGATGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcttc                                       <  2:1881995/100‑1 (MQ=255)
                                                     aGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAGTCTTCGCGGCTGACACAGTCATTTCATTTCTTCATGt                                   >  2:1757565/1‑100 (MQ=255)
                                                      gAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGtt                                  >  1:2451133/1‑100 (MQ=255)
                                                           gATACTTTCTTTAGGAATTTTTTTTTTTTTTTTTAAGGTGCGAATAAGCGGGGAAATTCTTTTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                              >  2:2333353/1‑99 (MQ=17)
                                                            aTACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGTTGCGAATAAGCGGGGAAATTCTTATCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGcc                            >  2:918858/1‑100 (MQ=255)
                                                             tACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCg                           >  2:1194356/1‑100 (MQ=255)
                                                               ctttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAt                         <  2:303384/100‑1 (MQ=255)
                                                               ctttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAt                         <  1:1311315/100‑1 (MQ=255)
                                                                tttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTAATTTCATTTCTTCATGTTTGAGCCGAtt                        >  2:407972/1‑100 (MQ=255)
                                                                   ctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAttttt                     >  2:932331/1‑100 (MQ=255)
                                                                    tttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                    <  2:2076706/100‑1 (MQ=255)
                                                                    tttAGGAATTTTTGATTTATGCACGAAGGTGCGAAGCATCGTTGAAATTCTTCACTACTGACTGATAAATTGCATGTCTGCAtgtgtg                                >  2:1141383/1‑88 (MQ=255)
                                                                     ttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtc                   <  1:1917393/100‑1 (MQ=255)
                                                                     ttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtc                   <  2:461166/100‑1 (MQ=255)
                                                                      tAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTATTCTCGGCTGACACAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                    >  2:1014788/1‑98 (MQ=255)
                                                                        ggAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTcc                >  2:1772621/1‑100 (MQ=255)
                                                                        ggAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTcc                <  1:1884891/100‑1 (MQ=255)
                                                                         gAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTccc               <  2:2090364/100‑1 (MQ=255)
                                                                          aaTTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTccc               <  2:591813/99‑1 (MQ=255)
                                                                               ttGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAt         <  1:1047395/100‑1 (MQ=40)
                                                                                   tttATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt     >  1:1209977/1‑100 (MQ=35)
                                                                                    ttattaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt    >  2:1533952/1‑100 (MQ=33)
                                                                                    ttattaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt    >  1:934694/1‑100 (MQ=33)
                                                                                     tattaAGGAAGGTGCGAATAAGCGGGTAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg   <  2:520508/100‑1 (MQ=17)
                                                                                     tattaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATTCCTTg   >  2:175798/1‑100 (MQ=17)
                                                                                      attaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGa  >  1:365450/1‑100 (MQ=31)
                                                                                           |                                                                                              
TCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGA  >  NC_000913/2101658‑2101843

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