New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 291401	NA (NA)	6 (0.190)	6/254	3.2	15.8%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 2787809	32 (0.990)					pseudogene (1413/2253 nt)	ygaQ	uncharacterized protein; putative enzyme
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GAAAATATGGCGTACGAGGTTTACGTTATGATGCGGCAAAACACTCAAAACATGAACAAATAGAAAGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGCCATTGATTTTTCATTGCTCATGAAAACGATAAAAGCCTTTT  >  NC_000913/2787673‑2787951                                                                                                                                         |                                                                                                                                               gAAAATATGGCGTACGAGGTTTACGTTATGATGCGGCAAAACACTCAAAACATGAACAAATAGAAAGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTACCTTATTTAACCCAAAATATCATAAAAAAGcc                                                                                                                                    >  1:154339/1‑149 (MQ=255)      tatGGCGTACGAGGTTTACGTTATGATGCGGCAAAACACTCAAAACATGAACAAATAGAAAGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTAt                                                                                                                               <  3:138858/149‑1 (MQ=255)                gAGGTTTACGTTATGATGCGGCAAAACACTCAAAACATGAACAAATAGAAAGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATg                                                                                                                     >  1:174849/1‑149 (MQ=255)                                 gCGGCAAAACACTCAAAACATGAACAAATAGAAAGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTg                                                                                                    >  4:132900/1‑149 (MQ=255)                                   ggCAAAACACTCAAAACATGAACAAATAGAAAGATCAATAACCCCACCTATTAAAAATTATAATGATCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTa                                                                                                                                <  5:12945/119‑1 (MQ=255)                                    gCAAAACACTCAAAACATGAACAAATAGAAAGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCa                                                                                                 <  2:154339/149‑1 (MQ=255)                                           acTCAAAACATGAACAAATAGAAAGATCAATAACCCCACAGCTTAAAAATTATAATGAGAGGTTACACAATACTAACTTATTTAACCCAAAATATAATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGAt                                                                                          <  5:23849/149‑1 (MQ=255)                                               aaaaCATGAACAAATAGAAAGATCAATAACCCCACCGCTTAAAAATTATAGTGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGaac                                                                                      <  8:6674/149‑1 (MQ=255)                                               aaaaCATGAACAAATAGAAAGATCAATAACCCCACCGCTTAAAAATTATAGTGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGaa                                                                                       >  7:6674/1‑148 (MQ=255)                                                      gAACAAATGGAAAGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAAt                                                                               >  5:51350/1‑149 (MQ=255)                                                             tggAAAGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGtcatca                                                                        <  6:51350/147‑1 (MQ=255)                                                                aaaGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGa                                                                                           >  8:94395/1‑127 (MQ=255)                                                                aaaGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGa                                                                                           <  7:94395/127‑1 (MQ=255)                                                                     tCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTa                                                                >  1:159052/1‑149 (MQ=255)                                                                        aTAACCCCACCTCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATg                                                              <  7:33601/148‑1 (MQ=255)                                                                        aTAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATg                                                              >  8:33601/1‑148 (MQ=255)                                                                           aCCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGaaa                                                           <  2:159052/148‑1 (MQ=255)                                                                             cccACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGaaa                                                           >  1:33292/1‑146 (MQ=255)                                                                              tatattcTTAAAAATTATAATGAGCTTTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAgaga                                                       <  5:64357/143‑1 (MQ=255)                                                                                    cTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGTAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAgaga                                                       >  6:64357/1‑143 (MQ=255)                                                                                     ttAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTAt                                                               <  8:91672/134‑1 (MQ=255)                                                                                     ttAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTAt                                                               >  7:91672/1‑134 (MQ=255)                                                                                                  tGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGCCATTGAtt                                   <  2:174849/149‑1 (MQ=255)                                                                                                     gCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGCCATTGAttttt                                <  2:33292/149‑1 (MQ=255)                                                                                                       ggTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGCCATTGATTTTTCa                              >  3:81974/1‑149 (MQ=255)                                                                                                                      aCTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTa                                                                >  7:131368/1‑100 (MQ=255)                                                                                                                      aCTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTa                                                                <  8:131368/100‑1 (MQ=255)                                                                                                                          aTTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGCCATTGATTTTTCATTGCTCATGAAAACGATa            >  1:154277/1‑148 (MQ=255)                                                                                                                              aaCCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGcc                                          <  6:6513/114‑1 (MQ=255)                                                                                                                              aaCCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGcc                                          >  5:6513/1‑114 (MQ=255)                                                                                                                               aCCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGCCATTGATTTTTCATTGCTCATGAAAACGATAAAAGcc      <  3:132900/149‑1 (MQ=255)                                                                                                                                   aaaaTATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGCCATTGATTTTTCATTGCTCATGAAAACGATAAAAGCCtttt  >  4:75178/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                               GAAAATATGGCGTACGAGGTTTACGTTATGATGCGGCAAAACACTCAAAACATGAACAAATAGAAAGATCAATAACCCCACCGCTTAAAAATTATAATGAGCGGTTACACAATACTAACTTATTTAACCCAAAATATCATAAAAAAGCCGTTATGAATTACATGGAATATCTGGTAACTTGTCAGTTGGATGAACAACAAATGTCATCACTGCTTTATGAAAGAGATGATTTAAGCGCCATTGATTTTTCATTGCTCATGAAAACGATAAAAGCCTTTT  >  NC_000913/2787673‑2787951

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.