New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1186541	10 (0.310)	7 (0.220)	7/486	1.6	41.2%	coding (698/1227 nt)	pepT	peptidase T
	?	NC_000913	= 1469240	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

GGTGATATTCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTATCATCTGGCGAGCATGAAAGGCACCGTTGAACGGGCCGATATGCACTACATCATCCGTGATTTCGACCGTAAACAGTTTGAAGCGCGTAAACGTAAAATGATGGAGATCGCCAAAAAAGTGGGCAAAGGGTTACATCCTGA  >  NC_000913/1186330‑1186766                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  ggTGATATTCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGTGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGGGATGGTAAAtagacgggccccctgaatctccagacaaccaatatcact                                                                                                                                                                                            <  2:323764‑M1/251‑40 (MQ=255) ggTGATATTCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTcc                                                                                                                                                                                            <  1:177294/251‑1 (MQ=255)         tCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGAt                                                                                                                                                                                     <  1:49654/250‑1 (MQ=255)              tCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGAGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAAtagactggccccctgaatctccagacaaccactatcacttaaataagtgat                                                                                                                                                                                >  1:221933‑M1/1‑199 (MQ=255)                                               gggCGAAACATGTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGAGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAAtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtccttgg                                                                                                                                               >  2:190491‑M1/1‑166 (MQ=255)                                                  cGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTatcat                                                                                                                                           >  2:76526/1‑251 (MQ=255)                                                       cATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAAtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagat                                                                                                                                      <  1:96974‑M1/251‑94 (MQ=255)                                                                              tGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTATCATCTGGCGAGCATGAAAGGCACCGTTGAAc                                                                                                               >  1:114927/1‑251 (MQ=255)                                                                                   gCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTATCATCTGGCGAGCATGAAAGGCACCGTTGAACGGGcc                                                                                                          <  2:111982/251‑1 (MQ=255)                                                                                                                            gAGTTTGAAATCTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCGGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTATCATCTGGCGAGCATGAAAGGCACCGTTGAACGGGCCGATATGCACTACATCATCCGTGATTTCGACCGTAAACAGtt                                                                 <  1:127458/251‑1 (MQ=255)                                                                                                                            gAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTATCATCTGGCGAGCATGAAAGGCACCGTTGAACGGGCCGATATGCACTACATCATCCGTGATTTCGACCg                                                                          >  1:176513/1‑242 (MQ=255)                                                                                                                            gAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTATCATCTGGCGAGCATGAAAGGCACCGTTGAACGGGCCGATATGCACTACATCATCCGTGATTTCGACCg                                                                          <  2:176513/242‑1 (MQ=255)                                                                                                                            gAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTATCATCTGGCGAGCATGAAAGGCACCGTTGAACGGGCCGATATGCACTACATCATCCGTGATTTCGACCGTAAACAGtt                                                                 <  1:260663/251‑1 (MQ=255)                                                                                                                                                                   aTTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAAtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccg                          >  1:144203‑M1/1‑50 (MQ=255)                                                                                                                                                                          gTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAAtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggcaggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacat                   >  2:204877‑M1/1‑43 (MQ=255)                                                                                                                                                                                       cccgGGCACGGCGAAAGGAGTGATGGTAAAtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtgtagcagcca      <  2:144203‑M1/250‑222 (MQ=255)                                                                                                                                                                                           ggaaatataggggcaaatccaggaaaTGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTATCATCTGGCGAGCATGAAAGGCACCGTTGAACGTGCCGATATGCACTACATCATCCGTGATTTCGACCGTAAACAGTTTGAAGCGCGTAAACGTAAAATGATGGAGATCGCCAAAAAAGTGGGCAAAGGGTTACATCCTGa  <  1:176010/228‑1 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  GGTGATATTCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCATCCGGGCACGGCGAAAGGAGTGATGGTAAATGCGCTGTCGCTGGCGGCACGTATTCATGCGGAAGTTCCGGCGGATGAAAGCCCGGAAATGACAGAAGGCTATGAAGGTTTCTATCATCTGGCGAGCATGAAAGGCACCGTTGAACGGGCCGATATGCACTACATCATCCGTGATTTCGACCGTAAACAGTTTGAAGCGCGTAAACGTAAAATGATGGAGATCGCCAAAAAAGTGGGCAAAGGGTTACATCCTGA  >  NC_000913/1186330‑1186766

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