New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 379832528 (1.080)14 (0.570) 4/126 2.0 32.7% coding (988/1074 nt) waaU lipopolysaccharide core biosynthesis
?NC_000913 = 3798516 31 (1.260)coding (797/1074 nt) waaU lipopolysaccharide core biosynthesis

TGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAAT  >  NC_000913/3798449‑3798580
                                                                   |                                                                
tGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATg                                                                  <  2:555742/68‑1 (MQ=255)
tGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATg                                                                  <  2:277396/68‑1 (MQ=255)
  tATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGgt                                                                >  1:433799‑M2/1‑66 (MQ=255)
       aCAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAgtgt                                                           <  2:216311/68‑1 (MQ=255)
       aCAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAgtgt                                                           <  2:216310/68‑1 (MQ=255)
        cAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGtt                                                          >  1:91392/1‑68 (MQ=255)
            aaTGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGa                                                      <  1:441018/68‑1 (MQ=255)
            aaTGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGa                                                      <  1:456146/68‑1 (MQ=255)
                   aaaTCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAg                                               >  2:595985/1‑68 (MQ=255)
                    aaTCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAACGGCAGTGTTTCGACTTCAAgt                                              <  2:694334/68‑1 (MQ=255)
                            taCTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      >  2:626684‑M2/1‑40 (MQ=255)
                            taCTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAgtaatggtat                                      >  2:957088/1‑68 (MQ=255)
                             aCTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAgtattggtatt                                     >  2:597973/1‑68 (MQ=255)
                              cTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGgt                                                                <  2:315099‑M2/40‑3 (MQ=255)
                              cTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGgt                                                                >  1:315099‑M2/1‑38 (MQ=255)
                               ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      <  1:399146‑M2/65‑29 (MQ=255)
                               ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      >  2:399146‑M2/1‑37 (MQ=255)
                               ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      <  1:768181‑M2/65‑29 (MQ=255)
                               ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      <  2:656035‑M2/65‑29 (MQ=255)
                               ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      >  1:656035‑M2/1‑37 (MQ=255)
                               ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      >  2:768178‑M2/1‑37 (MQ=255)
                               ttGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGt                                   >  1:783304/1‑68 (MQ=255)
                                tGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagatatcga                                  <  1:183870‑M2/68‑33 (MQ=255)
                                tGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagatatcga                                  <  2:591610‑M2/68‑33 (MQ=255)
                                tGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      >  2:50562‑M2/1‑36 (MQ=255)
                                tGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      >  2:413970‑M2/1‑36 (MQ=255)
                                tGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      <  1:413970‑M2/64‑29 (MQ=255)
                                tGGTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagata                                      <  1:50562‑M2/64‑29 (MQ=255)
                                 ggTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCa                                 <  1:597973/68‑1 (MQ=255)
                                  gTCAACGCAACGGTATAAATAAACTCATCAAATGgtatccccaacatatacggtaaaagatatcgata                                >  2:183870‑M2/1‑34 (MQ=255)
                                       cgcaacgGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAAt                           >  2:638307/1‑68 (MQ=255)
                                            cgGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATATATCTTGc                      <  1:681319/68‑1 (MQ=255)
                                             ggTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCg                     >  2:389406/1‑68 (MQ=255)
                                                   aataaaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAAcc               >  2:415204/1‑68 (MQ=255)
                                                   aataaaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAAcc               <  2:469176/68‑1 (MQ=255)
                                                    ataaaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAAccc              >  2:920374/1‑68 (MQ=255)
                                                     taaaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCg             >  2:272073/1‑68 (MQ=255)
                                                     taaaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCg             >  2:862303/1‑68 (MQ=255)
                                                     taaaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCg             >  1:283566/1‑68 (MQ=255)
                                                      aaaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCgg            <  1:365773/68‑1 (MQ=255)
                                                      aaaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCgg            <  1:908645/68‑1 (MQ=255)
                                                      aaaCTCATCAAATGGCAGTGTTTCGATTCCAAGTATTGGTATTGTCAATAAATCTCGCGGTAACCCgg            <  2:274415/68‑1 (MQ=255)
                                                       aaCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGt           >  2:623331/1‑68 (MQ=255)
                                                        aCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCgg            <  1:177360/66‑1 (MQ=255)
                                                        aCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCgg            >  2:177360/1‑66 (MQ=255)
                                                        aCTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTa          >  1:44031/1‑68 (MQ=255)
                                                        aCTCATCAAATGGTAGTGTATCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTa          >  1:616854/1‑68 (MQ=255)
                                                           catcaAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAat       <  1:735900/68‑1 (MQ=255)
                                                           catcaAATGGTAGTGTTTCGATTTCAAGTATAGGAATTGTCAATAAATCTTGCGGTAACCCGGTAAat       <  1:501116/68‑1 (MQ=255)
                                                                aaTGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAAt  <  2:417612/68‑1 (MQ=255)
                                                                   |                                                                
TGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAAT  >  NC_000913/3798449‑3798580

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