New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3867992	23 (1.130)	3 (0.160)	3/282	2.1	11.4%	coding (265/333 nt)	yidQ	DUF1375 family outer membrane protein
	?	NC_000913	= 3868005	25 (1.300)					coding (278/333 nt)	yidQ	DUF1375 family outer membrane protein
Rejected: Frequency below cutoff threshold.

TGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGACCATCCGGCGAGAAGACCAC  >  NC_000913/3867862‑3868139                                                                                                                                   |                                                                                                                                                    tGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCcaacgcg                                                                                                                                 >  1:123448‑M2/1‑144 (MQ=255)                          aCTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCAc                                                                                                        <  1:111801/151‑1 (MQ=255)                          aCTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCAc                                                                                                        <  1:21670/151‑1 (MQ=255)                                   cTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGGCTTCGCTTTTTTcaacgcgcgagcgcaccgagctgtctttgcggaacacatc                                                                                                >  1:8363‑M2/1‑110 (MQ=255)                                    tGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAAt                                                                                              <  2:166863/151‑1 (MQ=255)                                       cGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGaaaaaagcgaagccaacgcgcgagcgcaccgagctgtctttgcggaacacatcccacg                                                                                           <  1:75121‑M1/151‑59 (MQ=255)                                                 tatGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGaaaaaagcgaagccaacgcgcgagcgcaccgagctgtctttgcggaacacatcccacggcagtaaaag                                                                                 <  2:63008‑M1/151‑69 (MQ=255)                                                      cGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAg                                                                            <  2:237757/151‑1 (MQ=255)                                                       gTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAg                                                                            >  2:211021/1‑150 (MQ=255)                                                       gTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAg                                                                            >  1:249774/1‑150 (MQ=255)                                                             cGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAg                                                                     <  1:188569/151‑1 (MQ=255)                                                                 gtgATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGttt                                                                  >  1:266243/1‑150 (MQ=255)                                                                    aTGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCgg                                                               <  2:208713/150‑1 (MQ=255)                                                                     tGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTg                                                             <  1:270240/151‑1 (MQ=255)                                                                     tGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTg                                                             <  2:28384/151‑1 (MQ=255)                                                                       gACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGa                                                            <  2:157114/150‑1 (MQ=255)                                                                           cGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCca                                                        >  1:60605/1‑150 (MQ=255)                                                                           cGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCca                                                        >  1:140479/1‑150 (MQ=255)                                                                             cTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGGCTTCGCTTTTTTcaacgcgcgagcgcaccgagctgtctttgcggaacacatcccacggcagtaaaagcgtgtccatcacagcggtaaacggcata                                                     <  2:63008‑M2/151‑84 (MQ=255)                                                                              ttttACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGc                                                    <  2:121227/151‑1 (MQ=255)                                                                              ttttACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGc                                                    <  2:212808/151‑1 (MQ=255)                                                                                      ccGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGgccgc                                             >  2:65634/1‑150 (MQ=255)                                                                                       cGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGGCTTCGCTTTTTTcaacgcgcgagcgcaccgagctgtctttgcggaacacatcccacggcagtaaaagcgtgtccatcacagcggtaaacggcatatcgagaatcg                                           <  1:75121‑M2/151‑94 (MQ=255)                                                                                            gATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGaaaaaagcgaagccaacgcgcgagcgcaccgagctgtctttgcggaacacatcccacggcagtaaaagcgtgtccatcacagcggtaaacggcatatcgagaatcgccag                                       >  1:8363‑M1/1‑40 (MQ=255)                                                                                             aTGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTtccatc                                     <  1:164821/151‑1 (MQ=255)                                                                                              tgtgTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTtccatcc                                    <  1:70830/151‑1 (MQ=255)                                                                                                  ttCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACgc                                <  2:290378/151‑1 (MQ=255)                                                                                                                    gTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGACCATCCgg              >  2:95638/1‑151 (MQ=255)                                                                                                                             cgcgTTGgcttcgcttttttcaacgcgcgagcgcaccgagctgtctttgcggaacacatcccacggcagtaaaagcgtgtccatcacagcggtaaacggcatatcgagaatcgccagcgacttagtgccccagttggtttcatcgtcgcca     >  1:123448‑M1/1‑7 (MQ=255)                                                                                                                               cgTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGACCATCCGGCGAGAAGAcc    <  1:134650/150‑1 (MQ=255)                                                                                                                                 ttGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGACCATCCGGCGAGAAGACCAc  <  2:31987/150‑1 (MQ=255)                                                                                                                                   |                                                                                                                                                    TGGCGACGATGAAACCAACTGGGGCACTAAGTCGCTGGCGATTCTCGATATGCCGTTTACCGCTGTGATGGACACGCTTTTACTGCCGTGGGATGTGTTCCGCAAAGACAGCTCGGTGCGCTCGCGCGTTGAAAAAAGCGAAGCCAACGCGCAAGCCACCAATGCGGTTATTCCACCGGCCAGAATGCCTGATAACTGATTTAGCGATCAGTTTCGGTGATCCACAGCTGGCCGCCTTCCATCCACGCCAGCCACTGACCATCCGGCGAGAAGACCAC  >  NC_000913/3867862‑3868139

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