New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1411924	6 (0.130)	21 (0.490)	20/432	1.0	85.1%	coding (25/936 nt)	ttcA	tRNA cytosine(32) 2‑sulfurtransferase TtcA
	?	NC_000913	1434985 =	2 (0.050)					intergenic (‑579/+200)	ynaM/uspF	protein YnaM/universal stress protein F

TCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGCCGTCAACAGCCAGCAGGCAAAAGTATGC  >  NC_000913/1411692‑1412162                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               tCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgatttt                                                                                                                                                                                                                               <  2:319378‑M1/250‑18 (MQ=255)     ggAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGATGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgac                                                                                                                                                                                                                          <  2:346070‑M1/251‑23 (MQ=255)     ggAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgac                                                                                                                                                                                                                          <  1:456678‑M1/251‑23 (MQ=255)                    tttGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgact                                                                                                                                                                                                            <  2:335404‑M1/250‑37 (MQ=255)                    tttGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgac                                                                                                                                                                                                             <  1:238413‑M1/249‑36 (MQ=255)                     ttGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactc                                                                                                                                                                                                           <  2:417144‑M1/250‑38 (MQ=255)                                tAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTGCTTGCATTATTTCAGTTCTCTGGTACTAAAt                                                                                                                                                                                                <  1:160551/250‑1 (MQ=255)                                               aaaaaTTGATTGGCGCGCTTTGCTGCTAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGATGCAAACCATGATTCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGTCCACGTTACGATGCACACGTTTTTGTAATTTGTGCAGGTTGTATTGTTCTTTCTTcctaatttcttgattttgcgacttttaaaagcgactcaattcgttatatcgcatcagaagagta                                                                                                                                                                                <  2:281070‑M1/251‑65 (MQ=255)                                                      gattgGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGACAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcataagaagagtatgcgtc                                                                                                                                                                          <  2:184945‑M1/250‑71 (MQ=255)                                                          ggCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCATCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCaa                                                                                                                                                                     <  1:124078/251‑1 (MQ=255)                                                          ggCGCGCTTTGCTGCAAAATGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTTCGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTTTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcagcagaagagtatgcgtcatgcc                                                                                                                                                                     <  2:224906‑M1/251‑76 (MQ=255)                                                            cgcgcTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAACTTGGGGGCAAACTTTGCacct                                                                                                                                                                   >  1:400955/1‑248 (MQ=255)                                                               gcTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagtatgcgtcatgccggaac                                                                                                                                                                >  2:236509‑M1/1‑171 (MQ=255)                                                                    gctgcAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagtatgcgtcatgccggaacgccc                                                                                                                                                            >  1:388616‑M1/1‑166 (MQ=255)                                                                    gcngcAAATTGCGCAGAATCTCCAGCATGGTATAGCCGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagtatgcgtcatgccggaacgccc                                                                                                                                                            <  2:446159‑M1/250‑85 (MQ=255)                                                                                        tcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagtatgcgtcatgccggaacgcccagcataagaaatctgatata                                                                                                                                        >  1:29161‑M1/1‑146 (MQ=255)                                                                                         ccAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagtatgcgtcatgccggaacgcccagcataagaaatctggtataaa                                                                                                                                      <  2:300719‑M1/251‑107 (MQ=255)                                                                                                    taGCTGTCTTTACCACCGGAGAGGCAAACCAGGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGTGCAACCTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGtt                                                                                                                           >  2:188028/1‑251 (MQ=255)                                                                                                                cccccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagtatgcgtcatgccggaacgcccagcataagaaatctgatataaaaaactgtggcgtgtatggtacgg                                                                                                               <  1:363320‑M1/251‑130 (MQ=255)                                                                                                                  cccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagtatgcgtcatgccggaacgcccagcataagaaatctgatataaaaaactgtggcgggtgaggtacggct                                                                                                             >  2:31961‑M1/1‑120 (MQ=255)                                                                                                                    cGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAGTCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGGATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcataagagtatgcgtcatgcaggaacgcccatcataggcaatctgatataaaaaactgtggcgtgtagggtccggcttt                                                                                                           >  2:477945‑M1/1‑118 (MQ=255)                                                                                                                                         gATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagaatgcgtcatgccggaacgcccagcataagaaatctgatataaaaaactgtggcgtgtgtggtccggataaaaggagaaaatgatagcacat                                                                                      >  2:168803‑M1/1‑97 (MQ=255)                                                                                                                                                                  gTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagtatgcgtcatgccggaacgcccagcataagaaatctgatataaaaaactgtggcgtgtatggtacggattagaggggaaaatgtcagcacatttgcgaaatgaatcaaaaagcccgc                                                             <  1:470930‑M1/251‑180 (MQ=255)                                                                                                                                                                   tCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTactaatttcttgattttgcgacatttaaaagcgactcaattcgttatatggcatcagaagagtatgcgtcatgccggaacgcccagcataagaaatctgatataaaaaactgtggcgtgtatggtacggattagaggggaaaatgtcagcacatttgcgaaatgaatcaaaaagcccgc                                                             >  2:76487‑M1/1‑71 (MQ=255)                                                                                                                                                                                                   atggttttgtaaTTAGTTCAGGTTGTATTGTTCTTTCTTactaatgtcttgattgtgcgacatttaaaagcgacgcaagtcgttatatggcatcagaagagtatgcgtcatgccggaacgcccagcataagaaatctgatataaaaaactgtggcgtgtatggtacggattagaggggaaaatgtcagcacatttgcgaaatgaatcaaaaagcccgcagcaatgtgcgggcgttagtgtcagcgtacaac                            <  2:406796‑M1/247‑213 (MQ=255)                                                                                                                                                                                                     gTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGCCGTc                          >  2:252581/1‑251 (MQ=255)                                                                                                                                                                                                                             aTTGTTCTTTCTTTGTCATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGCCGTCAACAGCCAGCAGGCAAAAGTATGc  >  2:279220/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               TCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGCCGTCAACAGCCAGCAGGCAAAAGTATGC  >  NC_000913/1411692‑1412162

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.