New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	= 4542860	17 (0.980)	4 (0.240)	4/266	1.9	17.5%	coding (486/549 nt)	fimA	major type 1 subunit fimbrin (pilin)
	?	NC_000913_3_pae_pgi	= 4542880	21 (1.240)					coding (506/549 nt)	fimA	major type 1 subunit fimbrin (pilin)
Rejected: Frequency below/above cutoff threshold.

AAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATT  >  NC_000913_3_pae_pgi/4542734‑4543007                                                                                                                               |                                                                                                                                                    aaaCGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAcccggggttgcggccccggttgc                                                                                                                              >  2:497‑M1/1‑127 (MQ=255)          tgtgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGc                                                                                                                     <  3:66201/150‑1 (MQ=255)            tgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCgg                                                                                                                   >  2:90941/1‑150 (MQ=255)             gCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGa                                                                                                                  <  2:168730/150‑1 (MQ=255)                                     tgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAcccggggttgcggccccggttgcaaaataacgcgcctggaacggaatggtattggttcc                                                                                          <  1:497‑M1/150‑60 (MQ=255)                                     tgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAACCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                                                                                          <  4:116380/150‑1 (MQ=255)                                                tGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc                                                                               >  1:97010/1‑150 (MQ=255)                                                  gATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                                                                             <  2:125528/150‑1 (MQ=255)                                                    tGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt                                                                           <  1:159789/150‑1 (MQ=255)                                                      gTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc                                                                         <  4:43427/150‑1 (MQ=255)                                                       tGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTCTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa                                                                        <  2:88179/150‑1 (MQ=255)                                                                    tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAcccggggttgcggccccggttgcaa                                                                                                                            <  1:29476‑M1/85‑26 (MQ=255)                                                                    tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAcccggggttgcggccccggttgcaa                                                                                                                            >  2:29476‑M1/1‑60 (MQ=255)                                                                         aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCa                                                      <  1:122399/150‑1 (MQ=255)                                                                                   aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCa                                            >  4:52019/1‑150 (MQ=255)                                                                                         ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGtgcaaaataacgcgcctggaacggaatggtattggttccgttattcagggttgtttctgaactaaatgtcgcaccatccagcgtcagcgca                                      <  1:497‑M2/150‑92 (MQ=255)                                                                                                    cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTACTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGATCAGGTACGTCATTACAGGCAGGGATGCCCACCCTTGTGCGATaaaaa                           >  4:19227/1‑150 (MQ=255)                                                                                                                  cgcgTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                                                                                       >  2:13320/1‑60 (MQ=255)                                                                                                                  cgcgTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                                                                                       <  1:13320/60‑1 (MQ=255)                                                                                                                    cgTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGaa           <  1:90941/150‑1 (MQ=255)                                                                                                                          tttGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                                                                                      <  3:29026/53‑1 (MQ=255)                                                                                                                          tttGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                                                                                      >  4:29026/1‑53 (MQ=255)                                                                                                                           ttGCAACCGGGGCCGCAACCCCGGGtgcaaaataacgcgcctggaacggaatggtattggttccgttattcagggttgtttctga                                                                     >  2:29476‑M2/1‑25 (MQ=255)                                                                                                                           ttGCAACCGGGGCCGCAACCCCGGGtgcaaaataacgcgcctggaacggaatggtattggttccgttattcagggttgtttctga                                                                     <  1:29476‑M2/85‑61 (MQ=255)                                                                                                                             gCAACCGGGGCCGCAACCCCGGGtgcaaaataacgcgcctggaacggaatggtattggttccgttattcagggttgtttctgaactaaatgtcgcaccatccagcgtcagcgcagcacccgttctgtccaggatctgcacaccaacgttt  >  2:497‑M2/1‑23 (MQ=255)                                                                                                                               |                                                                                                                                                    AAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATT  >  NC_000913_3_pae_pgi/4542734‑4543007

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