Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 3,894,997 Δ6,991 bp IS150‑mediated rbsD[yieO] rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 3893603–3894995 3901987 6993–8385 11 [10] [0] 47 [insJ‑5]–[yieO] [insJ‑5],insK‑5,rbsD,rbsA,rbsC,rbsB,rbsK,rbsR,[yieO]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 3894996NA (NA)46 (1.210) 34/150 0.2 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901988 = 0 (0.000)coding (432/1428 nt) yieO predicted multidrug or homocysteine efflux system

AAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGTCGAGCA  >  REL606/3901914‑3902063
                                                                          |                                                                           
ctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaAc                                                                            >  2:538239‑M2/75‑76 (MQ=255)
ctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaAc                                                                            >  1:946819‑M2/75‑76 (MQ=255)
 tgactccanttgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaAcc                                                                           >  1:1159371‑M2/74‑76 (MQ=255)
 tgactccanttgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaAcc                                                                           >  1:1165855‑M2/74‑76 (MQ=255)
 tgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaAcc                                                                           >  1:395585‑M2/74‑76 (MQ=255)
 tgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaAcc                                                                           <  1:639289‑M2/3‑1 (MQ=255)
 tgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaAcc                                                                           >  2:610617‑M2/74‑76 (MQ=255)
 tgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttgggatcagtacaAcc                                                                           <  1:64599‑M2/3‑1 (MQ=255)
  gactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaAccc                                                                          <  1:584328‑M2/4‑1 (MQ=255)
           tgaatatcggaaacatacctatatgcctcgtgtgtaactgtccaactttttggggtcagtanaACCCGGCATGGCg                                                                 >  1:1038781‑M2/64‑76 (MQ=255)
            gaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttgggggcagtacaACCCGGCAGGGCGa                                                                >  1:233473‑M2/63‑76 (MQ=255)
             aatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGAc                                                               <  2:1023375‑M2/15‑1 (MQ=255)
               tatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACaa                                                             >  1:678910‑M2/60‑76 (MQ=255)
                 tcggaatcagacctatatgcctcgtgtttaactgtccanctatttggggtcagtacaACCCGGCATGGCGACaaaa                                                           <  1:1185466‑M2/19‑1 (MQ=255)
                  cggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAAt                                                          <  1:1075227‑M2/20‑1 (MQ=255)
                   ggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAAtt                                                         >  2:588690‑M2/56‑76 (MQ=255)
                       tcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAAt                                                     <  1:409636‑M2/25‑1 (MQ=255)
                         agacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATAc                                                   >  2:487981‑M2/50‑76 (MQ=255)
                            cctatatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAAAACTgg                                                <  1:79382‑M2/30‑1 (MQ=255)
                             ctatatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGga                                               >  1:410715‑M2/46‑76 (MQ=255)
                              tatatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGgaa                                              >  1:705530‑M2/45‑76 (MQ=255)
                               atatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGgaag                                             <  2:85418‑M2/33‑1 (MQ=255)
                               atatgcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGgaag                                             >  1:521698‑M2/44‑76 (MQ=255)
                                  tgcctcgtgtttaactgtccaacttcttcgggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGaagaag                                          <  1:184599‑M2/36‑1 (MQ=255)
                                   gcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGt                                         <  1:968541‑M2/37‑1 (MQ=255)
                                   gcctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGt                                         <  1:968532‑M2/37‑1 (MQ=255)
                                     ctcgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTc                                       <  1:552928‑M2/39‑1 (MQ=255)
                                       cgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCAt                                     <  2:405060‑M2/41‑1 (MQ=255)
                                       cgtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCAt                                     <  2:428672‑M2/41‑1 (MQ=255)
                                        gtgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCAtt                                    <  2:65191‑M2/42‑1 (MQ=255)
                                         tgtttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTa                                   <  1:53029‑M2/43‑1 (MQ=255)
                                           tttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACg                                 >  2:423337‑M2/32‑76 (MQ=255)
                                           tttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACg                                 <  1:588690‑M2/45‑1 (MQ=255)
                                           tttaactgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACg                                 >  1:630462‑M2/32‑76 (MQ=255)
                                                ctgtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGAt                            <  2:655097‑M2/50‑1 (MQ=255)
                                                  gtccaactttttggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATaa                          >  1:805444‑M2/25‑76 (MQ=255)
                                                   tccaactttttggggtcagtacaACCCGGCATGGCGACAAAAATCAATACTGGAAGAAGTTCATAACGAGGATAAg                         <  1:1111508‑M2/53‑1 (MQ=255)
                                                         tttttgggntcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGCATAAGCGCGCa                   >  1:1162125‑M2/18‑76 (MQ=255)
                                                             tggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTaa               <  1:557808‑M2/63‑1 (MQ=255)
                                                              ggggtcagtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAg              <  1:473195‑M2/64‑1 (MQ=255)
                                                              ggggtcactacaACCCGGCATGGCGACAAAATTCAATACNNGAAGAAGNTCATTACGAGGATAAGCGCGCAGTAAg              <  1:1191010‑M2/64‑1 (MQ=255)
                                                                 gtcagtacaACNCGGCATGGCGACAAAATTCAATTCTGGTAGACNTTCATTACGAGGATAAGGGCGTAGTAAGGcc           <  1:144282‑M2/67‑1 (MQ=255)
                                                                    agtacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGt        >  1:431065‑M2/7‑76 (MQ=255)
                                                                      tacaACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGTCg      <  2:483605‑M2/72‑1 (MQ=255)
                                                                        caACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGTCGAg    >  2:142110‑M2/3‑76 (MQ=255)
                                                                         aACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGTCGAGc   >  1:162720‑M2/2‑76 (MQ=255)
                                                                          aCCCGGCATGGCGACAAAATTCAATACTGGAAGGAGTTCATTACGAGGATAAGCGCGCAGTNAGGCCAGTCGAGCa  >  1:1039426/1‑76 (MQ=255)
                                                                          |                                                                           
AAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGTCGAGCA  >  REL606/3901914‑3902063

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: