New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4285004 =	31 (0.820)	5 (0.140)	5/244	3.1	13.7%	coding (210/315 nt)	yjcH	DUF485 family inner membrane protein
	?	NC_000913	4285020 =	33 (0.910)					coding (194/315 nt)	yjcH	DUF485 family inner membrane protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GAAAGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTT  >  NC_000913/4284861‑4285152                                                                                                                                                |                                                                                                                                                     gAAAGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCTGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAac                                                                                                                                                 <  2:99610/149‑1 (MQ=255)       gAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAAt                                                                                                                                           <  7:102455/149‑1 (MQ=255)                                 ctctTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGgtgt                                                                                                                 >  1:3881/1‑149 (MQ=255)                                         atgatgCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGAtt                                                                                                          >  5:63619/1‑148 (MQ=255)                                              gCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTc                                                                                                               >  3:144963/1‑138 (MQ=255)                                              gCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTc                                                                                                               <  4:144963/138‑1 (MQ=255)                                                          tGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTgg                                                                                                                    <  6:130222/121‑1 (MQ=255)                                                          tGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTgg                                                                                                                    >  5:130222/1‑121 (MQ=255)                                                           gCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCccagccag                                                                                       >  6:166918/1‑149 (MQ=255)                                                             aGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCgt                                                                                                                              <  6:36010/108‑1 (MQ=255)                                                             aGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCgt                                                                                                                              >  5:36010/1‑108 (MQ=255)                                                                  ttCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGgcg                                                                                >  8:8771/1‑149 (MQ=255)                                                                   tCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGgcgc                                                                               >  1:15879/1‑149 (MQ=255)                                                                     attattAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGcgcg                                                                              >  2:207627/1‑148 (MQ=255)                                                                         ttaAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAgg                                                                          <  6:63619/148‑1 (MQ=255)                                                                              aCGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCa                                                                    <  2:3881/149‑1 (MQ=255)                                                                              aCGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCa                                                                    <  7:8771/149‑1 (MQ=255)                                                                                 gTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGt                                                                  >  3:59709/1‑148 (MQ=255)                                                                                  tCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTa                                                                 >  7:128305/1‑148 (MQ=255)                                                                                        tCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATa                                                            >  4:119778/1‑147 (MQ=255)                                                                                         cGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATaaa                                                          <  8:128305/148‑1 (MQ=255)                                                                                          gCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGc                                                        <  4:59709/149‑1 (MQ=255)                                                                                             gTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGa                                                     <  2:15879/149‑1 (MQ=255)                                                                                                     gccAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACt                                             <  5:166918/149‑1 (MQ=255)                                                                                                         gATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGc                                           >  8:160453/1‑147 (MQ=255)                                                                                                                      gTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTg                                                                                                <  7:70936/81‑1 (MQ=255)                                                                                                                      gTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTg                                                                                                >  8:70936/1‑81 (MQ=255)                                                                                                                      gTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGa                            >  6:46714/1‑149 (MQ=255)                                                                                                                          gaACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAg                                                                           <  7:45756/98‑1 (MQ=255)                                                                                                                          gaACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAg                                                                           >  8:45756/1‑98 (MQ=255)                                                                                                                                              aCTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGcc     <  1:207627/148‑1 (MQ=255)                                                                                                                                                tCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGtt  >  6:76804/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     GAAAGGGAGTGTGGCGGCAAGCGCCGTCAGAACTCTTTTCATGATGCTTGTACCTCATGCAGGACTTCATTATTAAGACGGTCGAATTCGCCGTTCGCCCGCCAGATGTAGATACCGGTGAGAACAAAGGAGATCACAATCACTCCAACACCAATTGGAATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTT  >  NC_000913/4284861‑4285152

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