New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3598659	25 (0.550)	4 (0.090)	4/242	NT	12.0%	coding (1089/1521 nt)	PP_3175	ferredoxin reductase dioxygenase subunit
	?	NC_002947	= 3598684	34 (0.770)					coding (1064/1521 nt)	PP_3175	ferredoxin reductase dioxygenase subunit

GCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTCCAGGCCTTCGATCCGCTCGACCTCGGTGGGGCCGTGGAAGGTGA  >  NC_002947/3598535‑3598807                                                                                                                             |                                                                                                                                                     gcCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCgg                                                                                                                             >  1:257636/1‑150 (MQ=255)      tGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGcc                                                                                                                                                 >  1:1180122/1‑125 (MQ=255)      tGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGcc                                                                                                                                                 <  2:1180122/125‑1 (MQ=255)              gcgcACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGt                                                                                                                <  2:384994/150‑1 (MQ=255)                                 aGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCCCCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAgggg                                                                                             <  1:1028447/150‑1 (MQ=255)                                    gggAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTc                                                                                          >  1:243097/1‑150 (MQ=255)                                        aCGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCa                                                                                      <  1:906737/150‑1 (MQ=255)                                                             cccACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGtt                                                                 >  2:278357/1‑150 (MQ=255)                                                              ccACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCccgccaccgcgttccttcaaggcgtgcc                                                                                                                          <  2:518911‑M1/92‑29 (MQ=255)                                                              ccACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCccgccaccgcgttccttcaaggcgtgcc                                                                                                                          >  1:518911‑M1/1‑64 (MQ=255)                                                                acaAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCa                                                              >  1:83587/1‑150 (MQ=255)                                                                 caAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCccgccaccgcgttccttcaaggcgt                                                                                                                             <  2:129568‑M1/86‑26 (MQ=255)                                                                 caAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCccgccaccgcgttccttcaaggcgt                                                                                                                             >  1:129568‑M1/1‑61 (MQ=255)                                                                 caAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCAc                                                             <  1:278357/150‑1 (MQ=255)                                                                    aCCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTAGTCAGTTGCACCGc                                                          >  1:128174/1‑150 (MQ=255)                                                                       gTCAGCTGCGCGCCTTTCTTCGTAGATCCTCACTGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTc                                                       <  2:63062/150‑1 (MQ=255)                                                                        tCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTg                                                                        <  1:1188603/132‑1 (MQ=255)                                                                        tCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTg                                                                        >  2:1188603/1‑132 (MQ=255)                                                                          aGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGAc                                                    <  2:243097/150‑1 (MQ=255)                                                                                gcgcATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTc                                              >  2:231468/1‑150 (MQ=255)                                                                                         ttcgtcgtTCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTCCAGGCCTTc                                     <  2:257636/148‑1 (MQ=255)                                                                                                  ccTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCg                                                                                                            >  1:444950/1‑70 (MQ=255)                                                                                                  ccTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCg                                                                                                            <  2:444950/70‑1 (MQ=255)                                                                                                  ccTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTCCAGGCCTTCGATCCGCTc                            >  1:277497/1‑150 (MQ=255)                                                                                                   cTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGa                                                                                                     <  2:950708/76‑1 (MQ=255)                                                                                                   cTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGa                                                                                                     >  1:950708/1‑76 (MQ=255)                                                                                                    tCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGcc                                                                                   >  2:258655/1‑93 (MQ=255)                                                                                                    tCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGcc                                                                                   <  1:258655/93‑1 (MQ=255)                                                                                                            ccgatctCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTCCAGGCCTTCGATCCGCTCGACCTCGGTg                  <  1:861083/144‑1 (MQ=255)                                                                                                                    gCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTCCAGGCCTTCGATCCGCTCGACCTCGGTg                  >  2:120665/1‑142 (MQ=255)                                                                                                                    gCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTCCAGGCCTTCGATCCGCTCGACCTCGGTg                  <  1:120665/142‑1 (MQ=255)                                                                                                                           ggCACGCCTTGAAGGAACGCGGTGGCGGgccgcttggcgaagacaaatcggtgagggtcgacgcagaaatgcgcgcagctgacggtttgtgg                                                             >  1:518911‑M2/1‑28 (MQ=255)                                                                                                                           ggCACGCCTTGAAGGAACGCGGTGGCGGgccgcttggcgaagacaaatcggtgagggtcgacgcagaaatgcgcgcagctgacggtttgtgg                                                             <  2:518911‑M2/92‑65 (MQ=255)                                                                                                                            gCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAgg                                                                                               >  1:776277/1‑57 (MQ=255)                                                                                                                            gCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAgg                                                                                               <  2:776277/57‑1 (MQ=255)                                                                                                                            gCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTCCAGGCCTTCGATCCGCTCGACCTCGGTGGGGCCGTGGAAGGTGa  >  1:727838/1‑150 (MQ=255)                                                                                                                            gCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTCCAGACCTTCGATCCGCTCGACCTCGGTGGGGCCGTGGAAGGTGa  <  1:231468/150‑1 (MQ=255)                                                                                                                             |                                                                                                                                                     GCCAGTGCTCGATGCGCACCGGGCGCCCCGTGAGCGGGAACGTGGCGATATCGCCGGCGGCCCACAAACCGTCAGCTGCGCGCATTTCTGCGTCGACCCTCACCGATTTGTCTTCGCCAAGCGGCACGCCTTGAAGGAACGCGGTGGCGGGCCTTACCCCGGTGCCGAGCAGGACCAAGGGGGTCTGCAGCCGTTCACCATTGGTCAGTTGCACCGCTTCGACCTTGTCCAGGCCTTCGATCCGCTCGACCTCGGTGGGGCCGTGGAAGGTGA  >  NC_002947/3598535‑3598807

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