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 42 [0] [0] 43 [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)40 (1.520) 30/92 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

GTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTAT  >  NC_000913/2994944‑2995038
                                               |                                               
gTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCa                                                >  1:1764828‑M1/1‑48 (MQ=255)
 tAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCC                                                 <  1:1842776‑M1/47‑1 (MQ=255)
  aGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCact                                              >  1:43864‑M1/1‑46 (MQ=255)
   gATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCactt                                             <  1:127024‑M1/49‑5 (MQ=255)
   gATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCact                                              >  1:216679‑M1/1‑45 (MQ=255)
   gATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCact                                              >  1:375823‑M1/1‑45 (MQ=255)
   gATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCact                                              >  1:596537‑M1/1‑45 (MQ=255)
    aTTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCactta                                            >  1:1176876‑M1/1‑44 (MQ=255)
       tGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatt                                          <  1:1145333‑M1/48‑8 (MQ=255)
        gTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttattt                                         >  1:668610‑M1/1‑40 (MQ=255)
        gTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttattt                                         <  1:2035936‑M1/48‑9 (MQ=255)
         ttACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaa                                       <  1:355329‑M1/49‑11 (MQ=255)
         ttACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCacttattta                                        <  1:644975‑M1/48‑10 (MQ=255)
                 cGCTATAGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttga                               <  1:108892‑M1/49‑19 (MQ=255)
                    tataGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaaaa                            >  1:543940‑M1/1‑28 (MQ=255)
                    tataGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaaaa                            >  1:1283723‑M1/1‑28 (MQ=255)
                      taGCTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaaaata                          <  1:488975‑M1/49‑24 (MQ=255)
                         cTTTCTTCAGATCTTTATAGCCCacttatttaacaacttgaaaataagg                       <  1:120152‑M1/49‑27 (MQ=255)
                            tcttcAGATCTTTATAGCCCacttatttaacaacttgaaaataaggt                      <  1:290134‑M1/47‑28 (MQ=255)
                            tcttcAGATCTTTATAGCCCacttatttaacaacttgaaaataaggt                      <  1:264402‑M1/47‑28 (MQ=255)
                             cttcAGATCTTTATAGCCCacttatttaacaacttgaaaataaggtgt                    <  1:427386‑M1/48‑30 (MQ=255)
                              ttcAGATCTTTATAGCCCacttatttaacaacttgaaaataaggtgttt                  >  1:1148277‑M1/1‑18 (MQ=255)
                                cAGATCTTTATAGCCCacttatttaacaacttgaaaataaggtgtttt                 <  1:399796‑M1/48‑33 (MQ=255)
                                cAGATCTTTATAGCCCacttatttaacaacttgaaaataaggtgtttt                 <  1:1144365‑M1/48‑33 (MQ=255)
                                 aGATCTTTATAGCCCacttatttaacaacttgaaaataaggtgtttta                <  1:949032‑M1/48‑34 (MQ=255)
                                      tttATAGCCCacttatttaacaacttgaaaataaggtgttttacctggg          >  1:537102‑M1/1‑10 (MQ=255)
                                       ttATAGCCCacttatttaacaacttgaaaataaggtgttttacctgggt         >  1:561987‑M1/1‑9 (MQ=255)
                                        tataGCCCacttatttaacaacttgaaaataaggtgttttacctgggtt        >  1:1782695‑M1/1‑8 (MQ=255)
                                        tataGCCCacttatttaacaacttgaaaataaggtgttttacctgggtt        >  1:1556188‑M1/1‑8 (MQ=255)
                                        tataGCCCacttatttaacaacttgaaaataaggtgttttacctgggtt        >  1:2027770‑M1/1‑8 (MQ=255)
                                        tataGCCCacttatttaacaacttgaaaataaggtgttttacctggg          <  1:1371820‑M1/47‑40 (MQ=255)
                                        tataGCCCacttatttaacaacttgaaaataaggtgttttacctg            <  1:1865043‑M1/45‑38 (MQ=255)
                                        tataGCCCacttatttaacaacttgaaaataaggtgttttacct             <  1:1617391‑M1/44‑37 (MQ=255)
                                        tataGCCCacttatttaacaacttgaaaataaggtgtt                   <  1:1474357‑M1/38‑31 (MQ=255)
                                         ataGCCCacttatttaacaacttgaaaataaggtgttttacctgggtt        <  1:1857369‑M1/48‑42 (MQ=255)
                                         ataGCCCacttatttaacaacttgaaaataaggtgttttacctgggtt        <  1:23352‑M1/48‑42 (MQ=255)
                                         ataGCCCacttatttaacaacttgaaaataaggtgttttacctgggtt        <  1:1192571‑M1/48‑42 (MQ=255)
                                          taGCCCacttatttaacaacttgaaaataaggtgttttacctgggttgt      <  1:2014047‑M1/49‑44 (MQ=255)
                                            gCCCacttatttaacaacttgaaaataaggtgttttacctgggttgtta    >  1:477019‑M1/1‑4 (MQ=255)
                                            gCCCacttatttaacaacttgaaaataaggtgttttacctgggttgtta    >  1:890157‑M1/1‑4 (MQ=255)
                                            gCCCacttatttaacaacttgaaaataaggtgttttacctgggttgtt     <  1:1400657‑M1/48‑45 (MQ=255)
                                               cacttatttaacaacttgaaaataaggtgttttacctgggttgttaca  >  1:765039‑M1/1‑1 (MQ=255)
                                               |                                               
GTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTAT  >  NC_000913/2994944‑2995038

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: