Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,994,992 Δ3,877 bp [ygeL]ygeQ [ygeL], ygeN, ygeO, insD‑4, insC‑4, ygeP, ygeQ

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2994992 2998868 3877 37 [0] [0] 39 [ygeL]–ygeQ [ygeL],ygeN,ygeO,insD‑4,insC‑4,ygeP,ygeQ

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 29949910 (0.000)36 (1.250) 31/94 0.0 100% pseudogene (102/156 nt) ygeL protein YgeL
?NC_000913 2998869 = 0 (0.000)intergenic (‑41/+115) ygeQ/glyU protein YgeQ/tRNA‑Gly

TAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATC  >  NC_000913/2994945‑2995039
                                              |                                                
tAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCac                                                <  1:1446821‑M1/49‑3 (MQ=255)
tAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCa                                                 <  1:1823506‑M1/48‑2 (MQ=255)
  gATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCactta                                             <  1:1709762‑M1/50‑6 (MQ=255)
  gATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCactt                                              <  1:1074742‑M1/49‑5 (MQ=255)
   aTTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttat                                            <  1:2083491‑M1/50‑7 (MQ=255)
   aTTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttat                                            <  1:206239‑M1/50‑7 (MQ=255)
       gTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttattta                                         <  1:2215967‑M1/49‑10 (MQ=255)
       gTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttattta                                         <  1:1834986‑M1/49‑10 (MQ=255)
        ttACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaac                                       <  1:95896‑M1/50‑12 (MQ=255)
        ttACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaa                                        <  1:482635‑M1/49‑11 (MQ=255)
          acacTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaa                                     >  1:1753102‑M1/1‑37 (MQ=255)
           cacTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaac                                    <  1:1897680‑M1/50‑15 (MQ=255)
            acTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaact                                   <  1:991488‑M1/50‑16 (MQ=255)
            acTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaact                                   <  1:841009‑M1/50‑16 (MQ=255)
             cTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaactt                                  <  1:2246596‑M1/50‑17 (MQ=255)
             cTCCGCTATAGCTTTCTTAAGATCTTTATAGCCCacttatttaacaactt                                  <  1:2068408‑M1/50‑17 (MQ=255)
                cGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttg                                 <  1:921620‑M1/48‑18 (MQ=255)
                 gCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaaa                              <  1:1021295‑M1/50‑21 (MQ=255)
                 gCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaaa                              <  1:1802631‑M1/50‑21 (MQ=255)
                 gCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaa                               >  1:1815157‑M1/1‑30 (MQ=255)
                  cTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaaaa                             <  1:1310000‑M1/50‑22 (MQ=255)
                     taGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaaaata                           <  1:1627753‑M1/49‑24 (MQ=255)
                        cTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaaaataaggt                       >  1:1170498‑M1/1‑23 (MQ=255)
                           tcttcAGATCTTTATAGCCCacttatttaacaacttgaaaataaggtgt                     <  1:1255250‑M1/49‑30 (MQ=255)
                              tcAGATCTTTATAGCCCacttatttaacaacttgaaaataaggtgtttta                 >  1:1799149‑M1/1‑17 (MQ=255)
                               cAGATCTTTATAGCCCacttatttaacaacttgaaaataaggtgtttta                 <  1:2032678‑M1/49‑34 (MQ=255)
                                aGATCTTTATAGCCCacttatttaacaacttgaaaataaggtgttttac                <  1:2154213‑M1/49‑35 (MQ=255)
                                 gATCTTTATAGCCCacttatttaacaacttgaaaataaggtgttttacct              <  1:519289‑M1/50‑37 (MQ=255)
                                   tCTTTATAGCCCacttatttaacaacttgaaaataaggtgttttacctgg            >  1:677900‑M1/1‑12 (MQ=255)
                                    cTTTATAGCCCacttatttaacaacttgaaaataaggtgttttacctggg           >  1:1012788‑M1/1‑11 (MQ=255)
                                      ttATAGCCCacttatttaacaacttgaaaataaggtgttttacctgggtt         >  1:90315‑M1/1‑9 (MQ=255)
                                       tataGCCCacttatttaacaacttgaaaataaggtgttttacctgggtt         <  1:1652773‑M1/49‑42 (MQ=255)
                                       tataGCCCacttatttaacaacttgaaaataaggtgttttacctg             <  1:50826‑M1/45‑38 (MQ=255)
                                        ataGCCCacttatttaacaacttgaaaataaggtgttttacctgggttgt       >  1:1776414‑M1/1‑7 (MQ=255)
                                           gCCCacttatttaacaacttgaaaataaggtgttttacctgggttgttac    <  1:1967258‑M1/50‑47 (MQ=255)
                                           gCCCacttatttaacaacttgaaaataaggtgttttacctgggttgtta     >  1:2002492‑M1/1‑4 (MQ=255)
                                              cacttatttaacaacttgaaaataaggtgttttacctgggttgttacaa  >  1:1662511‑M1/1‑1 (MQ=255)
                                              |                                                
TAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATC  >  NC_000913/2994945‑2995039

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