New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4110307 =	29 (1.240)	4 (0.180)	4/244	NT	14.9%	coding (2085/2991 nt)	PP_3616	membrane protein
	?	NC_002947	4110336 =	18 (0.800)					coding (2114/2991 nt)	PP_3616	membrane protein

TCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACGGTACTAAGGCGCTACGA  >  NC_002947/4110159‑4110436                                                                                                                                                     |                                                                                                                                   tCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTg                                                                                                                                   >  2:42552/1‑150 (MQ=255)               aacaGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGc                                                                                                                     <  1:465523/150‑1 (MQ=255)               aacaGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGc                                                                                                                     >  1:563281/1‑150 (MQ=255)                acaGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCa                                                                                                                    <  1:343800/150‑1 (MQ=255)                 caGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCAt                                                                                                                   <  2:516149/150‑1 (MQ=255)                    ccTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGc                                                                                                                <  1:508034/150‑1 (MQ=255)                    ccTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCAGCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGc                                                                                                                >  2:30972/1‑150 (MQ=255)                      tGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑Ct                                                                                                             >  2:575285/1‑150 (MQ=255)                           aagaGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGc                                                                                                        >  2:532126/1‑150 (MQ=255)                              agCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAAt                                                                                                     <  1:207925/150‑1 (MQ=255)                               gCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATg                                                                                                    <  2:610788/150‑1 (MQ=255)                                cGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATgg                                                                                                   <  2:492496/150‑1 (MQ=255)                                     cACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTa                                                                                              >  1:377644/1‑150 (MQ=255)                                      aCCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTAt                                                                                             >  2:10422/1‑150 (MQ=255)                                      aCCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTAt                                                                                             >  1:55509/1‑150 (MQ=255)                                               cacCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGttttt                                                                                    <  2:126137/150‑1 (MQ=255)                                                acCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAAAGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTg                                                                                   >  2:134084/1‑150 (MQ=255)                                                 cccTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGa                                                                                  <  1:532126/150‑1 (MQ=255)                                                 cccTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGa                                                                                  <  2:509542/150‑1 (MQ=255)                                                  gtagcgccttagtaccgtccgcgaaaaaacccagcgcaatcgccgggtatccggctaaacgcgcaaacagattgctactcaaaaacttgataaccccatTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATggg                                                                                                  >  2:377320‑M1/100‑133 (MQ=255)                                                  gtagcgccttagtaccgtccgcgaaaaaacccagcgcaatcgccgggtatccggctaaacgcgcaaacagattgctactcaaaaacttgataaccccatTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATggg                                                                                                  <  1:377320‑M1/34‑1 (MQ=255)                                                  ccTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAg                                                                                 <  2:377644/150‑1 (MQ=255)                                                                     gtgtTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCCGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGt                                                              <  2:37461/150‑1 (MQ=255)                                                                       gtTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGa                                                                                  <  2:102018/128‑1 (MQ=255)                                                                       gtTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGa                                                                                  >  1:102018/1‑128 (MQ=255)                                                                                            ggcggcAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGCGTAGCACTCTGTTCTCGCGTTTAGCCTTAGACCCGGCGATTgc                                       <  2:529403/150‑1 (MQ=255)                                                                                                 cAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTgg                                  <  1:133766/150‑1 (MQ=255)                                                                                                          tCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTcg                         <  1:575285/150‑1 (MQ=255)                                                                                                           cAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATCCCCGGCGATTGCGCTGGGTTTTTTcgc                        <  2:55509/150‑1 (MQ=255)                                                                                                                tACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACg                   <  1:30972/150‑1 (MQ=255)                                                                                                                            aaaaGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACGGTACTAAGgcgc       >  1:446316/1‑150 (MQ=255)                                                                                                                                 gCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACGGTACTAAGGCGCTACGa  <  1:134084/150‑1 (MQ=255)                                                                                                                                      acttgataaccccatTGCCAAAG‑GCCATGCCCGGCGCTGTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGttt                              <  1:427614‑M1/102‑1 (MQ=255)                                                                                                                                      acttgataaccccatTGCCAAAG‑GCCATGCCCGGCGCTGTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGttt                              >  2:427614‑M1/16‑117 (MQ=255)                                                                                                                                                     |                                                                                                                                   TCAAGACGGCCTACAACAGCCTGCAGAAGAGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACGGTACTAAGGCGCTACGA  >  NC_002947/4110159‑4110436

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