Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,635,363 IS5 (–) +4 bp intergenic (‑53/‑155) rob ← / → creA DNA‑binding transcriptional dual regulator Rob/PF05981 family protein CreA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2288919 =NA (NA)55 (1.320) 27/128 0.1 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 4635366 0 (0.000)intergenic (‑56/‑155) rob/creA DNA‑binding transcriptional dual regulator Rob/PF05981 family protein CreA
* ? NC_000913 = 3653230NA (NA)70 (1.680) 33/128 0.0 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 4635363 = 0 (0.000)intergenic (‑53/‑158) rob/creA DNA‑binding transcriptional dual regulator Rob/PF05981 family protein CreA

TCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGG  >  NC_000913/3653175‑3653289
                                                       |                                                           
tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGac                                                    >  1:961093/1‑65 (MQ=35)
tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGac                                                    >  1:1063006/1‑65 (MQ=35)
tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGac                                                    >  1:3033198/1‑65 (MQ=35)
tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGac                                                    >  1:2882156/1‑65 (MQ=35)
tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGac                                                    >  1:2612685/1‑65 (MQ=35)
tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGac                                                    >  1:1449868/1‑65 (MQ=35)
 cTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGaca                                                   <  1:580662/65‑1 (MQ=35)
   aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACt                                                 >  1:2780921/1‑65 (MQ=39)
   aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACt                                                 >  1:1882913/1‑65 (MQ=39)
     ggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAt                                               >  1:2496888/1‑65 (MQ=255)
     ggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAt                                               >  1:1354223/1‑65 (MQ=255)
     ggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGGGACACTAt                                               >  1:2725455/1‑65 (MQ=255)
           aGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAttttttt                                         >  1:169524/1‑65 (MQ=255)
             aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTAc                                       >  1:1829828/1‑65 (MQ=255)
              tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACa                                      >  1:2010673/1‑65 (MQ=255)
                 cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATaa                                   <  1:3140580/65‑1 (MQ=255)
                  aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATaaa                                  <  1:1684660/65‑1 (MQ=255)
                  aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATaaa                                  <  1:2835788/65‑1 (MQ=255)
                  aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATaaa                                  <  1:21876/65‑1 (MQ=255)
                  aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATaaa                                  >  1:2031089/1‑65 (MQ=255)
                    aCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGa                                <  1:2487310/65‑1 (MQ=255)
                      aTGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAac                              <  1:2321939/65‑1 (MQ=255)
                      aTGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATCCTTTTACATAAAGAac                              >  1:1993489/1‑65 (MQ=255)
                       tGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAacg                             <  1:285301/65‑2 (MQ=255)
                        gATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAacac                            <  1:1069183/65‑1 (MQ=255)
                         aTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACAcc                           <  1:1539217/65‑1 (MQ=255)
                         aTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACAcc                           <  1:866191/65‑1 (MQ=255)
                         aTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACAcc                           <  1:2513160/65‑1 (MQ=255)
                           ctcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCt                         >  1:772033/1‑65 (MQ=255)
                           ctcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCt                         >  1:293170/1‑65 (MQ=255)
                            tcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTa                        >  1:1269655/1‑65 (MQ=255)
                            tcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTa                        >  1:101424/1‑65 (MQ=255)
                              atatCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAcc                      <  1:639025/65‑1 (MQ=255)
                                atCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAcccc                    >  1:1725652/1‑65 (MQ=255)
                                atCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAcccc                    >  1:64656/1‑65 (MQ=255)
                                 tCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccccc                   >  1:1382397/1‑65 (MQ=255)
                                 tCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccccc                   >  1:3154903/1‑65 (MQ=255)
                                 tCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccccc                   >  1:2925422/1‑65 (MQ=255)
                                    gggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                >  1:1577031/1‑65 (MQ=255)
                                       aCTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATta             >  1:3119141/1‑65 (MQ=255)
                                        cTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATtat            >  1:1180093/1‑65 (MQ=255)
                                        cTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATtat            <  1:748972/65‑1 (MQ=255)
                                        cTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATtat            >  1:1036592/1‑65 (MQ=255)
                                            ttCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATaa        >  1:642520/1‑65 (MQ=255)
                                             tCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATaaa       >  1:1704710/1‑65 (MQ=255)
                                                cACCTTCCTTAGGTGACACTGTTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAt    <  1:3136203/65‑1 (MQ=255)
                                                 aCCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATg   >  1:1120606/1‑65 (MQ=255)
                                                 aCCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATg   >  1:1978808/1‑65 (MQ=255)
                                                 aCCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATg   >  1:934965/1‑65 (MQ=255)
                                                 aCCTTCCTTAGGCGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATg   >  1:3225506/1‑65 (MQ=255)
                                                 aCCTTCCTTAGGCGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATg   >  1:2821766/1‑65 (MQ=255)
                                                  ccttccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATgg  <  1:2772043/65‑1 (MQ=255)
                                                  ccttccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATgg  <  1:2000427/65‑1 (MQ=255)
                                                       |                                                           
TCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGG  >  NC_000913/3653175‑3653289

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