Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 257,908 Δ776 bp insB9[crl] insB9, insA9, [crl]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 257908 258683 776 41 [0] [0] 40 insB9–[crl] insB9,insA9,[crl]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2579070 (0.000)32 (1.030) 23/82 0.1 100% intergenic (+8/+16) crl/insB9 RNA polymerase holoenzyme assembly factor Crl/IS1 transposase B
?NC_000913 258684 = 0 (0.000)pseudogene (9/331 nt) crl RNA polymerase holoenzyme assembly factor Crl

CAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCC  >  NC_000913/257858‑257955
                                                 |                                                
cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTA                                                  <  1:92549‑M1/50‑1 (MQ=255)
    ttGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtg                                              >  1:2401954‑M1/1‑46 (MQ=255)
    ttGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtg                                              >  1:483880‑M1/1‑46 (MQ=255)
     tGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtg                                              >  1:1297245‑M1/1‑45 (MQ=255)
       aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaa                                           <  1:2184935‑M1/50‑8 (MQ=255)
       aTCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaa                                           <  1:1293780‑M1/50‑8 (MQ=255)
        tCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaa                                          <  1:85398‑M1/50‑9 (MQ=255)
          aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaag                                         <  1:2358165‑M1/49‑10 (MQ=255)
          aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaag                                         <  1:2536364‑M1/49‑10 (MQ=255)
          aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagcgcaaag                                         <  1:1608778‑M1/49‑10 (MQ=255)
           aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagat                                       <  1:801624‑M1/50‑12 (MQ=255)
              aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataat                                    >  1:784929‑M1/1‑36 (MQ=255)
                tttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcg                                  <  1:295594‑M1/50‑17 (MQ=255)
                tttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcg                                  <  1:2036834‑M1/50‑17 (MQ=255)
                 ttACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcga                                 >  1:2507886‑M1/1‑33 (MQ=255)
                  tACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgat                                >  1:2286596‑M1/1‑32 (MQ=255)
                    ccGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattc                              <  1:543711‑M1/50‑21 (MQ=255)
                        aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                          <  1:2521263‑M1/50‑25 (MQ=255)
                        aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                          >  1:112174‑M1/1‑26 (MQ=255)
                        aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                          <  1:901942‑M1/50‑25 (MQ=255)
                        aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                          <  1:1670389‑M1/50‑25 (MQ=255)
                        aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                          >  1:1315253‑M1/1‑26 (MQ=255)
                           aGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                        >  1:704274‑M1/1‑23 (MQ=255)
                            ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcga                      <  1:684381‑M1/50‑29 (MQ=255)
                            ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                        <  1:325333‑M1/48‑27 (MQ=255)
                              cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                    <  1:1411738‑M1/50‑31 (MQ=255)
                               ccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattg                   >  1:783281‑M1/1‑19 (MQ=255)
                                cGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                  <  1:2670335‑M1/50‑33 (MQ=255)
                                 gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtc                 >  1:2115312‑M1/1‑17 (MQ=255)
                                 gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtc                 >  1:122930‑M1/1‑17 (MQ=255)
                                 gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtc                 <  1:126643‑M1/50‑34 (MQ=255)
                                    tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctg               <  1:1255014‑M1/49‑36 (MQ=255)
                                      ttCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggc             <  1:1731591‑M1/49‑38 (MQ=255)
                                          tGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat        <  1:627876‑M1/50‑43 (MQ=255)
                                           gAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatg       >  1:267715‑M1/1‑7 (MQ=255)
                                            aaGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc      >  1:540279‑M1/1‑6 (MQ=255)
                                            aaGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc      >  1:1269400‑M1/1‑6 (MQ=255)
                                             aGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcg     <  1:11486‑M1/50‑46 (MQ=255)
                                              ggTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgt    <  1:316888‑M1/50‑47 (MQ=255)
                                               gTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtc   <  1:578123‑M1/50‑48 (MQ=255)
                                                tAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtca  <  1:45225‑M1/50‑49 (MQ=255)
                                                 |                                                
CAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCC  >  NC_000913/257858‑257955

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