New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 454343836 (1.040)4 (0.120) 4/252 3.9 12.4% coding (324/549 nt) fimA major type 1 subunit fimbrin (pilin)
?NC_000913 = 4543474 23 (0.710)coding (360/549 nt) fimA major type 1 subunit fimbrin (pilin)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGTTCTGCTGTCGGTTTTAACATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGG  >  NC_000913/4543328‑4543619
                                                                                                                                                  |                                                                                                                                                 
aGTTCTGCTGTCGGTTTTAACATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACaaa                                                                                                                                                 <  8:84690/149‑1 (MQ=255)
     tgctgTCGGTTTTAACATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTg                                                                                                                                            >  2:64698/1‑149 (MQ=255)
       ctgTCGGTTTTAACATTCAGCTGAATGATTGCGATCCCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTgg                                                                                                                                           <  4:113871/148‑1 (MQ=255)
       ctgTCGGTTTTAACATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTgg                                                                                                                                           <  5:159609/148‑1 (MQ=255)
       ctgTCGGTTTTAACATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTg                                                                                                                                            <  2:53080/147‑1 (MQ=255)
         gTCGGTTTTAACATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGgtgt                                                                                                                                        >  3:84796/1‑149 (MQ=255)
          tCGGTTTTAACATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGgtgt                                                                                                                                        <  5:133968/148‑1 (MQ=255)
            ggTTTTAACATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGtgtg                                                                                                                                       >  2:171867/1‑147 (MQ=255)
              ttttAACATTCAGCTGAATGATTGCGACACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGa                                                                                                                                   <  8:108843/149‑1 (MQ=255)
                  aaCATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCt                                                                                                                               >  8:102996/1‑149 (MQ=255)
                    cATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTgg                                                                                                                             <  7:102996/149‑1 (MQ=255)
                     aTTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGc                                                                                                                                      >  6:49171/1‑139 (MQ=255)
                     aTTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGc                                                                                                                                      <  5:49171/139‑1 (MQ=255)
                     aTTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGa                                                                                                                            <  4:84796/149‑1 (MQ=255)
                                                     aTCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCTGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTg                                                                                            >  4:84755/1‑149 (MQ=255)
                                                          aaGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTg                                                                                                                                            <  7:88065/96‑1 (MQ=255)
                                                          aaGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTg                                                                                                                                            >  8:88065/1‑96 (MQ=255)
                                                                     gCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACAGTTAGTTCAg                                                                             >  8:50570/1‑148 (MQ=255)
                                                                       cTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGaaa                                                                          <  1:64698/149‑1 (MQ=255)
                                                                                     cGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAaacaa                                                                       >  7:37859/1‑138 (MQ=255)
                                                                                     cGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAaacaa                                                                       <  8:37859/138‑1 (MQ=255)
                                                                                                             ttCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGaa                                                        <  4:142235/129‑1 (MQ=255)
                                                                                                             ttCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGaa                                                        >  3:142235/1‑129 (MQ=255)
                                                                                                                    tctGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGt                              >  5:62365/1‑148 (MQ=255)
                                                                                                                       gCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTAtt                          >  6:52758/1‑149 (MQ=255)
                                                                                                                         agagTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGAGCCTGGACAGCACGGGGGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCCGAAACAACCCTGAATAACGGAACCAATACCAATCCCTTCCAGGCGCGGTAtttt                        <  3:84755/149‑1 (MQ=255)
                                                                                                                          gagTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACACTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTg                       >  1:17023/1‑149 (MQ=255)
                                                                                                                                 gCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCgg                <  5:52758/149‑1 (MQ=255)
                                                                                                                                 gCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCgg                <  6:62365/149‑1 (MQ=255)
                                                                                                                                        gTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGc          >  5:102020/1‑148 (MQ=255)
                                                                                                                                            cgcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc      >  6:102020‑M2/1‑146 (MQ=255)
                                                                                                                                               aacaaACGTGGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata  >  1:116326‑M2/1‑143 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGTTCTGCTGTCGGTTTTAACATTCAGCTGAATGATTGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGG  >  NC_000913/4543328‑4543619

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.