New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2290113	NA (NA)	3 (0.120)	3/486	2.8	12.5%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	2421265 =	21 (0.820)					coding (645/648 nt)	yfcG	disulfide bond oxidoreductase YfcG
Rejected: Frequency below/above cutoff threshold.

TACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCC  >  NC_000913/2289898‑2290346                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          tACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGa                                                                                                                                                                                                         >  1:65092/1‑250 (MQ=255)        cTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGAAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaa                                                                                                                                                                                                  >  1:233708/1‑250 (MQ=255)            gggTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGCCAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaaa                                                                                                                                                                                                >  2:235512/1‑247 (MQ=255)                    cgtcTCGCTGACTTCCACCATGGTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTGCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTTATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TGGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAAc                                                                                                                                                                                     <  2:84720/248‑1 (MQ=255)                             gACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTTACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGAt                                                                                                                                                                            <  2:65092/251‑1 (MQ=255)                                          atttGCCATGGCAGAATATGCTCCATGCGAGAACAGACAATCTCTTTTCTGGTCTGACGGCGATTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAATCCTATTCTATGCCTCCAGATGACAACCATGATCTACTATCAGG‑ACATTGTTACCACCTTCAATAACGCTGCAGCCAACAGACCAGATATTTCGGAACGACAAAGAGAAACAAAACACGATTGCAGCATTGacatttt                                                                                                                                                               >  2:282215/2‑245 (MQ=9)                                                             gCTCCATGCGGGACAAGAAAATCTCTTTTCTGTTCTGACGGCGCTTACTGCTGACTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCTCCTTCCATCACGCTGTGGCCACCAGAACCGATATTGAGGAACGACAATGCGAAACAGAACCAGATTGATGCATTAAGCTTTCATCCTATTAAATTAATTGc                                                                                                                                            >  2:33901/1‑251 (MQ=18)                                                                      gggACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTaaa                                                                                                                                     >  1:11727/1‑249 (MQ=255)                                                                          cAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAtt                                                                                                                               <  1:3175/251‑1 (MQ=255)                                                                                                                            tggCGAAGGTAAGATTATGACTCATGATGATTCCCGTTCTATGGCCTCATATGAAAAAAATGTTCTTATATCAGGGAC‑TTGTTCGCACCTTCCATATCGCTATAGCCACCAGAAAATATATTGCGGAACTAAAGAGAGAATCAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAGTGGGTACAGAAAATACCCATAGCTCCATACACGGAGTCAGTTTTTaataa                                                                              <  2:163820/249‑4 (MQ=255)                                                                                                                                     tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCATATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTaa                                                                      >  1:101322/1‑249 (MQ=255)                                                                                                                                               aCTCATTATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGttatt                                                                                   >  2:202488/1‑223 (MQ=255)                                                                                                                                               aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGttttt                                                                                   <  1:202488/226‑1 (MQ=255)                                                                                                                                                                          ccAGATGACAAACAGGATCTCAGAGTAGGGAC‑TTGTTCGCACCGTCCATAACGCTGTAGCCACCAGAACTGATATGGCGGAACGACAAAGAGAAACAGTACCAGATTGATGCATTGAGCTTGCGCCCTATGAAATTAATTGCGGTGAAAAGCATTGGGTACAGAAAGTGCCCATAGCTCCATACCCGGAGTCAGTTTTGAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGttt                               <  1:181334/251‑1 (MQ=255)                                                                                                                                                                                 aCAAACATGATCCCATATCAGGGAC‑TTGTGCGCACCGTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAGTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTtga                        <  2:11727/251‑1 (MQ=255)                                                                                                                                                                                  cAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTtgat                       <  1:204063/251‑1 (MQ=255)                                                                                                                                                                                    aaCATGATCTCATATCAGGTAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCATTTTTTAAAAACTGTTTAAAGGAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTtgatga                     <  2:233708/251‑1 (MQ=255)                                                                                                                                                                                      cATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGGTGTTTTTCTTGATgaga                   <  1:233875/251‑1 (MQ=255)                                                                                                                                                                                                       ggc‑TTGTTCGAACCTACATTAACGCTGTAGCCAACAGAACAGATATTGCTGAACGACAAAGAGAAACAGAACCAGATTTATGTATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATATCTCCATACCCGGAGTCAGTGTTTAAAAACTGTTTAAAGAAATTCACAAGTATTGTGATTTATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAATATTcc  <  2:222512/249‑1 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          TACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAC‑TTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCC  >  NC_000913/2289898‑2290346

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