New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 1899908	31 (0.780)	4 (0.110)	4/314	NT	14.3%	coding (459/4026 nt)	PP_1703	assimilatory nitrate reductase/sulfite reductase
	?	AE015451	= 1899959	19 (0.510)					coding (510/4026 nt)	PP_1703	assimilatory nitrate reductase/sulfite reductase

ACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCCGCGGCGCAGCGCAACCGCCGACAAGGCCGACCTGTTTCTGCAGGTGCGCCCCGGCAGTGACATGGCCTTGCTCAACGGCCTTTTGCACCTGCTGCATCGCAACGGCCACACCAACCCCGATTTCATCGCCCGCCACACCG  >  AE015451/1899739‑1900191                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         aCGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACAt                                                                                                                                                                                                                            <  1:730570/235‑1 (MQ=255)      ccGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGa                                                                                                                                                                                                                      <  1:649802/236‑1 (MQ=255)                  gcgcTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGatcggcgccaagcgacagcttgtagcca                                                                                                                                                                                                              >  1:279100‑M2/1‑204 (MQ=255)                  gcgcTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGa                                                                                                                                                                                                           >  1:452717/1‑235 (MQ=255)                       gTACGTGGCCGGGCAAAGGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGatcggcgccaagc                                                                                                                                                                                                                             >  2:732436‑M2/1‑199 (MQ=255)                                     aaaTGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGcc                                                                                                                                                                                            >  1:729595/1‑231 (MQ=255)                                     aaaTGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTg                                                                                                                                                                                          >  2:661347/1‑233 (MQ=255)                                      aaTGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCtt                                                                                                                                                                                       <  1:694234/235‑1 (MQ=255)                                                                      aGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGCGCACGCCGCCGTATTTCTGGTGAGCGGCGCCAACATGGCCGACTGCCACCCGCTc                                                                                                                                                                                                         >  2:105575/1‑185 (MQ=255)                                                                            ccAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGAc                                                                                                                                                <  1:690641/236‑1 (MQ=255)                                                                             cAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCgg                                                                                                                                                                      <  1:539119/213‑1 (MQ=255)                                                                                                                              aTGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGCTCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCcg                                                                                                                                             >  2:165762/1‑189 (MQ=255)                                                                                                                              aTGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCCGCGGCGCAGCGCAACCGCCGACAAGGCCGACCTGTTTCTGCAGGTgcg                                                                                              >  1:605793/1‑236 (MQ=255)                                                                                                                                   cagcgccagcAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCCGCGGCGCAGCGCAACCGCCGACAAGGCCGACCTGTTTCTGCAGGTGCGCCCCg                                                                                         <  1:661348/236‑1 (MQ=255)                                                                                                                                                         cTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCCGCGGCGCAGCGCAACCGCCGACAAGGCCGACCTGTTTCTGCAGGTGCGCCCCGGCAGTGACATGGCCTTGCTCa                                                                    <  1:265761/235‑1 (MQ=255)                                                                                                                                                              gCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCCAGGCCGGTGCGAAGCTgct                                                                                                                                                          >  2:92332/1‑142 (MQ=255)                                                                                                                                                                   gcgcCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCCGCGGCGCAGCGCAACCGCCGACAAGGCCGACCTGttt                                                                                                         >  2:661827/1‑188 (MQ=255)                                                                                                                                                                        gATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGatcggcgccaagcgacagcttgtagccactgctggcgctggccatgcacaggcgcgagttggactcgatgtggcgcgtacggataaagcccttggccaacttgttggccaggtactgggcctccagcgacatttgcccggacacgtacagcgccacggcatccgggccgtgctggtcgat                                                      <  2:534009‑M2/234‑181 (MQ=255)                                                                                                                                                                            gCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGatcggcgccaagcgacagcttgtagccactgatggcgctggccatgcacaggcgcgagttggactcgatgtggcgcgtacggataaagcccttggccaacttgttggccaggtactgggcctccagcgacatttgcccggacacgtacagcgccacggcatccgggccgtgctggtcgatgatgc                                                 <  1:602125‑M2/235‑186 (MQ=255)                                                                                                                                                                                  cTGGCAGTTACCAGGCCTTCGAGCCCGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCaa                                                                                                                                                                            <  2:229665/106‑1 (MQ=255)                                                                                                                                                                                             cAGGACTTCTAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCCGCGGCGCAGCGCAACCGCCGACAAGGCCGACCTGTTTCTGCAGGTGCGCCCCGGCAGTGACATGGCCTTGCTCAACGGCCTTTTGCACCTGCTGCATCGCAACGGCca                                  <  2:408158/233‑1 (MQ=255)                                                                                                                                                                                                        gCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCCGCGGCGCAGCGCAACCGCCGACAAGCCCGAc                                                                                                               >  2:119761/1‑145 (MQ=255)                                                                                                                                                                                                         cACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCCGCGGCGCAGCGCAACCGCCGACAAGGCCGACCTGTTTCTGCAGGTGCGCCCCGGCAGTGACATGGCCTTGCTCAACGGCCTTTTGCACctg                                                   >  2:500490/1‑204 (MQ=255)                                                                                                                                                                                                         cACGCCGACGTATTGCTGGTGATCGGCGCCACCATGGCCGACTGCCCCCCGATCCTGTTCCTGCGCCTGCTTGATCgg                                                                                                                                                                                 <  2:514824/78‑1 (MQ=255)                                                                                                                                                                                                               gACGTATTTCTGGTGatcggcgccaagcgacagcttgtagccactgctggcgctggccatgcacaggcgcgagttggactcgatgtggcgcgtacggataaagcccttggccaacttgttggccaggtactgggcctccagcgaca                                                                                                       <  2:92654‑M2/146‑132 (MQ=255)                                                                                                                                                                                                                           gTGatcggcgccaagcgacagcttgtagccactgctggcgctggccatgcacaggcgcgagttggactcgatgtggcgcgtacggataaagcccttggccaacttgttggccaggtactgggcctccagcgacatttgcccggacacgtacagcgccacggcatccgggccgtgctggtcgatgatgccgcgcagacgctcggcggtggcagcgatggcctggtccatgctggtg  >  1:534010‑M2/1‑3 (MQ=255)                                                                                                                                                                                                                           gTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGACCGGGTCAAGGCCGGTGCGAAGCTGAtcgt                                                                                                                                                       >  2:338392/1‑86 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         ACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTGGTGATCGGCGCCAACATGGCCGACTGCCACCCGATCCTGTTCCTGCGCCTGCTTGATCGGGTCAAGGCCGGTGCGAAGCTGATCGTCGTCGACCCGCGGCGCAGCGCAACCGCCGACAAGGCCGACCTGTTTCTGCAGGTGCGCCCCGGCAGTGACATGGCCTTGCTCAACGGCCTTTTGCACCTGCTGCATCGCAACGGCCACACCAACCCCGATTTCATCGCCCGCCACACCG  >  AE015451/1899739‑1900191

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