New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3643087 =	22 (0.760)	4 (0.150)	4/236	3.2	14.9%	coding (46/753 nt)	rsmJ	16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
	?	NC_000913	3643107 =	25 (0.930)					coding (26/753 nt)	rsmJ	16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAA  >  NC_000913/3642940‑3643208                                                                                                                                                    |                                                                                                                          ccATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAg                                                                                                                          >  3:119249/1‑149 (MQ=255) ccATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAg                                                                                                                          <  3:93916/149‑1 (MQ=255)                                  cgccAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAg                                                                                         <  5:33975/149‑1 (MQ=255)                                   gccAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCC‑‑AACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAg                                                                                         >  5:89145/1‑146 (MQ=255)                                       aGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGAt                                                                                    <  4:119249/149‑1 (MQ=255)                                                      cGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCC‑‑AACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAg                                                                   <  6:89145/149‑1 (MQ=255)                                                       gCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCTCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAAtt                                                                                            <  8:22552/125‑1 (MQ=255)                                                       gCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAAtt                                                                                            >  7:22552/1‑125 (MQ=255)                                                                                               cATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAAt                                                                                             >  3:115518/1‑84 (MQ=255)                                                                                               cATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAAt                                                                                             <  4:115518/84‑1 (MQ=255)                                                                                                     gTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCa                      >  4:106336/1‑149 (MQ=255)                                                                                                               tCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACgg             >  7:61746/1‑148 (MQ=255)                                                                                                                  tGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCt                                                               <  8:2845/95‑1 (MQ=255)                                                                                                                  tGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCt                                                               >  7:2845/1‑95 (MQ=255)                                                                                                                     tCCAGCCGCCAGCGGGCCGCCAGGACGGATAAGGCACGGGCTGCGGTGCCTGTGTCAGCAATTAAGCAGATTTTCACTGAGTGGTCGGCCCTTAGTGCCGTAATGCTCCAGCATCGCATCCGGCTGCGGGTCAGGACCAGGGAAGCGtg      >  8:85982/1‑148 (MQ=255)                                                                                                                     tCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGtt      >  1:97156/1‑149 (MQ=255)                                                                                                                       cAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTACGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGttt     <  6:119792/148‑1 (MQ=255)                                                                                                                        aGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGa   <  8:61746/149‑1 (MQ=255)                                                                                                                         gCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCTGTTCACGACCACGGAAGCGTTTGaa  <  7:85982/149‑1 (MQ=255)                                                                                                                         gCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGaa  <  3:106336/149‑1 (MQ=255)                                                                                                                                                gATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGttt     >  7:16841/1‑123 (MQ=255)                                                                                                                                                gATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGttt     <  8:16841/123‑1 (MQ=255)                                                                                                                                                    |                                                                                                                          CCATCGCTCCGCCAACAAAATCAACAAAGATGCCGCCAAGTTTTGGCTCATCACGCTTGCGCAATTCCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAA  >  NC_000913/3642940‑3643208

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