New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 526508	NA (NA)	3 (0.040)	3/260	NT	38.7%	coding (3248/4281 nt)	rhsD	protein RhsD
	?	NC_000913	= 526599	5 (0.060)					coding (3339/4281 nt)	rhsD	protein RhsD

AAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGG  >  NC_000913/526493‑526746                                                                                                           |                                                                                                                                                    aaaaaGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGcgccgctgcgctttttcccgctcgccgttctctgtctcgaccc                                                                                                          <  2:118349‑M2/150‑44 (MQ=255)               cgcAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGcgccgctgcgctttttcccgctcgccgttctctgtctcgacccggatgagcggcgtg                                                                                            <  1:1167966‑M2/150‑58 (MQ=255)                                                 gTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGcgccgctgcgctttttcccgctcgccgttctctgtctcgacccggatgagcggcgtgaagcttcccggctcgtataccgtctggatacgtg                                                          >  2:1167966‑M2/1‑59 (MQ=255)                                                                                tGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGc                           >  1:783730/1‑150 (MQ=11)                                                                                        tGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCttt                   <  2:661392/150‑1 (MQ=34)                                                                                           ggCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATc                <  2:93619/150‑1 (MQ=35)                                                                                                 tgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGAGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCg          <  2:475526/150‑1 (MQ=37)                                                                                                  ggACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGa         >  1:94237/1‑150 (MQ=37)                                                                                                     cAGGCTGcgccgctgcgctttttcccgctcgccgttctctgtctcgacccggatgagcggcgtgaagcttcccggctcgtataccgtctggatacgtgtggtgtcag                                                 <  1:1345117‑M2/107‑101 (MQ=255)                                                                                                     cAGGCTGcgccgctgcgctttttcccgctcgccgttctctgtctcgacccggatgagcggcgtgaagcttcccggctcgtataccgtctggatacgtgtggtgtcag                                                 >  2:1345117‑M2/1‑7 (MQ=255)                                                                                                      aGGCTGcgccgctgcgctttttcccgctcgccgttctctgtctcgacccggatgagcggcgtgaagcttcccggctcgtataccgtctggatacgtgtggtgtcagtctgcaccgtcgtcagcctgtctccgtcccagccataccacgtc     >  2:1398703‑M2/1‑6 (MQ=255)                                                                                                       ggCTGcgccgctgcgctttttcccgctcgccgttctctgtctcgacccggatgagcggcgtgaagcttcccggctcgtataccgtctggatacgtgtggtgtcagtctgcaccgtcgtcagcctgtctccgtcccagccataccacgtca    <  2:1303366‑M2/150‑146 (MQ=255)                                                                                                       ggCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGaccacc    >  2:945576/1‑150 (MQ=37)                                                                                                         cTGcgccgctgcgctttttcccgctcgccgttctctgtctcgacccggatgagcggcgtgaagcttcccggctcgtataccgtctggatacgtgtggtgtcagtctgcaccgtcgtcagcctgtctccgtcccagccataccacgtcacc  >  2:1358372‑M2/1‑3 (MQ=255)                                                                                                         cTGcgccgctgcgctttttcccgctcgccgttctctgtctcgacccggatgagcggcgtgaagcttcccggctcgtataccgtctggatacgtgtggtgtcagtctgcaccgtcgtcagcctgtctccgtcccagccataccacgtcacc  <  2:657321‑M2/150‑148 (MQ=255)                                                                                                         cTGcgccgctgcgctttttcccgctcgccgttctctgtctcgacccggatgagcggcgtgaagcttcccggctcgtataccgtctggatacgtgtggtgtcagtctgcaccgtcgtcagcctgtctccgtcccagccataccacgtcacc  >  2:883212‑M2/1‑3 (MQ=255)                                                                                                           |                                                                                                                                                    AAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAAGGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGG  >  NC_000913/526493‑526746

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