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

AGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACGGTACTAAGGCGCTACGACAGTACGAAAATGGTGATCTAACTGCAGCCGGTTACACTGTG  >  NC_002947/4110188‑4110478                                                                                                                                                      |                                                                                                                                               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                                                                                                                            <  2:509542/150‑1 (MQ=255)                    cccTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGa                                                                                                                            <  1:532126/150‑1 (MQ=255)                     gtagcgccttagtaccgtccgcgaaaaaacccagcgcaatcgccgggtatccggctaaacgcgcaaacagattgctactcaaaaacttgataaccccatTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATggg                                                                                                                                            <  1:377320‑M1/34‑1 (MQ=255)                     gtagcgccttagtaccgtccgcgaaaaaacccagcgcaatcgccgggtatccggctaaacgcgcaaacagattgctactcaaaaacttgataaccccatTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATggg                                                                                                                                            >  2:377320‑M1/100‑133 (MQ=255)                     ccTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAg                                                                                                                           <  2:377644/150‑1 (MQ=255)                                        gtgtTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCCGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGt                                                                                                        <  2:37461/150‑1 (MQ=255)                                          gtTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGa                                                                                                                            >  1:102018/1‑128 (MQ=255)                                          gtTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGa                                                                                                                            <  2:102018/128‑1 (MQ=255)                                                aaacccagcgcaatcgccgggtatccggctaaacgcgcaaacagattgctactcaaaaacttgataaccccattgccacagcgccgggcatggcctttggcaATGGGGTTATCAAGt                                                                                                                                  <  1:427614‑M2/15‑1 (MQ=255)                                                aaacccagcgcaatcgccgggtatccggctaaacgcgcaaacagattgctactcaaaaacttgataaccccattgccacagcgccgggcatggcctttggcaATGGGGTTATCAAGt                                                                                                                                  >  2:427614‑M2/103‑117 (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                                                                        >  2:427614‑M1/16‑117 (MQ=255)                                                                                                         acttgataaccccatTGCCAAAG‑GCCATGCCCGGCGCTGTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGttt                                                                        <  1:427614‑M1/102‑1 (MQ=255)                                                                                                                    cccattgccaaaggccatgcccggcgctgtggcaATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACGGTACTAAGGCGCTAc                                              >  2:377320‑M2/35‑133 (MQ=255)                                                                                                                    cccattgccaaaggccatgcccggcgctgtggcaATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACGGTACTAAGGCGCTAc                                              <  1:377320‑M2/99‑1 (MQ=255)                                                                                                                                              cTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACGGTACTAAGGCGCTACGACAGTACGAAAATGGTGATCTAACTGCAGCCGGTTACACtgt   <  1:579299/150‑1 (MQ=255)                                                                                                                                               tttGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGCTGCTTCATGGACGGTACTAAGGCGCTACGACAGTACGCAAATGGTGATCTAGTTGCAGCGGTTTACACtgtg  >  1:167275/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                               AGCGATGCACCGGAATACACCCTGGGGTTCGCTTCATCTGTGTTCGGGGTGATCGGTGCAGCGGCGGCAACGTTGGTCAGCGTACGGGCAACGCAAAAGGCCGTCATGCTGAGATTAAGTGCCACAGCGCCGGGCATGGC‑CTTTGGCAATGGGGTTATCAAGTTTTTGAGTAGCAATCTGTTTGCGCGTTTAGCCGGATACCCGGCGATTGCGCTGGGTTTTTTCGCGGACGGTACTAAGGCGCTACGACAGTACGAAAATGGTGATCTAACTGCAGCCGGTTACACTGTG  >  NC_002947/4110188‑4110478

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