New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4559490 =	45 (1.410)	6 (0.200)	5/248	2.9	12.9%	coding (37/1173 nt)	iadA	isoaspartyl dipeptidase
	?	NC_000913	4559512 =	38 (1.240)					coding (15/1173 nt)	iadA	isoaspartyl dipeptidase
Rejected: Frequency below/above cutoff threshold.

GGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGCAGACCGTAATGATATGCGGGTAATATTTGGCATTCACTTCGGCTGTC  >  NC_000913/4559349‑4559635                                                                                                                                              |                                                                                                                                                  ggATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCTTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTACCTGCAGCAgg                                                                                                                                            <  2:98168/149‑1 (MQ=255) ggATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAgg                                                                                                                                            <  2:100677/149‑1 (MQ=255)   aTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAggg                                                                                                                                           >  2:38248/1‑148 (MQ=255)           aCTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgcagc                                                                                                                                               <  8:107139/136‑1 (MQ=255)           aCTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgcagc                                                                                                                                               >  7:107139/1‑136 (MQ=255)                 aTCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTa                                                                                                                                         <  1:146511/136‑1 (MQ=255)                 aTCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTa                                                                                                                                         >  2:146511/1‑136 (MQ=255)                            tGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATaatcaatc                                                                                                                 <  1:38248/149‑1 (MQ=255)                            tGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCACAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTa                                                                                                                                         <  1:113617/125‑1 (MQ=255)                            tGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCACAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTa                                                                                                                                         >  2:113617/1‑125 (MQ=255)                                  gccaaatattacccgcatatcattacggtctgcgtgattaacgcattgctgtcgatctgggttatgcagcttattgtttaacaaggagttaccatgattgattataccGCAGCAGGGTAAAACCGGCTGCGGt                                                                                                                           <  4:71193‑M1/25‑1 (MQ=255)                                  gccaaatattacccgcatatcattacggtctgcgtgattaacgcattgctgtcgatctgggttatgcagcttattgtttaacaaggagttaccatgattgattataccGCAGCAGGGTAAAACCGGCTGCGGt                                                                                                                           >  3:71193‑M1/109‑133 (MQ=255)                                      tACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGAGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTc                                                                                                       >  7:42976/1‑149 (MQ=255)                                       aCAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGt                                                                                                            >  8:72240/1‑143 (MQ=255)                                       aCAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGt                                                                                                            <  7:72240/143‑1 (MQ=255)                                               aGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAatcaatca                                                                                                                <  6:87442/131‑1 (MQ=255)                                               aGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAatcaatca                                                                                                                >  5:87442/1‑131 (MQ=255)                                                     gATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTg                                                                                                   >  6:93891/1‑138 (MQ=255)                                                     gATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTg                                                                                                   <  5:93891/138‑1 (MQ=255)                                                      atatTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAg                                                                                       >  3:107504/1‑149 (MQ=255)                                                             tGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATAGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCAt                                                                                 >  1:151970/1‑148 (MQ=255)                                                               gCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAAc                                                                              <  4:107504/149‑1 (MQ=255)                                                                     gattaacgcattgctgtcgatctgggttatgcagcttattgtttaacaaggagttaccatgattgattataccGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGtt                                                                                                 <  4:68028‑M1/51‑1 (MQ=255)                                                                     gattaacgcattgctgtcgatctgggttatgcagcttattgtttaacaaggagttaccatgattgattataccGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGtt                                                                                                 >  3:68028‑M1/74‑124 (MQ=255)                                                                       gATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATaat                                                                                                                      >  1:134216/1‑101 (MQ=255)                                                                       gATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATaat                                                                                                                      <  2:134216/101‑1 (MQ=255)                                                                       gATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATc                                                                       <  8:102869/148‑1 (MQ=255)                                                                       gATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATc                                                                       >  7:102869/1‑148 (MQ=255)                                                                        aTAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGa                                                                     >  3:77687/1‑149 (MQ=255)                                                                             tttGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTcc                                                                                                      <  3:103499/111‑1 (MQ=255)                                                                             tttGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTcc                                                                                                      >  4:103499/1‑111 (MQ=255)                                                                                   tgtcgatctgggttatgcagcttattgtttaacaaggagttaccatgattgattataccGCAGCAGGGTAAAACCGGCTgca                                                                                                                             <  1:57516‑M1/23‑2 (MQ=255)                                                                                   tgtcgatctgggttatgcagcttattgtttaacaaggagttaccatgattgattataccGCAGCAGGGTAAAACCGGCTgca                                                                                                                             >  2:57516‑M1/60‑81 (MQ=255)                                                                                   gTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAAcc                                                                             <  2:16833/130‑1 (MQ=255)                                                                                   gTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAAcc                                                                             >  1:16833/1‑130 (MQ=255)                                                                                            gaGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGc                                                 <  4:77687/149‑1 (MQ=255)                                                                                                   tCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGc                                                                                   >  4:82986/1‑108 (MQ=255)                                                                                                   tCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGc                                                                                   <  3:82986/108‑1 (MQ=255)                                                                                                   tCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGCAGACCg                                           >  2:36632/1‑148 (MQ=255)                                                                                                     gCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGCa                                                <  5:149367/141‑1 (MQ=255)                                                                                                     gCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGCa                                                >  6:149367/1‑141 (MQ=255)                                                                                                      cAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGCAGACCGTAAt                                       <  8:42976/149‑1 (MQ=255)                                                                                                             cccGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGCAGACCGTAATGATATGc                                <  1:36632/149‑1 (MQ=255)                                                                                                             cccGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGCAGACCGTAATGATATGc                                >  2:37498/1‑149 (MQ=255)                                                                                                              ccGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAg                                                                  >  1:113999/1‑114 (MQ=255)                                                                                                              ccGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAg                                                                  <  2:113999/114‑1 (MQ=255)                                                                                                                    ttCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAAc                                                                                             >  2:105879/1‑81 (MQ=255)                                                                                                                    ttCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAAc                                                                                             <  1:105879/81‑1 (MQ=255)                                                                                                                                 aaaTGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCTTTAATCACGCAGACCGTAATGATATGCGGGTAATATTTGGCATTCAc            <  2:151970/149‑1 (MQ=255)                                                                                                                                 aaaTGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGCAGACCGTAATGATATGCGGGTAATATTTGGCATTCAc            <  1:37498/149‑1 (MQ=255)                                                                                                                                           ccTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCGTGTTAAACAATAAGCTGCATAACCCAGCTCGACAGCAATGCGTTAATCACGCAGACCGTAATGATATGCGGGTAATATTTGGCATTCACTTCGGCTGTc  >  3:97873/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                  GGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGTAAAACCGGCTGCGGTATAATCAATCATGGTAACTCCTTGTTAAACAATAAGCTGCATAACCCAGATCGACAGCAATGCGTTAATCACGCAGACCGTAATGATATGCGGGTAATATTTGGCATTCACTTCGGCTGTC  >  NC_000913/4559349‑4559635

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