New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb 1209367 =16 (0.690)30 (1.810) 5/80 NT 76.5% coding (290/630 nt) ycfK hypothetical protein
?W3110S.gb 1211196 = 7 (0.420)pseudogene (37/537 nt) stfE ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment

GTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  W3110S.gb/1211135‑1211260
                                                             |                                                                
gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaa                                                         <  1:899484/71‑1 (MQ=255)
gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaa                                                         <  1:988213/71‑1 (MQ=255)
    cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     <  1:1782673/71‑1 (MQ=255)
    cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     <  1:1505160/71‑1 (MQ=255)
    cTGAAGGCCACGGAACCGGGACACCAACACGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     <  1:1733833/71‑1 (MQ=255)
        aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                  >  1:1065332/1‑70 (MQ=255)
              cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     <  1:878144/61‑1 (MQ=255)
                   tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      <  1:956523‑M2/29‑1 (MQ=255)
                   tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      <  1:1303124‑M2/29‑1 (MQ=255)
                   tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      <  1:1191276‑M2/29‑1 (MQ=255)
                   tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      <  1:1730971‑M2/29‑1 (MQ=255)
                   tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      <  1:736470‑M2/29‑1 (MQ=255)
                   tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      <  1:1007164‑M2/29‑1 (MQ=255)
                            ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                              <  1:1130449‑M2/37‑1 (MQ=255)
                            ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                              <  1:1073073‑M2/37‑1 (MQ=255)
                            ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                              <  1:336250‑M2/37‑1 (MQ=255)
                            ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                              <  1:515894‑M2/37‑1 (MQ=255)
                            ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                              <  1:1232248‑M2/37‑1 (MQ=255)
                            ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                              <  1:1428833‑M2/37‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:430773‑M2/40‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:933170‑M2/40‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:295934‑M2/40‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:32130‑M2/40‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:1236503‑M2/40‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:1240293‑M2/40‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:1702159‑M2/40‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:1275363‑M2/40‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:1350414‑M2/40‑1 (MQ=255)
                               cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:139034‑M2/40‑1 (MQ=255)
                                 ttgagttttgctccttctcccaaaccaaCGTTTATTAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:48006‑M2/29‑71 (MQ=255)
                                 ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                                 >  1:492012‑M2/29‑46 (MQ=255)
                                 ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                    >  1:464954‑M2/29‑43 (MQ=255)
                                 ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:1766889‑M2/29‑71 (MQ=255)
                                 ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:505400‑M2/29‑71 (MQ=255)
                                 ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:1492396‑M2/29‑71 (MQ=255)
                                   gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                      >  1:743192‑M2/27‑71 (MQ=255)
                                      ttttgctccttcacccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      <  1:522909‑M2/29‑1 (MQ=255)
                                        cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                  >  1:652290/1‑70 (MQ=255)
                                        cAGAGCCTTCTCCCAAACCAACGTATATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                  >  1:634643/1‑70 (MQ=255)
                                         tgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                 <  1:337228‑M2/50‑1 (MQ=255)
                                          gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc               <  1:711300/71‑1 (MQ=255)
                                          gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc               <  1:997683/71‑1 (MQ=255)
                                           ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc               <  1:245718‑M2/52‑1 (MQ=255)
                                           ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc               <  1:1817850‑M2/52‑1 (MQ=255)
                                           agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc               <  1:1700364/70‑1 (MQ=255)
                                              cttcacccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc               <  1:423594‑M2/52‑1 (MQ=255)
                                                       aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  <  1:1557065‑M2/65‑1 (MQ=255)
                                                       aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  <  1:707826‑M2/65‑1 (MQ=255)
                                                       aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  <  1:1635591‑M2/65‑1 (MQ=255)
                                                       aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  <  1:411842‑M2/65‑1 (MQ=255)
                                                       aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  <  1:369460‑M2/65‑1 (MQ=255)
                                                       aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  <  1:366145‑M2/65‑1 (MQ=255)
                                                       aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  <  1:205205‑M2/65‑1 (MQ=255)
                                                            aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                              <  1:539238‑M2/37‑1 (MQ=255)
                                                             |                                                                
GTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  W3110S.gb/1211135‑1211260

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