New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 233137033 (1.120)4 (0.140) 4/246 2.9 10.9% pseudogene (3033/4605 nt) yfaS pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
?NC_000913 = 2331375 34 (1.220)pseudogene (3028/4605 nt) yfaS pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Rejected: Frequency below/above cutoff threshold.

GTTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTG  >  NC_000913/2331254‑2331515
                                                                                                                    |                                                                                                                                                 
gTTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAActgtcttttgtggataacagctggcccgttg                                                                                                                    >  2:85441‑M1/1‑117 (MQ=255)
 tttCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttg                                                                                                                                         <  4:124095‑M2/126‑6 (MQ=255)
 tttCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttg                                                                                                                                         >  3:124095‑M2/1‑121 (MQ=255)
 tttCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGACCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGc                                                                                                                  >  4:79581/1‑149 (MQ=255)
       tactTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCAttt                                                                                                            >  5:65435/1‑149 (MQ=255)
               aGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAActgtcttttgtggataacagctggcccgttgaacagcagaaaaatgt                                                                                                    <  1:85441‑M1/149‑48 (MQ=255)
                   ggATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCt                                                                                                >  5:17511/1‑149 (MQ=255)
                    gATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgtg                                                                                                                                       <  4:28313‑M2/109‑8 (MQ=255)
                    gATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgtg                                                                                                                                       >  3:28313‑M2/1‑102 (MQ=255)
                        ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCAc                                                                                            <  6:17511/148‑1 (MQ=255)
                        ttGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGACCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCAc                                                                                            <  3:79581/148‑1 (MQ=255)
                          gCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACtt                                                                                          >  8:9486/1‑148 (MQ=255)
                          gCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACtt                                                                                          >  5:22276/1‑148 (MQ=255)
                            tCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTaa                                                                                        >  5:131783/1‑148 (MQ=255)
                                cTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCAc                                                                                                                      >  1:163139/1‑114 (MQ=255)
                                cTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCAc                                                                                                                      <  2:163139/114‑1 (MQ=255)
                                         cAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCt                                                                          <  7:9486/149‑1 (MQ=255)
                                           aCATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTg                                                                         >  2:68486/1‑148 (MQ=255)
                                              tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgt                                                                                                                                        <  3:12188‑M2/82‑7 (MQ=255)
                                              tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgt                                                                                                                                        >  4:12188‑M2/1‑76 (MQ=255)
                                              tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTttgttg                                                                                                       <  6:16522/115‑1 (MQ=255)
                                              tCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTttgttg                                                                                                       >  5:16522/1‑115 (MQ=255)
                                                                caccGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAActgtcttttgtggataac                                                                                                                                 >  7:116965‑M1/1‑53 (MQ=255)
                                                                caccGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAActgtcttttgtggataac                                                                                                                                 <  8:116965‑M1/71‑19 (MQ=255)
                                                                        aTGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGc                                           <  6:65435/149‑1 (MQ=255)
                                                                           aCATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctgatgttgcccgagcaggcgagcaatatccggct                                        <  1:85441‑M2/149‑103 (MQ=255)
                                                                             aTTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCtt                                      >  6:68935/1‑149 (MQ=255)
                                                                              tttttCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCttt                                     <  5:68935/149‑1 (MQ=255)
                                                                               ttttCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTtg                                    >  5:44935/1‑149 (MQ=255)
                                                                                tttCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTtg                                    <  6:22276/148‑1 (MQ=255)
                                                                                 ttCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTtgt                                   <  1:68486/148‑1 (MQ=255)
                                                                                      ctgTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagc                              <  6:118617/148‑1 (MQ=255)
                                                                                         ttCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtca                          <  6:44935/149‑1 (MQ=255)
                                                                                         ttCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgt                            <  6:131783/147‑1 (MQ=255)
                                                                                           cAACGGGCCAGCTGTTATCCACAAAAGACAGttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctgatgttgcccgagcaggcgagcaatatccggctacaaagtagtgaaac                         >  2:85441‑M2/1‑31 (MQ=255)
                                                                                            aaCGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtcagcg                       >  5:54013/1‑149 (MQ=255)
                                                                                            aaCGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtcagc                        >  1:35004/1‑148 (MQ=255)
                                                                                            aaCGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGcagcgtcagc                        <  2:35004/148‑1 (MQ=255)
                                                                                                     gCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGaaa                                                     <  4:124207/110‑1 (MQ=255)
                                                                                                       tGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCAtc            >  5:13266/1‑149 (MQ=255)
                                                                                                        gTTATCCACAAAAGACAGttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtg                                                                                         <  8:116965‑M2/71‑54 (MQ=255)
                                                                                                        gTTATCCACAAAAGACAGttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtg                                                                                         >  7:116965‑M2/1‑18 (MQ=255)
                                                                                                          tATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCAtct           >  1:113/1‑147 (MQ=255)
                                                                                                             ccACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAAtt                                              <  5:51713/109‑1 (MQ=255)
                                                                                                             ccACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAAtt                                              >  6:51713/1‑109 (MQ=255)
                                                                                                              cacaAAActgtcttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctgatgttgcccgagcaggcgagcaatatc                                             <  4:28313‑M1/109‑103 (MQ=255)
                                                                                                              cacaAAActgtcttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctgatgttgcccgagcaggcgagcaatatc                                             >  3:28313‑M1/1‑7 (MQ=255)
                                                                                                               acaAAActgtcttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctga                                                                       >  4:12188‑M1/1‑6 (MQ=255)
                                                                                                               acaAAActgtcttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctga                                                                       <  3:12188‑M1/82‑77 (MQ=255)
                                                                                                                caAAActgtcttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctgatgttgcccgagcaggcgagcaatatccggctacaaagtagtgaaa                          <  4:124095‑M1/126‑122 (MQ=255)
                                                                                                                caAAActgtcttttgtggataacagctggcccgttgaacagcagaaaaatgtcatgctcggtggtggcgataacgcgctgatgttgcccgagcaggcgagcaatatccggctacaaagtagtgaaa                          >  3:124095‑M1/1‑5 (MQ=255)
                                                                                                                caAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCt                    >  2:31991/1‑132 (MQ=255)
                                                                                                                 aaaaGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTg  >  3:128880/1‑149 (MQ=255)
                                                                                                                    |                                                                                                                                                 
GTTTCACTACTTTGTAGCCGGATATTGCTCGCCTGCTCGGGCAACATCAGCGCGTTATCGCCACCACCGAGCATGACATTTTTCTGCTGTTCAACGGGCCAGCTGTTATCCACAAAAGACAGTTTTGTGCTAATGCTGTCCTGCACTTGCCCATTTTGTTGCAGTTCTGCACTTAACAAGCCAGATTGCTGAATATTCTGCGTCAGCGAAATATAATTCGCCCCTTTGTGCAGCGTCAGCGTCTGGCGCATCTCAGCGCCTG  >  NC_000913/2331254‑2331515

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