New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	23 (0.600)	7 (0.210)	7/252	3.2	20.5%	coding (290/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	1209619 =	34 (1.000)					pseudogene (1/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.

GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGC  >  NC_000913/1209481‑1209758                                                                                                                                           |                                                                                                                                            gATAAGCCTTTGCCAGTTCCGGTTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCAAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                                                         <  1:543721/143‑1 (MQ=255)          ttGTCAGTTTCGGGTATTCTTCAGCAGAAAAAGCGGCACCATGGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                                                                                                <  1:282166/143‑1 (MQ=255)            gCCAGTTCCGGGTATTCTTCAGCATAACAAGCGGCACCCTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGTCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAGGAAAATgaa                                                                                                                              <  2:52413/143‑1 (MQ=255)            gCCAGTTCCGGGTAGTCTTCAGCAGAAAAAGCGGCATCTTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                              <  1:190090/143‑1 (MQ=255)                   ccGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATaa                                                                                                                       <  2:211436/143‑1 (MQ=255)                      ggTATTCTTCAGCAGAAAAAGCCGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                                                    <  1:118992/143‑1 (MQ=255)                       gTATTCTGCAGGAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                                                                                                   <  2:725008/143‑1 (MQ=255)                        tATGCTTCAGCAGAAAAGGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCATCGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                  <  2:618020/143‑1 (MQ=255)                        gaTTCTGCAGCCGAAAAAGCGGCACCATTGCGTTTCTGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                  <  2:527298/142‑1 (MQ=255)                             ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGCAACTGACGAAATACCAAGCaaaag                                                                                                             >  1:138291/1‑141 (MQ=255)                             ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGACAACGGAACCGCGACACCAACAGGTAATGCAGAGCCTTCACCCAAACCACAGACAATGAAAATGAAGAAATAACAAGAAAATg                                                                                                             >  2:671276/1‑143 (MQ=255)                                 gcagAGCAAGCGGCGCGATGGCCTGTGAGCCAGGCGGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                         <  2:325776/143‑1 (MQ=255)                                  cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCCGAAGGCCACGGAACCGAGACAACAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAAGAAGAAATAACAAGCAAATGGcatc                                                                                                        >  1:183430/1‑143 (MQ=255)                                         ggCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCACATCCTTCTCCCAAACCAACGTTTATGAAAATGAATAAATAACAAGCAAATGGCATCATTCCTg                                                                                                 <  2:675717/142‑1 (MQ=255)                                                   ttGCATTTCATCCATCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATTAAAATGAATTAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                       <  1:166244/143‑1 (MQ=255)                                                   ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTATTACCAg                                                                                       >  2:186616/1‑143 (MQ=255)                                                               ccgccTGTTGGCGGAGTGGCTGCAGGCCACGGAACTGCTACACCATCAGTTATTGCAGAGCCTTCTCCCAAACCACCTTTCATGTAAATGAAAAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                           <  2:10375/141‑1 (MQ=255)                                                                       tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                                   <  1:682338/143‑1 (MQ=255)                                                                          cGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGATTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                >  2:660238/1‑143 (MQ=255)                                                                                 gCTGAAGGCCACGGAACCGGGACACGAACAGGTACTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                         <  2:239696/143‑1 (MQ=255)                                                                                  cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTGTACAAGGGAGAATTAACATGCTTAATGGCTATGTAcg                                                        >  1:515540/1‑143 (MQ=255)                                                                                      aGGCCACGGACCGGGGACAGCAACAGGTAATGCATAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa                                                    <  2:509868/143‑1 (MQ=255)                                                                                         ccACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATca                                                 >  2:586372/1‑143 (MQ=255)                                                                                         ccACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTCAACATGCTTATTGGCTATGTACGCGTATca                                                 >  2:236353/1‑143 (MQ=255)                                                                                           aCGGAAACGAGAAACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGCAAATAAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAGGCACATTGGCTATGAACGCGTATcaac                                               >  2:718530/1‑143 (MQ=255)                                                                                              gAACCGGGACACTAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAATAAGTTAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                            <  1:201616/143‑1 (MQ=255)                                                                                                   gggACACCAACAGGTAAGGGAGAGCTTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTATGCGTATCAACAAATGAcc                                       <  2:312919/143‑1 (MQ=255)                                                                                                      acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAGGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTAATAGACAATGTACGCGTATCAACAAATGACCAGa                                    >  1:29187/1‑143 (MQ=255)                                                                                                       cacCAACAGGTAATGCAGAGCCTTCTCGCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                   <  2:720751/143‑1 (MQ=255)                                                                                                          cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGATTATTGGGTACGTAAGAGTATCAGCAAAGGAGTAGAacac                                >  2:171769/1‑143 (MQ=255)                                                                                                            aCAGGTAATGCAGAGCATTCTACCAAAACAACGATTAGGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGATTTTACAAGGGGGATTTAACATGCATATTGGCTATGTAGGCGGATCAACAACTGACCCGAACACAg                              >  2:573352/1‑143 (MQ=255)                                                                                                                gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGCCCAGAACACAGATCt                          >  2:215872/1‑143 (MQ=255)                                                                                                                 tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAACTa                         >  2:372592/1‑143 (MQ=255)                                                                                                                    tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                      >  2:189086/1‑143 (MQ=255)                                                                                                                              tctcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcgc            <  1:284862/143‑1 (MQ=255)                                                                                                                                tcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg          <  1:586372/143‑1 (MQ=255)                                                                                                                                  ccAAACCAACGTTTATGACAATGAATAAATATCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCACCAAATGACCAGAACCCAGATCTACAACGTAATGCGCTGaa        <  2:102110/143‑1 (MQ=255)                                                                                                                                  ccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTCTTGGCTATGTACGAGTATCAACAAATGACCAGAAACCAGATCTACAACGTAATGCGATGaa        >  1:155269/1‑143 (MQ=255)                                                                                                                                  ccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGACATCATTACTGCTTTTACAAGAGGAATTTAACATGCTTATTGACTATGTACACGTAACAACAAATGACAAGAACACAGATCTACAACATAATGAGCAGaa        >  2:405087/1‑143 (MQ=255)                                                                                                                                   cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACGAGGGGGATTTACCATGCTTATTGGCTATGTAGNCGTATCAACAAATGACAAGAACACAGATCTCCAACGTAATgcgaggaag       >  1:377868/1‑137 (MQ=255)                                                                                                                                        cgaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc  <  1:516429/141‑1 (MQ=255)                                                                                                                                           |                                                                                                                                            GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGC  >  NC_000913/1209481‑1209758

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