New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 379832521 (1.050)8 (0.420) 4/126 1.6 26.7% coding (988/1074 nt) waaU lipopolysaccharide core biosynthesis
?NC_000913 = 3798516 24 (1.260)coding (797/1074 nt) waaU lipopolysaccharide core biosynthesis

TAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGA  >  NC_000913/3798260‑3798391
                                                                 |                                                                  
tAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTa                                                                  <  1:30115/68‑1 (MQ=255)
tAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTa                                                                  <  2:368002/68‑1 (MQ=255)
    ctctTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTg                                                              >  1:74698/1‑68 (MQ=255)
     tctTTTCTGAATTAGTTAAAGTTTCAGTATCGATAGCTTTTACCGTATATGTTGGGGATACTATTTGt                                                             <  1:677854/68‑1 (MQ=255)
     tctTTAGTGTATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTTTATGGTTGTGATACTAttttt                                                             <  1:677853/68‑3 (MQ=255)
          aCTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGTGATACTATTTGTATAGa                                                        <  2:16045/68‑1 (MQ=255)
                          tttCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTa                                                                  >  2:12386/1‑42 (MQ=255)
                          tttCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTa                                                                  <  1:12386/42‑1 (MQ=255)
                                  tCGATATCTTTTACCGTATATGTTGGGGATACcatttgatgagtttatttataccgttgcgttgacca                                <  2:222668‑M1/68‑37 (MQ=255)
                                   cGATATCTTTTACCGTATATGTTGGGGATACcatttgatgagtttatttataccgttgcgttgaccaa                               <  1:14336‑M1/68‑38 (MQ=255)
                                    gATATCTTTTACCGTATATGTTGGGGATACcatttgatgagtttatttataccgttgcgttgaacaag                              <  1:195890‑M1/68‑39 (MQ=255)
                                      tatCTTTTACCGTATATGTTGGGGATACcatttgatgagtttatttataccgttgcgttgaccaagta                            >  2:358834‑M1/1‑28 (MQ=255)
                                      tatCTTTTACCGTATATGTTGGGGATACcatttgatgagtttatttataccgttgcgttgaccaa                               >  1:183526‑M1/1‑28 (MQ=255)
                                      tatCTTTTACCGTATATGTTGGGGATACcatttgatgagtttatttataccgttgcgttgaccaa                               <  2:183526‑M1/65‑38 (MQ=255)
                                      tatCTTTTACCGTATATGTTGGGGATACaatttgatgagattaattataacgttgcgttgaccaagta                            >  2:195890‑M1/1‑28 (MQ=255)
                                      tatCTTTTACCGTATATGTTGGGGATACTATTTGTATAGa                                                        >  1:558712/1‑40 (MQ=255)
                                      tatCTTTTACCGTATATGTTGGGGATACTATTTGTATAGa                                                        <  2:558713/40‑1 (MQ=255)
                                      tatCTTATACCGTATATGTTGGGGATACcatttgatgagtttatttatacagttgagatgaacaagta                            >  2:14336‑M1/1‑28 (MQ=255)
                                              aCCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAggg                    <  2:407215/68‑1 (MQ=255)
                                              aCCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAggg                    <  2:111088/68‑1 (MQ=255)
                                                cGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGat                  >  1:494461/1‑68 (MQ=255)
                                                cGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGat                  >  1:280815/1‑68 (MQ=255)
                                                  tataTGTTGTGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGatat                <  1:18689/68‑1 (MQ=255)
                                                          ggggATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATTTTCAGGAGt        >  2:58219/1‑68 (MQ=255)
                                                          ggggATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGt        >  2:58218/1‑68 (MQ=255)
                                                            ggATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTAc      >  1:93934/1‑68 (MQ=255)
                                                            ggATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTAc      >  1:601186/1‑68 (MQ=255)
                                                             gATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACg     >  2:626234/1‑68 (MQ=255)
                                                                aCTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGa  <  2:280815/68‑1 (MQ=255)
                                                                 |                                                                  
TAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGA  >  NC_000913/3798260‑3798391

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