New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 454286017 (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.

TGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTG  >  NC_000913_3_pae_pgi/4542743‑4543023
                                                                                                                                         |                                                                                                                                               
tgtgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGc                                                                                                                                     <  3:66201/150‑1 (MQ=255)
  tgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCgg                                                                                                                                   >  2:90941/1‑150 (MQ=255)
   gCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGa                                                                                                                                  <  2:168730/150‑1 (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)
                                                               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                                                                                                                       <  1:13320/60‑1 (MQ=255)
                                                                                                        cgcgTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                                                                                                       >  2:13320/1‑60 (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)
                                                                                                                       ccGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACCCCCTCACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTa              >  2:94085/1‑150 (MQ=255)
                                                                                                                         ggggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTAtt            >  4:118540/1‑150 (MQ=255)
                                                                                                                               gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGcgtc      >  2:11910/1‑150 (MQ=255)
                                                                                                                                   ccccGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTtg  >  4:34588/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                               
TGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTG  >  NC_000913_3_pae_pgi/4542743‑4543023

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.