New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 1145450100 (0.970)8 (0.100) 4/36 NT 10.0% intergenic (‑15/‑140) flgA/flgB flagellar basal body P‑ring biosynthesis protein A/flagellar basal‑body rod protein B
?REL606 = 3613817 69 (0.900)coding (1255/1950 nt) yhjK predicted diguanylate cyclase

Only 100 of 159 total aligned reads displayed.

GCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCC  >  REL606/3613783‑3613851
                                  |                                  
gCACTGACCAGTTTACCGCTGGTCATCTAGACCtg                                    <  1:6162882/35‑1 (MQ=255)
 cACTGACCAGTTTACCGCTGGTCATCTCGACCTGt                                   <  1:17822454‑M2/35‑2 (MQ=255)
 cACTGACCAGTTTACCGCTGGTCATCTCGACCTGt                                   <  1:2730403‑M2/35‑2 (MQ=255)
  aCTGACCAGTTTACCGCTGGTCATCTCGACCTGtg                                  <  1:16560975‑M2/35‑3 (MQ=255)
  aCTGACCAGTTTACCGCTGGTCATCTCGACCTGtg                                  <  1:5629024‑M2/35‑3 (MQ=255)
  aCTGACCAGTTTACCGCTGGTCATCTCGACCTGtg                                  <  1:13664452‑M2/35‑3 (MQ=255)
   cTGACCAGTTTACCGCTGGTCATCTCGACCTGtgg                                 <  1:5996164‑M2/35‑4 (MQ=255)
   cTGACCAGTTTACCGCTGGTCATCTCGACCTGtgg                                 <  1:15419607‑M2/35‑4 (MQ=255)
   cTGACCAGTTTACCGCTGGTCATCTCGACCTGtgg                                 <  1:4176129‑M2/35‑4 (MQ=255)
    tGACCAGTTTACCGCTGGTCATCTCGACCTGtggc                                <  1:18945157‑M2/35‑5 (MQ=255)
    tGACCAGTTTACCGCTGGTCATCTCGACCTGtggc                                <  1:20185828‑M2/35‑5 (MQ=255)
    tGACCAGTTTACCGCTGGTCATCTCGACCTGtggc                                <  1:15296884‑M2/35‑5 (MQ=255)
     gcccATTTTCCCGCTGGTCATCTCGACCTGTGGCt                               <  1:16071771/33‑1 (MQ=25)
     gACCAGTTTCCCGCTGGTCATCTCGACCTGtggct                               >  1:385226‑M2/1‑30 (MQ=255)
     gACCAGTTTCCCGCTGGTCATCTCGACCTGtggct                               >  1:11420197‑M2/1‑30 (MQ=255)
     gACCAGTTTACCGCTGGTCATCTCGACCTGtggct                               <  1:3723687‑M2/35‑6 (MQ=255)
     gACCAGTTTACAGCTGGTCATCTCGACCTGtggct                               >  1:7027816‑M2/1‑30 (MQ=255)
      aCCATTTTACCGCTGGTCATCTCGACCTGtggctg                              >  1:6897450‑M2/1‑29 (MQ=255)
      aCCATTTTACCGCTGGTCATCTCGACCTGtggctg                              >  1:13628533‑M2/1‑29 (MQ=255)
       ccAGTTTCCCGCTGGTCATCTCGACCTGTGGCTGt                             <  1:16364023/35‑1 (MQ=255)
       ccAGTTTACCGCTGGTCATCTCGACCTGtggctgg                             <  1:9601908‑M2/35‑8 (MQ=255)
        cAGTTTCCCGCTGGTCATCTCGACCTGTGGCTGTa                            <  1:10568966/35‑1 (MQ=255)
        cAGTTTCCCGCTGGTCATCTCGACCTGTGGCTGTa                            <  1:17356736/35‑1 (MQ=255)
        cAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTa                            <  1:3437096/35‑1 (MQ=255)
         aGTTTACCGCTGGTCATCTCGACCTGtggctgggg                           <  1:8356968‑M2/35‑10 (MQ=255)
         aGTTTACCGCTGGTCATCTCGACCTGTGGCTGTaa                           <  1:14642500/35‑1 (MQ=255)
          ttttCCCGCTGGTCATCTCGACCTGTGGCTGTAAc                          <  1:11748022/34‑1 (MQ=255)
          gTTTCCCGCTGGTCATCTCGACCTGTGGCTGTAAc                          <  1:17092226/35‑1 (MQ=255)
          gTTTCCCGCTGGTCATCTCGACCTGTGGCTGTAAc                          <  1:5806044/35‑1 (MQ=255)
          gTTTACCGCTGGTCATCTCGACCTGtggctggggg                          <  1:8867781‑M2/35‑11 (MQ=255)
          gTTTACCGCTGGTCATCTCGACCTGtggctggggg                          <  1:8329763‑M2/35‑11 (MQ=255)
          gTTTACCGCTGGTCATCTCGACCTGtggctggggg                          <  1:4234132‑M2/35‑11 (MQ=255)
           tttCCCGCTGGCCATCTCGACCTGTGGCTGTAAcc                         <  1:5324973/35‑1 (MQ=37)
           tttACCGCTGGTCATCTCGACCTGtggctgggggt                         <  1:10148076‑M2/35‑12 (MQ=255)
           tttACCGCTGGTCATCTCGACCTGtggctgggggg                         <  1:1139372‑M2/35‑12 (MQ=255)
           tttACCGCTGGTCATCTCGACCTGtggctgggggg                         <  1:9341412‑M2/35‑12 (MQ=255)
           tttACCGCTGGTCATCTCGACCTGtggctgggggg                         <  1:6472871‑M2/35‑12 (MQ=255)
           tttACCGCTGGTCATCTCGACCTGTGGCTGTAAcc                         <  1:19759909/35‑1 (MQ=255)
            ttcccGCTGGTCATCTCGACCTGTGGCTGTAACCa                        <  1:15283915/32‑1 (MQ=255)
            ttcccGCTGGTCATCTCGACCTGTGGCTGTAACCa                        <  1:10573693/32‑1 (MQ=255)
            ttACCGCTGGTCATCTCGACCTGtggctggggggg                        <  1:1663014‑M2/35‑13 (MQ=255)
            ttACCGCTGGTCATCTCGACCTGtggctgggagcg                        <  1:5267305‑M2/35‑13 (MQ=255)
            ttACCGCTGGTCATCTCGACCTGTGGCTGTAACCa                        <  1:3502054/35‑1 (MQ=255)
             tACCGCTGGTCATCTCGACCTGtggctggggggga                       <  1:11139360‑M2/35‑14 (MQ=255)
              cccGCTGGTCATCTCGACCTGTGGCTGTAACCaaa                      <  1:13275193/34‑1 (MQ=255)
              cccGCTGGTCATCTCGACCTGTGGCTGTAACCaaa                      <  1:17002802/34‑1 (MQ=255)
              aCCGCTGGTCATCTCGACCTGtggctgggggcgaa                      <  1:3871715‑M2/35‑15 (MQ=255)
              aCCGCTGGTCATCTCGACCTGTTGCTGTAACTaaa                      >  1:569179/1‑35 (MQ=37)
               ccGCTGGTCATCTCGCCCTGTGGCTGTAACCAAAt                     <  1:12369012/35‑1 (MQ=255)
               ccGCTGGTCATCTCGACCTGTGGCTGTAACCAAAt                     <  1:11786727/35‑1 (MQ=255)
               ccGCTGGTCATCTCGACCTGTGGCTGTAACCAAAt                     <  1:12810923/35‑1 (MQ=255)
               ccGCTGGTCATCTCGACCTGTGGCTGTAACCAAAt                     <  1:9557591/35‑1 (MQ=255)
               ccGCTGGTCATCTCGACCTGTGGCTGTAACCAAAt                     <  1:9337051/35‑1 (MQ=255)
               ccGCTGGTCATCTCGACCTGTGGCTGTAACCAAAt                     <  1:14525974/35‑1 (MQ=255)
               ccGCAGGTCATCTCGACCTGTGGCTGTAACCAAAt                     <  1:6101809/35‑1 (MQ=255)
                cGCTGGTCATCTCGACCTGTGGCTGTAACCAAATa                    <  1:2346857/35‑1 (MQ=255)
                cGCTGGTCATCTCGACCTGTGGCTGTAACCAAATa                    <  1:228184/35‑1 (MQ=255)
                 gCTGGTCATCTCGACCTGtggctgggggggaaaaa                   <  1:18138525‑M2/35‑18 (MQ=255)
                  cTGGTCATCTCGACCTGtggctgggggcgaaaagc                  <  1:2976972‑M2/35‑19 (MQ=255)
                  cTGGTCATCTCGACCTGTGGCTGTAACCAAATAGc                  <  1:10119005/35‑1 (MQ=255)
                  cTGGTCATCTCGACCTGTGGCTGTAACCAAATAGc                  <  1:13526058/35‑1 (MQ=255)
                  cTGGTCATCTCGACCTGTGGCTGTAACCAAATAGc                  <  1:17755312/35‑1 (MQ=255)
                  cTGGTCATCTCGACCTGTGGCTGTAACCAAATAGc                  <  1:5932291/35‑1 (MQ=255)
                   tGGTCATCTCGACCTGTGGCTGTAACCAAATAGCa                 <  1:734186/35‑1 (MQ=255)
                   tGGTCATCTCGACCTGTGGCTGTAACCAAATAGCa                 <  1:7091030/35‑1 (MQ=255)
                   tGGTCATCTCGACCTGTGGCTGTAACCAAATAGCa                 <  1:17975941/35‑1 (MQ=255)
                   tGGTCATCTCGACCTGTGGCTGTAACCAAATAGCa                 <  1:11583355/35‑1 (MQ=255)
                    ggTCATCTCGACCTGTGGCTGTAACCAAATAGCaa                <  1:12112526/35‑1 (MQ=255)
                    ggTCATCTCGACCTGTGGCTGTAACCAAATAGCaa                <  1:18194058/35‑1 (MQ=255)
                     gTCATCTCGACCTGTGGCTGTAACCAAATAGCaaa               <  1:9946625/35‑1 (MQ=255)
                     gTCATCTCGACCTGTGGCTGTAACCAAATAGCaaa               <  1:223107/35‑1 (MQ=255)
                     gTCATCTCGACCTGTGGCTGTAACCAAATAGCaaa               <  1:20534358/35‑1 (MQ=255)
                     gTCATCTCGACCTGTGGCTGTAACCAAATAGCaaa               <  1:15774585/35‑1 (MQ=255)
                     gTCATCTCGACCTGTGGCTGTAACCAAATAGCaaa               <  1:19322848/35‑1 (MQ=255)
                      tCATCTCGACCTGTGGCTGTAACCAAATAGCAAAc              <  1:3549401/35‑1 (MQ=255)
                      tCATCTCGACCTGTGGCTGTAACCAAATAGCAAAc              <  1:3404122/35‑1 (MQ=255)
                      tCATCTCGACCTGTGGCTGTAACCAAATAGCAAAc              <  1:4430932/35‑1 (MQ=255)
                        aTCTCGACCTGTGGCTGTAACCAAATAGCAAACtg            <  1:7547498/35‑1 (MQ=255)
                        aTCTCGACCTGTGGCTGTAACCAAATAGCAAACtg            <  1:13617555/35‑1 (MQ=255)
                        aTCTCGACCTGTGGCTGTAACCAAATAGCAAACtg            <  1:19711532/35‑1 (MQ=255)
                         tctcGCCCTGTGGCTGTAACCAAATAGCAAACtga           <  1:19174706/35‑1 (MQ=255)
                         tctcGACCTGTGGCTGTAACCAAATAGCAAACtga           <  1:12911418/35‑1 (MQ=255)
                         tctcGACCTGTGGCTGTAACCAAATAGCAAACtga           <  1:6443229/35‑1 (MQ=255)
                         tctcGACCTGTGGCTGTAACCAAATAGCAAACtga           <  1:20249172/35‑1 (MQ=255)
                         tcgcGACCTTTGGCTGTAACCAAATAGCAAACtga           <  1:20119364/32‑1 (MQ=37)
                          ctcGACCTGTGGCTGTAACCAAATAGCAAACtgat          <  1:18405826/35‑1 (MQ=255)
                          ctcGACCTGTGGCTGTAACCAAATAGCAAACtgat          <  1:16591456/35‑1 (MQ=255)
                          ctcGACCTGTGGCTGTAACCAAATAGCAAACtgat          <  1:7126201/35‑1 (MQ=255)
                           tcGCCCTGTGGCTGTAACCAAATAGCAAACtgatg         <  1:16090885/35‑1 (MQ=255)
                           tcGACCTGTGGCTGTAACCAAATAGCAAACtgatg         <  1:1113706/35‑1 (MQ=255)
                             gcccTGTGGCTGTAACCAAATAGCAAACTgatgat       <  1:9115695/33‑1 (MQ=255)
                              aCCTGTGGCTGTAACCAAATAGCAAACTGATGAtt      <  1:14977675/35‑1 (MQ=255)
                               ccTGTGGCTGTAACCAAATAGCAAACTGATGAttt     <  1:8501577/35‑1 (MQ=255)
                               ccTGTGGCTGTAACCAAATAGCAAACTGATGAttt     <  1:19019820/35‑1 (MQ=255)
                               ccTGTGGCTGTAACCAAATAGCAAACTGATGAttt     <  1:19086871/35‑1 (MQ=255)
                               ccTGTGGCTGTAACCAAATAGCAAACTGATGAttt     <  1:19717513/35‑1 (MQ=255)
                                 tgtgGCTGTAACCAAATAGCAAACTGATGATTTTc   <  1:16307725/35‑1 (MQ=255)
                                 tgtgGCTGTAACCAAATAGCAAACTGATGATTTTc   <  1:15910992/35‑1 (MQ=255)
                                 tgtgGCTGTAACCAAATAGCAAACTGATGATTTTc   <  1:18311357/35‑1 (MQ=255)
                                 tgtgGCTGTAACCAAATAGCAAACTGATGATTTTc   <  1:12474249/35‑1 (MQ=255)
                                  |                                  
GCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCC  >  REL606/3613783‑3613851

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