New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 603799	0 (0.000)	16 (0.600)	15/490	0.6	100%	coding (865/1248 nt)	ybdG	miniconductance mechanosensitive channel YbdG
	?	NC_000913	3720633 =	NA (NA)					noncoding (1/1443 nt)	IS150	repeat region

AAACCGTTATCACCTGGAGCCAGTTGGCGTACCATTAAGGTCATATCTTTACGAATCCGCGGATGGTTACGCAGATATTCGTTCAGATAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGCGCTTTATTCAGACGTTGCATTTCGTCTTCATCAAGAAAACGAATACTGGTGACATCAATACTGATACTGCGCTTAATACGTCGCCCGCCAGATGCTGACATCCCGCTCCAGTTTTTAAAGGAGTCAGAAACCAGAGACCAGGTGGGAATAGTGGTAATGGTATTGTCCCAGTTACGCACTTTGACGGTGGTTAACCCAATATCGATCACCGCGCCATCCGCGCCGTATTTCGGCATCTCC  >  NC_000913/603588‑604040                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  aaaCCGTTATCACCTGGAGCCAGTTGGCGTACCATTAAGGTCATATCTTTACGAATCCGCGGATGGTTACGCAGATATTCGTTCAGATAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaat                                                                                                                                                                                                              >  2:71821‑M1/1‑212 (MQ=255)            ccctGGAGCCAGTTGGCGTACCATTAAGGTCATATCTTTACGAATCCGCGGATGGTTACGCAGATATTCGTTCAGATAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatg                                                                                                                                                                                                 <  1:239034‑M1/250‑51 (MQ=255)               gggaGCCAGTTTGCGTACCATTAAGGTCATATCCTTACGAATCCGCGGATGGTTACGCAGATATTCGTTCAGATAGGCACGAAAAGTTCCAATATTGTTCAGTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACACATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTGTAACAAATGtgtactgcacccattttgtgggacgatgaaatggaatagcccctaatgtggat                                                                                                                                                                                              <  2:55096‑M1/250‑54 (MQ=255)                ggAGCCAGTTGGCGTACCATTAAGGTCATATCTTTACGAATCCGCGGATGGTTACGCAGATATTCGTTCAGATAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtca                                                                                                                                                                                              >  1:183803‑M1/1‑197 (MQ=255)                 gAGCCAGTTGGCGTACCATTAAGGTCATATCTATACGAATCCGAGGATGGTTACGCAGATATTCGTTCCGATAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAACTAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaa                                                                                                                                                                                             >  1:161981‑M1/1‑196 (MQ=255)                           gCGTACCATTAAGGTCATATCTTTACGAATCCGCGGATGGTTACGCAGATATTCGTTCAGATAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatac                                                                                                                                                                                  <  1:192513‑M1/251‑66 (MQ=255)                                                               aTGGTTACGCAGATATTCGTTCAGATAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactac                                                                                                                                              >  2:143854‑M1/1‑150 (MQ=255)                                                                 ggTTACGCAGATATTCGTTCAGATAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactactt                                                                                                                                            >  1:240543‑M1/1‑148 (MQ=255)                                                                                       aTAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacagg                                                                                                                      >  2:115759‑M1/1‑126 (MQ=255)                                                                                                    gTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcac                                                                                                         <  1:115759‑M1/251‑139 (MQ=255)                                                                                                                  cATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttgatgtccc                                                                                            <  1:71821‑M1/250‑152 (MQ=255)                                                                                                                                                              gATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagaccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtaaggccatgggttgccctctatgaaaaacatgg                                               >  1:88784‑M1/1‑55 (MQ=255)                                                                                                                                                                                             ggTTAAATAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctctatgaaaaacatggagaaaaaggtttaattcccaaacctaaaggc                <  2:32033‑M1/251‑228 (MQ=255)                                                                                                                                                                                                  aaTAAGGCTTTAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatgcccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacgtgatcatcgcggcacgttttggtgtcccccgaacccaggtcaggacatggtttgccctctatgaaaaacatggagaaaaaggtttaattcccaaacctaaaggcgttag           <  1:83220‑M1/251‑233 (MQ=255)                                                                                                                                                                                                          tttAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccataatatgtcaaagccaaaatacccttttgaaaagcgccttgaagacgtgaataactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctctatgaaaaacatggagaaaaaggtttaattcccaaacctaaaggcgttagtgctgatc   >  1:113891‑M1/1‑11 (MQ=255)                                                                                                                                                                                                            tAACAAATGtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtgtcccccgaacccaggtcaggacatgggttgccctctatgaaaaacatggagaaaaaggtttaattcccaaacctaaaggcgttagtgctgatcc  >  2:164128‑M1/1‑9 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  AAACCGTTATCACCTGGAGCCAGTTGGCGTACCATTAAGGTCATATCTTTACGAATCCGCGGATGGTTACGCAGATATTCGTTCAGATAGGCACGAAAGGTTCCAATATTGGTCATTCGGCGCAGATTTAATACCGACTCCGTAGAACCTTGCTGGCGATTCCACTCATTAATTTCCTGATGGCGGCTGGTTAAATAAGGCTTTAACAAATGCGCTTTATTCAGACGTTGCATTTCGTCTTCATCAAGAAAACGAATACTGGTGACATCAATACTGATACTGCGCTTAATACGTCGCCCGCCAGATGCTGACATCCCGCTCCAGTTTTTAAAGGAGTCAGAAACCAGAGACCAGGTGGGAATAGTGGTAATGGTATTGTCCCAGTTACGCACTTTGACGGTGGTTAACCCAATATCGATCACCGCGCCATCCGCGCCGTATTTCGGCATCTCC  >  NC_000913/603588‑604040

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