New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5821859 =16 (0.590)6 (0.250) 5/258 2.3 32.4% intergenic (+40/+89) PP_5092/pilT NLP/P60 protein/twitching motility protein
?NC_002947 5821922 = 11 (0.460)intergenic (+103/+26) PP_5092/pilT NLP/P60 protein/twitching motility protein

AGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCA  >  NC_002947/5821721‑5822006
                                                                                                                                                |                                                                                                                                                   
ggaggggcggttggcgcagttgtcttcagtgatgcagggcggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctttcgcgggcgCGCCCg                                                                                                                                                <  1:163828‑M1/6‑1 (MQ=255)
 gTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCGTGCCCGcgaaggg                                                                                                                                         <  1:2206/150‑7 (MQ=255)
          ttggcgcagttgtcttcagtgatgcagggcggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctttcgcgggcgCGCCCGCTCCCACAgg                                                                                                                                      >  2:410764‑M1/135‑150 (MQ=255)
           cTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGcgcg                                                                                                                               >  1:452933/1‑150 (MQ=255)
           cTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCCCTCCCACAGGTAGcgcg                                                                                                                               >  1:498493/1‑150 (MQ=255)
                              cAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTc                                                                                                            <  1:34694/150‑1 (MQ=255)
                                gAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGCAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTcac                                                                                                          <  1:522710/150‑1 (MQ=255)
                                gAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTcac                                                                                                          <  1:253135/150‑1 (MQ=255)
                                     ggcggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctttcgcgggcaCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTca                                                                                                           >  2:9466‑M1/108‑150 (MQ=255)
                                       cggggcggcggcggggtcgctgacgagggagggggggttgcggagggtgagggagcgagggttgatcataggcttgtagtacctgcaaaggccccttcgcgggcaCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTcaca                                                                                                         >  1:207186‑M1/106‑150 (MQ=255)
                                       cggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccccttcgcgggcaCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTcaca                                                                                                         >  1:207180‑M1/106‑150 (MQ=255)
                                          gtggtgGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                                <  2:102132/150‑1 (MQ=255)
                                           gcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctttcgcgggcgCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATAcc                                                                                                     >  2:163828‑M1/102‑150 (MQ=255)
                                               ggCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCTGCAAAGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgtg                                                                                                 <  1:298328/150‑1 (MQ=255)
                                                ggcggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctttcgcgggcgCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                                >  1:155507‑M1/97‑150 (MQ=255)
                                                   aaGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                                       <  1:264777/102‑1 (MQ=255)
                                                   aaGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                                       >  2:264777/1‑102 (MQ=255)
                                                              gcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGGTCCCACAGGTAGCGCGTCGCTGTGGCGCGGAGCTCACATACCTGTTGGTGCGGGGGGGCCCGCGaaa                                                                            >  2:37106/1‑150 (MQ=255)
                                                                        cagcaggTGAGTCTGGCCAAGCTTTAAGCTTATCGGCGCTATCCGTTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCCCGCCCGCTCCCCCAGGTAGCGCGTCGCTCTCGCGCGAAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGc                                                                  <  1:168931/150‑1 (MQ=38)
                                                                                   tCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCGCCCACAGGTAGCGCCTCGCGCTCGCGCGCAGCGCACAAACCGGTGGGAGCGGGCGCGCCCCCGCAAGGGGCTTTTCAGGGACTACaa                                                       >  2:260460/1‑150 (MQ=255)
                                                                                    gctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccccttcgcgggcaCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGtac                                                            >  1:396768‑M1/61‑150 (MQ=255)
                                                                                                 aaGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAAcc                                         <  1:455870/150‑1 (MQ=255)
                                                                                                             gcTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCa                             <  2:338782/150‑1 (MQ=255)
                                                                                                                                              cACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCa  <  1:452129/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACC‑‑‑‑‑‑GGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCA  >  NC_002947/5821721‑5822006

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