New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4421693 =	38 (0.690)	4 (0.080)	4/240	NT	10.6%	coding (1443/1662 nt)	PP_3912	DNA‑cytosine methyltransferase
	?	NC_002947	4421750 =	32 (0.620)					coding (1500/1662 nt)	PP_3912	DNA‑cytosine methyltransferase

TCTGCGCGAAAACCCAGCGAGCGGTCCAGATGGCGTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACTT  >  NC_002947/4421606‑4421897                                                                                                                                                 |                                                                                                                                                    atcgccttgtagcgcggcccgtcagggcattcgtctgcaggcttgccgcgataggggatgcgggtgtggtcgccgggaaagccctggAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAgg                                                                                                                                                 >  4:69941‑M1/88‑149 (MQ=255)    gcgcGAAAACCCAGCGAGCGGTCCAGATGGCGTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCt                                                                                                                                              >  3:152481/1‑149 (MQ=255)       cgAAAACCCAGCGAGCGGTCCAGATGGCGTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTc                                                                                                                                           >  6:84394/1‑149 (MQ=255)        gAAAACCCAGCGAGCGGTCCAGATGGCGTTGGTGTGCAGGCCAATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTcc                                                                                                                                          <  8:36110/149‑1 (MQ=255)            aCCCAGCGAGCGGTCCAGATGGCGTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGc                                                                                                                                      <  7:279875/149‑1 (MQ=255)              ccAGCGAGCGGTCCAGATGGCGTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGa                                                                                                                                    >  1:110386/1‑149 (MQ=255)              ccAGCGAGCGGTCCAGATGGCGTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGa                                                                                                                                    >  6:201310/1‑149 (MQ=255)               cAGCGAGCGGTCCAGATGGCGTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGAc                                                                                                                                   <  4:152481/149‑1 (MQ=255)                                  cGTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATc                                                                                                                <  5:84394/149‑1 (MQ=255)                                   gTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATcg                                                                                                               >  5:85972/1‑149 (MQ=255)                                            cAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGcc                                                                                                      >  6:86076/1‑149 (MQ=255)                                              ggCCGATCACGCATACACATTGGAAGCACGGGCAGCAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCCCGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTg                                                                                                    >  2:250258/1‑149 (MQ=255)                                              ggCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTg                                                                                                    <  2:324589/149‑1 (MQ=255)                                                    tCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGAc                                                                                               >  7:211265/1‑148 (MQ=255)                                                                    gAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGcc                                                                              >  5:294584/1‑149 (MQ=255)                                                                       gCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCcgcg                                                                           <  7:334623/149‑1 (MQ=255)                                                                         aCGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCgcgc                                                                          >  3:269392/1‑148 (MQ=255)                                                                         aCGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCt                                                                         <  5:86076/149‑1 (MQ=255)                                                                              cAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAg                                                                    <  6:294584/149‑1 (MQ=255)                                                                                gAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGc                                                                  <  2:110386/149‑1 (MQ=255)                                                                                   gTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTACAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACatcgc                                                                  >  1:181033/1‑142 (MQ=255)                                                                                   cctggagccattcgcactcccgtggagtgaggcggcgcacctgcgcattggtctgaatggctCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGAt                                                               >  4:69941‑M2/63‑149 (MQ=255)                                                                                      cAAGCCATTCAGACCAATGCGCAGGTGCGCCGACTCACTCCACGGTAGTCCGAATGGCTCCACGGCTATCCCGGCGACCACACCCGCATCCACTATCGCGGCAAGCCTACAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCgg                                                            >  1:118631/1‑149 (MQ=255)                                                                                       aaGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCgg                                                            >  5:352064/1‑148 (MQ=255)                                                                                       aaGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCgg                                                            >  5:120103/1‑148 (MQ=255)                                                                                        agccattcgcactcccgtggagtgaggcggcgcacctgcgcattggtctgaatggctCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAgg                                                                   <  8:267239‑M2/83‑1 (MQ=255)                                                                                        agccattcgcactcccgtggagtgaggcggcgcacctgcgcattggtctgaatggctCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAgg                                                                   >  7:267239‑M2/58‑140 (MQ=255)                                                                                        AGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCa                                                          >  4:269392‑M1/1‑149 (MQ=255)                                                                                           CATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCTGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACa                                                       >  1:250258‑M1/1‑149 (MQ=255)                                                                                             TTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAg                                                      <  2:281906‑M1/148‑1 (MQ=255)                                                                                                 gACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAgg                                                 <  6:352064/149‑1 (MQ=255)                                                                                                   ccAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGcc                                               <  5:201310/149‑1 (MQ=255)                                                                                                     aaTGCGCAGGTGCGCCGCCTCACTCCTCGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGGGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGAGCTACAAGGCGATCGGCAACAGCAAGGCCgt                                             <  2:118631/149‑1 (MQ=255)                                                                                                          gcAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTtcg                                        <  2:109927/149‑1 (MQ=255)                                                                                                            aGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTtcgtc                                      <  8:211265/149‑1 (MQ=255)                                                                                                              gTGCGCCGCCTCACTCCACGGGAGTGCGAATGGATCCAGGGCTTTCCCGGCGACCACACCCGAATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTAAAAGGCGATCGGCAACAGCAAGGACGTGTtagtcg                                     >  6:213889/1‑148 (MQ=255)                                                                                                                  gcggcgcacctgcgcattggtctgaatggctCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGc                                <  3:69941‑M2/118‑1 (MQ=255)                                                                                                                               aCGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACACGCATCCCCTATCGCGGCCAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATcggccggc                   >  5:252737/1‑149 (MQ=255)                                                                                                                                 gggAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATccc                                                                                                                     >  3:118494/1‑49 (MQ=255)                                                                                                                                 gggAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATccc                                                                                                                     <  4:118494/49‑1 (MQ=255)                                                                                                                                        cGAATGGCGCCAGTGCTTTCCCGGCGACCACCCCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCcaa           <  6:252737/148‑1 (MQ=255)                                                                                                                                              gctCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCTCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTTTATCGGCCGGCGCATCCaacaacaac    <  5:213889‑M2/146‑1 (MQ=255)                                                                                                                                              gctCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCaacaacaac    >  1:128231‑M2/4‑149 (MQ=255)                                                                                                                                                tCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACtt  >  2:155762‑M2/2‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TCTGCGCGAAAACCCAGCGAGCGGTCCAGATGGCGTTGGTGTGCAGGCCGATCACGCATACACATTGGAAGCACGGGCAGAAGTCCAAGCCATTCAGACCAATGCGCAGGTGCGCCGCCTCACTCCACGGGAGTGCGAATGGCTCCAGGGCTTTCCCGGCGACCACACCCGCATCCCCTATCGCGGCAAGCCTGCAGACGAATGCCCTGACGGGCCGCGCTACAAGGCGATCGGCAACAGCAAGGCCGTGTTCGTCGTCCGCTGGATCGGCCGGCGCATCCAACAACAACTT  >  NC_002947/4421606‑4421897

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