New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1900904 =	20 (1.060)	5 (0.270)	5/240	NT	22.5%	coding (1413/4026 nt)	PP_1703	assimilatory nitrate reductase/sulfite reductase
	?	NC_002947	1900952 =	15 (0.820)					coding (1461/4026 nt)	PP_1703	assimilatory nitrate reductase/sulfite reductase

CCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATG  >  NC_002947/1900756‑1901046                                                                                                                                                     |                                                                                                                                               aggccactggcgcggctgtgtccgcaattcggcatagccgataccgcgcaggtcataaccggttgaggggttgtggaagcggcgtatttcctcatacaccgcttcggcattggggtagtcgaaggcgtcggcatagcccatggcacacGc                                                                                                                                               >  2:202863‑M1/149‑150 (MQ=255)             ttACCCAGGACGCATTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCaaa                                                                                                                                   <  2:117501/150‑1 (MQ=255)                  cAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCAGACATCGTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTg                                                                                                                              >  1:139251/1‑150 (MQ=255)                     gACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCg                                                                                                                           <  2:105736/150‑1 (MQ=255)                                   cacCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCAt                                                                                                             >  2:508720/1‑150 (MQ=255)                                            caaccGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCgcgt                                                                                                    >  2:92724/1‑150 (MQ=255)                                              gcgcaggtcataaccggttgaggggttgtggaagcggcgtatttcctcatacaccgcttcggcattggggtagtcgaaggcgtcggcatagcccatggcacacGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCgcgtgg                                                                                                  <  1:53230‑M1/47‑1 (MQ=255)                                                   ggtcataaccggttgaggggttgtggaagcggcgtatttcctcatacaccgcttcggcattggggtagtcgaaggcgtcggcatagcccatggcacacGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCgtgt                                                                                             >  2:53230‑M1/99‑150 (MQ=255)                                                      gCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGtgtg                                                                                            >  1:152088/1‑148 (MQ=255)                                                      gCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGtgtg                                                                                            <  2:152088/148‑1 (MQ=255)                                                          accggttgaggggttgtggaagcggcgtatttcctcatacaccgcttcggcattggggtagtcgaaggcgtcggcatagcccatggcacacGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCTTGGCGTGTGCCATgg                                                                                      >  2:33492‑M1/92‑150 (MQ=255)                                                                                     ccGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCa                                                           >  1:352629/1‑150 (MQ=255)                                                                                              aGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAg                                                  <  2:352629/150‑1 (MQ=255)                                                                                                      atgatCAACAGCGAGCGCAACCTGACCCTCATGCCACGCGCTGTCGAGCCCCCCGGTCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTAt                                          >  1:140367/1‑150 (MQ=255)                                                                                                            aaCAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGaa                                    <  1:302867/150‑1 (MQ=255)                                                                                                                  gagcgCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGAt                                                                                                                <  2:220293/68‑1 (MQ=255)                                                                                                                  gagcgCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGAt                                                                                                                >  1:220293/1‑68 (MQ=255)                                                                                                                                  tGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAAccccg              <  2:140367/150‑2 (MQ=255)                                                                                                                                      gCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGc                                                                                                                             <  2:162186/35‑1 (MQ=255)                                                                                                                                      gCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGc                                                                                                                             >  1:162186/1‑35 (MQ=255)                                                                                                                                        cACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGc                                                                                           <  2:16242/67‑1 (MQ=255)                                                                                                                                        cACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGc                                                                                           >  1:16242/1‑67 (MQ=255)                                                                                                                                         cccatggcacacGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCAtt                                                                                                            <  1:361783‑M1/37‑1 (MQ=255)                                                                                                                                         cccatggcacacGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCAtt                                                                                                            >  2:361783‑M1/13‑49 (MQ=255)                                                                                                                                          cgcgCTGTCGAGCGCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCAGGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGTGGAAATACGCCGCTTCCACATCCCCTCAACCGGt      <  2:139251/150‑1 (MQ=255)                                                                                                                                              cTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATg  <  1:508720/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               CCGAACTGGTCATTACCCAGGACGCCTTCCTCGACACCGAAACCAACCGTTACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATG  >  NC_002947/1900756‑1901046

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