Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 729,365 Δ9,598 bp rhsCybfD rhsC, ybfB, ybfO, ybfC, ybfQ, ybfL, ybfD

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 729365 738962 9598 38 [0] [0] 36 rhsC–ybfD rhsC,ybfB,ybfO,ybfC,ybfQ,ybfL,ybfD

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 7293640 (0.000)33 (1.150) 27/94 0.0 100% intergenic (+24/‑219) ybfA/rhsC DUF2517 domain‑containing protein YbfA/rhs element protein RhsC
?NC_000913 738963 = 0 (0.000)intergenic (+110/‑38) ybfD/ybgA H repeat‑associated putative transposase YbfD/DUF1722 domain‑containing protein YbgA

AAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAG  >  NC_000913/729315‑729413
                                                 |                                                 
aaaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAA                                                   <  1:489063‑M1/50‑1 (MQ=255)
aaaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAA                                                   <  1:1161151‑M1/50‑1 (MQ=255)
   aGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAACACATTGAAcgt                                                >  1:1859350‑M1/1‑47 (MQ=255)
     gCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAcgtct                                              <  1:152676‑M1/50‑6 (MQ=255)
       aaCTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAcgtctac                                            <  1:119853‑M1/50‑8 (MQ=255)
           ggTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAcgtctacact                                         >  1:2138943‑M1/1‑39 (MQ=255)
            gTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAcgtctacactt                                        <  1:1695975‑M1/49‑12 (MQ=255)
            gTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAcgtctacactt                                        <  1:1662019‑M1/49‑12 (MQ=255)
             tGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAcgtctacacttac                                      <  1:347054‑M1/50‑14 (MQ=255)
               aTTTTTATTGATTGCTGAGTTTTGAATACATTGAAcgtctacacttactc                                    >  1:1426788‑M1/1‑35 (MQ=255)
                 ttttATTGATTGCTGAGTTTTGAATACATTGAAcgtctacacttactcct                                  >  1:224596‑M1/1‑33 (MQ=255)
                  tttATTGATTGCTGAGTTTTGAATACATTGAAcgtctacacttactcct                                  <  1:1554798‑M1/49‑18 (MQ=255)
                   ttATTGATTGCTGAGTTTTGAATACATTGAAcgtctacacttactcctg                                 >  1:1799192‑M1/1‑31 (MQ=255)
                   ttATTGATTGCTGAGTTTTGAATACATTGAAcgtctacacttactcctg                                 >  1:962075‑M1/1‑31 (MQ=255)
                   ttATTGATTGCTGAGTTTTGAATACATTGAAcgtctacacttactcctg                                 >  1:1007535‑M1/1‑31 (MQ=255)
                    tATTGATTGCTGAGTTTTGAATACATTGAAcgtctacacttactcctgt                                <  1:1005308‑M1/49‑20 (MQ=255)
                    tATTGATTGCTGAGTTTTGAATACATTGAAcgtctacacttactcctgt                                <  1:733978‑M1/49‑20 (MQ=255)
                       tgattgCTGAGTTTTGAATACATTGAAcgtctacacttactcctgtaaac                            >  1:493143‑M1/1‑27 (MQ=255)
                        gattgCTGAGTTTTGAATACATTGAAcgtctacacttactcctgtaaacc                           >  1:223303‑M1/1‑26 (MQ=255)
                         attgCTGAGTTTTGAATACATTGAAcgtctacacttactcctgtaaaccg                          >  1:1337400‑M1/1‑25 (MQ=255)
                            gCTGAGTTTTGAATACATTGAAcgtctacacttactcctgtaaaccgct                        >  1:2170403‑M1/1‑22 (MQ=255)
                             cTGAGTTTTGAATACATTGAAcgtctacacttactcctgtaaaccgctca                      <  1:1413336‑M1/50‑30 (MQ=255)
                                aGTTTTGAATACATTGAAcgtctacacttactcctgtaaaccgctcag                     <  1:1814632‑M1/48‑31 (MQ=255)
                                     tGAATACATTGAAcgtctacacttactcctgtaaaccgctcaggagcag               >  1:466035‑M1/1‑13 (MQ=255)
                                      gAATACATTGAAcgtctacacttactcctgtaaaccgctcaggagcagt              >  1:47071‑M1/1‑12 (MQ=255)
                                      gAATACATTGAAcgtctacacttactcctgtaaaccgctcaggagcagt              >  1:1012958‑M1/1‑12 (MQ=255)
                                       aaTACATTGAAcgtctacacttactcctgtaaaccgctcaggagcagtaa            <  1:1287706‑M1/50‑40 (MQ=255)
                                         tACATTGAAcgtctacacttactcctgtaaaccgctcaggagcagtaatg          <  1:1608978‑M1/50‑42 (MQ=255)
                                         tACATTGAAcgtctacacttactcctgtaaaccgctcaggagcagtaat           >  1:1127380‑M1/1‑9 (MQ=255)
                                          aCATTGAAcgtctacacttactcctgtaaaccgctcaggagcagtaatg          >  1:1302320‑M1/1‑8 (MQ=255)
                                          aCATTGAAcgtctacacttactcctgtaaaccgctcaggagcagtaatg          >  1:922856‑M1/1‑8 (MQ=255)
                                           cATTGAAcgtctacacttactcctgtaaaccgctcaggagcagtaatga         <  1:1660898‑M1/49‑43 (MQ=255)
                                            aTTGAAcgtctacacttactcctgtaaaccgctcaggagcagtaatgaat       <  1:100192‑M1/50‑45 (MQ=255)
                                            aTTGAAcgtctacacttactcctgtaaaccgctcaggagcagtaatgaat       >  1:934184‑M1/1‑6 (MQ=255)
                                              tGAAcgtctacacttactcctgtaaaccgctcaggagcagtaatgaatc      >  1:770175‑M1/1‑4 (MQ=255)
                                                aacgtctacacttactcctgtaaaccgctcaggagcagtaatgaatctac   <  1:2137447‑M1/50‑49 (MQ=255)
                                                aacgtctacacttactcctgtaaaccgctcaggagcagtaatgaatcta    <  1:764644‑M1/49‑48 (MQ=255)
                                                 acgtctacacttactcctgtaaaccgctcaggagcagtaatgaatctaca  <  1:1719615‑M1/50‑50 (MQ=255)
                                                 |                                                 
AAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAG  >  NC_000913/729315‑729413

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: