New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2947467	15 (0.370)	5 (0.190)	5/330	1.7	38.5%	intergenic (+4/‑30)	metZ/metW	tRNA‑Met/tRNA‑Met
	?	NC_000913	2947688 =	6 (0.220)					intergenic (+5/+69)	metV/amiC	tRNA‑Met/N‑acetylmuramoyl‑L‑alanine amidase C

ACCGCCGCGTTTAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCA‑GTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTT  >  NC_000913/2947225‑2947713                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        aCCGCCGCGTTTAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                                                                                                                                                                                 <  2:237148/251‑1 (MQ=255)  ccgccgCGTTTAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                                                                                                                                >  1:57969/1‑251 (MQ=255)     ccgagTTTAATGAAATGAGCAGACTCTCCCTGACGCAAATTGTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTGGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAT‑CAAATTTGATGAAg                                                                                                                                                                                                                                            <  2:94870/247‑1 (MQ=255)                       gcagcCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCA‑GTACGGTGAc                                                                                                                                                                                                                          >  1:226750/1‑251 (MQ=255)                           ccTCGCCCTGACGCAAATTTTGCACAAAAAATAGGTTTTAGTGATTTGTTTTTGTTCAAAATCATGCCTAAGCGGTGATCGGGGTAAAAAAAAGGTGGCATGAAAGCGCGATCGGAGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGGCTGGTAGTTCGTCGGGCGCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAGTTGATGAAGTAAAGCA‑GTACGTTGACGCgg                                                                                                                                                                                                                      <  2:379035/251‑1 (MQ=25)                                 ccTGACGCAAAATTTGCACAAAAAATATGCTATATTGATTTGTGTTTTTTCAAAATCAGGCCAAATCCGTGATCGGGATAAAAAAAAGGATGCATGAAAGCGCGAGCTGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTTGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCA‑GTACGGTGACGCTGGGTGGa                                                                                                                                                                                                                <  2:65011/251‑1 (MQ=25)                                                                          gtgtttgttcaaAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACACGAGCGGAGTATAGGGCGCGTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAttgaacaccctaacgggtgtttttttgtttctggtctcccataaaaaagcgccattcagcgcctttttatcatccccttc                                                                                                                                                                         <  2:108491‑M1/247‑81 (MQ=255)                                                                                           tGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGCTGAAGTAAAGCA‑GTACGGTGACGCGGGGTGGAGCAGACTGGTAGCTCGCCGGGCTCATAACCCGAAGGACGTCGGTTAAAATCCGGcccc                                                                                                                                                      >  1:368079/1‑251 (MQ=255)                                                                                                 aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCA‑GTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAc                                                                                                                                                >  1:211927/1‑251 (MQ=255)                                                                                                     gTGATGGGGGTAAAAGAAAGGTTGAACGAAACCGCGAGCGTGGTATAGTGCGAATCCGCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAttgaacaccctaacgggtgtttttttgtttctggtctcccataaaaaagcgccattcagcgcctttttatcatccccttctcgccagcgtcgccccatcggcaaaata                                                                                                                                             <  2:386705‑M1/251‑109 (MQ=255)                                                                                                     gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAttgaacaccctaacgggtgtttttttgtttctggtctcccataaaaaagcgccattcagcgcctttttatcatccccttctcgccagcgtcgccccatcggcaaaat                                                                                                                                              >  1:400620‑M1/1‑143 (MQ=255)                                                                                                               tAAAAAAAAGGTTGCATGAAAACGCGAGCGGCGTATAGTGCGCATCCACGGACGCGGGGCGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATATTTGCTGAAGTAAAGCA‑GTACGGTGACACGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCACATAACCCGCAGGTCGTCGGTACAACTCCGGCCCCCGCAACCACTCAAATTtgat                                                                                                                                  >  1:181313/1‑251 (MQ=18)                                                                                                               tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAGTTTGATGAAGTAAAGCA‑GTACGGTGTCGCGGGGTGGAGCAGCCTGGCAGCTCTTCGGGCTCATAACCCGAAGGCCGTCGGTTCAAATCCCGCCCCCGCCAACAATCactttagac                                                                                                                                  >  1:7063/1‑243 (MQ=18)                                                                                                                       atgTTGCATGAAAACGCGAGCGGAGTATAGTGCGCTTCCACGGACGCGGGGTGGAGCAGTCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATGTGATGAATGAAAGCA‑GTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaa                                                                                                                          <  1:370993/249‑1 (MQ=18)                                                                                                                                         aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAttgaacaccctaacgggtgtttttttgtttctggtctcccataaaaaagcgccattcagcgcctttttatcatccccttctcgccagcgtcgccccatcggcaaaatacgctttaatccccgcaagaatagactccgcaacttc                                                                                                         >  1:3929‑M1/1‑107 (MQ=255)                                                                                                                                                 ataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCA‑GTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATCAATTTAGATCCCGTAAAAGGAGTACggggatctggtgtggagc                                                                                                >  1:393018/1‑235 (MQ=14)                                                                                                                                                                ggACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAttgaacaccctaacgggtgtttttttgtttctggtctcccataaaaaagcgccattcagcgcctttttatcatccccttctcgccagcgtcgccccatcggcaaaatacgctttaatccccgcaagaatagactccgcaacttcctgctgtaaagtcgccgttttca                                                                                  <  2:226750‑M1/251‑168 (MQ=255)                                                                                                                                                                                 gtctggtatgtcggcgggattaTAAGCCGAATGTCGTCGTGTCATAACCGGCCTCCGCAACAAATTAAGCTTTGAAGAAGTAAAGCATGT‑CGGTGACGCGGGTTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCTCGCAACCAATCAAATTTGATGAAGTAAAAGCATTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGaa                                                                <  1:408832/230‑1 (MQ=14)                                                                                                                                                                                                                                        cGCAACCAATTAAAATTTGATGAAGTAAAGCA‑GTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACTCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTCTTGAACACCCCAACGGgtg         >  1:319554/1‑251 (MQ=255)                                                                                                                                                                                                                                               aaTTAAAATTTGATGAAGTAAAGCA‑GTACGGTGACGCGGGGTGGAGCAGCCGGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt  <  1:353289/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        ACCGCCGCGTTTAATGAAATGAGCAGCCTCTCCCTGACGCAAATTTTGCACAAAAAATAGGCTTTAGTGATTTGTTTTTGTTCAAAATCATGCCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCA‑GTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTT  >  NC_000913/2947225‑2947713

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