New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5200454	27 (0.910)	5 (0.180)	4/232	NT	16.5%	coding (597/774 nt)	yafJ	amidotransferase
	?	NC_002947	= 5200465	25 (0.900)					coding (586/774 nt)	yafJ	amidotransferase

AGGGCCCACTGACCCGGCGCGTAGCGGTGCCAGGTCTCGTTCTCGGTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAG  >  NC_002947/5200317‑5200591                                                                                                                                                       |                                                                                                                               aGGGCCCACTGACCCGGCGCGTAGCGGTGCCCGGTCACGT‑‑TAGTTCAGACGTGTGCTGTTGCAATCACCGTGACAACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGgc                                                                                                                               <  1:178636/150‑1 (MQ=255)       cACTGACCCGGCGCGTAGCGGTGCCAGGTCTCGTTCTCGGTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGcgcacgg                                                                                                                         <  2:42256‑M2/150‑8 (MQ=255)          tGACCCGGCGCGTAGCGGTGCCAGGTCTCGTTCTCGGTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCAcggcgg                                                                                                                      <  1:401787/150‑1 (MQ=255)                  cgcgTAGCGGTGCCAGGTCTCGTTCTCGGTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAAtgtg                                                                                                              <  1:55092/150‑1 (MQ=255)                   gcgTAGCGGTGCCAGGTCTCGTTCTCGGTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGc                                                                                                             <  2:476657/150‑1 (MQ=255)                   gcgTAGCGGTGCCAGGTCTCGTTCTCGGTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGc                                                                                                             <  1:42256/150‑1 (MQ=255)                            tGCCAGGTCTCGTTCTCGGTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCttt                                                                                                    <  1:111726/150‑1 (MQ=255)                            tGCCAGGTCTCGTTCTCGGTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCttt                                                                                                    <  1:47403/150‑1 (MQ=255)                                              gTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAAc                                                                                  >  1:233770/1‑150 (MQ=255)                                                    ggCCT‑GGGTTGCAATCACCGTGAGCACGTCGTTGGGGGTTGGTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCACCCGGGCGGCAGCGAAGGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCa                                                                             <  2:76861/150‑1 (MQ=255)                                                    ggCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAATCAGAACAGCCa                                                                             >  1:130244/1‑150 (MQ=255)                                                        t‑CGGTGGCAATCACCGTGACCACTTCGTGGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCg                                                                         >  2:212754/1‑150 (MQ=255)                                                                 aaTCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGTACCAGCTTTGTCGa                                                                                               >  1:111468/1‑120 (MQ=255)                                                                 aaTCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGTACCAGCTTTGTCGa                                                                                               <  2:111468/120‑1 (MQ=255)                                                                      ccGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCa                                                            >  1:298138/1‑150 (MQ=255)                                                                      ccGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCa                                                            >  1:298288/1‑150 (MQ=255)                                                                             cACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGtt                                                     <  2:661167/150‑1 (MQ=255)                                                                                     tGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCATTTGAACAccc                                             <  2:298288/150‑1 (MQ=255)                                                                                     tGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACAccc                                             <  2:298138/150‑1 (MQ=255)                                                                                      gggggTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGc                                                         <  1:380738/137‑1 (MQ=255)                                                                                      gggggTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGc                                                         >  2:380738/1‑137 (MQ=255)                                                                                       ggggTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACAccccc                                           <  1:212754/150‑1 (MQ=255)                                                                                                          aaTCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAGTGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGcac                        >  1:295309/1‑150 (MQ=255)                                                                                                           aTCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGAAGAGCACCCCCAGATCCCGGTAGCCGGcaca                       >  1:599498/1‑150 (MQ=255)                                                                                                            tCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAg                      <  2:233770/150‑1 (MQ=255)                                                                                                                   tCAGTTCGACATCTTTCATCCGTGCGGCTCCGAATGGCGCACGGCGGGTCATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCCTGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGa               <  2:130244/150‑1 (MQ=255)                                                                                                                     aGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGAcc             >  1:46987/1‑150 (MQ=255)                                                                                                                          gACATCCTTCAACCGTGCGGCACCGAGTGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCAc        <  2:295309/150‑1 (MQ=255)                                                                                                                                cTTCAACCGTGCGGCACCGAATGGcgcacggttgaaggatgtcgacctgatcgtcgattttcataccgaaaccacccccaacgacgtggtcacggtgattgccaccgaggccctgaccgagaacgagacctggcaccgctacgcgccggg  <  1:367881‑M2/150‑127 (MQ=255)                                                                                                                                cTTCAACCGTGCGGCACCGAATGGcgcacggttgaaggatgtcgacctgatcgtcgattttcataccgaaaccacccccaacgacgtggtcacggtga                                                      <  2:334484‑M2/98‑75 (MQ=255)                                                                                                                                cTTCAACCGTGCGGCACCGAATGGcgcacggttgaaggatgtcgacctgatcgtcgattttcataccgaaaccacccccaacgacgtggtcacggtga                                                      >  1:334484‑M2/1‑24 (MQ=255)                                                                                                                                    aaCCGTGCGGCACCGAATGGcgcacggttgaaggatgtcgacctgatcgtcgattttcataccgaaaccacccccaacgacgtggtcacggtgattg                                                   >  1:599538‑M2/1‑20 (MQ=255)                                                                                                                                    aaCCGTGCGGCACCGAATGGcgcacggttgaaggatgtcgacctgatcgtcgattttcataccgaaaccacccccaacgacgtggtcacggtgattg                                                   <  2:599538‑M2/97‑78 (MQ=255)                                                                                                                                                       |                                                                                                                               AGGGCCCACTGACCCGGCGCGTAGCGGTGCCAGGTCTCGTTCTCGGTCAG‑GGCCT‑CGGTGGCAATCACCGTGACCACGTCGTTGGGGGTGGTTTCGGTATGAAAATCGACGATCAGGTCGACATCCTTCAACCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAG  >  NC_002947/5200317‑5200591

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