New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2330390 =28 (0.770)6 (0.180) 6/242 3.6 16.0% pseudogene (4013/4605 nt) yfaS pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
?NC_000913 2330449 = 37 (1.100)pseudogene (3954/4605 nt) yfaS pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCTGCACATTTTGCGGCGATAATTCGTCACCAAAAGAGAGAATGTCCGGCACGCCCTGACCAACCCAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTATTAAGCAGCACAACCGTGTGCGCCAGCGGCTGA  >  NC_000913/2330307‑2330596
                                                                                                                                              |                                                                                                                                                   
ccTGCACATTTTGCGGCGATAATTCGTCACCAAAAGAGAGAATGTCCGGCACGCCCTGACCAACCCAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCGAAACAACTGGAGGCATTGTc                                                                                                                                               >  5:40476/1‑149 (MQ=255)
            gcggagATAATTCGTCACCAAAAGAGAGAATGTCCTGCACGCCCTGACCAACCCAACGCCAGTCTTCGCCCCCTCAAGTTAACTTATGTTTAGCCCATGCGCCCTCAGGCGCAGGCAAAACAACTTGAGGCATTGTCGCCATATAttt                                                                                                                                    <  3:23939/148‑1 (MQ=255)
                      ttCGTCACCAAAAGAGAGAATGTCCGGCACGCCCTGACCAACCCAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCat                                                                                                                                          >  4:56471/1‑132 (MQ=255)
                      ttCGTCACCAAAAGAGAGAATGTCCGGCACGCCCTGACCAACCCAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCat                                                                                                                                          <  3:56471/132‑1 (MQ=255)
                                      agaATGTCCGGCACGCCCTGACCAACCCAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAgcgcgcg                                                                                                         >  2:159142/1‑149 (MQ=255)
                                                   cGCCCTGACCAACCCAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGAt                                                                                            <  1:159142/149‑1 (MQ=255)
                                                       cTGACCAACCCAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCACCAGCCAGTTCATgg                                                                                                                   <  6:80657/122‑1 (MQ=255)
                                                       cTGACCAACCCAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCACCAGCCAGTTCATgg                                                                                                                   >  5:80657/1‑122 (MQ=255)
                                                          cccaacccaacGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCACGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAg                                                                                     >  2:123464/2‑149 (MQ=255)
                                                          aCCAACACAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAg                                                                                     >  2:22162/1‑149 (MQ=255)
                                                            caacccaacGCCAGTCTTCGCCCCCTCCAGTTAATTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGcc                                                                                                                                            >  1:72633/1‑92 (MQ=255)
                                                            caacccaacGCCAGTCTTCGCCCCCTCCAGTTAATTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGcc                                                                                                                                            <  2:72633/92‑1 (MQ=255)
                                                            caacccaacGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCGAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCg                                                                                   <  6:40476/149‑1 (MQ=255)
                                                            caacccaacGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCg                                                                                   >  8:104766/1‑149 (MQ=255)
                                                                  aacGCCAGTCTTCGCCACCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAAATGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCaaa                                                                             >  5:105318/1‑149 (MQ=255)
                                                                       cAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACt                                                                        >  5:3352/1‑149 (MQ=255)
                                                                                cccccTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGc                                                               <  7:104766/149‑1 (MQ=255)
                                                                                 ccccTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTg                                                                                          >  3:30474/1‑121 (MQ=255)
                                                                                 ccccTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTg                                                                                          <  4:30474/121‑1 (MQ=255)
                                                                                 ccccTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTc                                                                                >  7:127618/1‑131 (MQ=255)
                                                                                 ccccTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTc                                                                                <  8:127618/131‑1 (MQ=255)
                                                                                          ttAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATg                                                                                                                    <  4:138264/86‑1 (MQ=255)
                                                                                          ttAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATg                                                                                                                    >  3:138264/1‑86 (MQ=255)
                                                                                               ttATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGt                                                                                                <  6:162827/101‑1 (MQ=255)
                                                                                               ttATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGt                                                                                                >  5:162827/1‑101 (MQ=255)
                                                                                                        gCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCAcc                                       >  1:133840/1‑149 (MQ=255)
                                                                                                          ccATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGc                                     >  2:106655/1‑149 (MQ=255)
                                                                                                           cATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAAc                                                                            <  7:134728/109‑1 (MQ=255)
                                                                                                           cATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAAc                                                                            >  8:134728/1‑109 (MQ=255)
                                                                                                                gcCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTATTaa                               <  1:22162/149‑1 (MQ=255)
                                                                                                                  ccGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGGCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTATTAagc                             <  6:105318/149‑1 (MQ=255)
                                                                                                                   cgcaggcgcaggCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGcc                                                           >  6:151367/1‑118 (MQ=255)
                                                                                                                              gcAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTATTAagca                            >  2:60461/1‑138 (MQ=255)
                                                                                                                              gcAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTATTAagca                            <  1:60461/138‑1 (MQ=255)
                                                                                                                                aaaaCAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTATTAAGCAGCACAACCGtgtg               >  1:24956/1‑149 (MQ=255)
                                                                                                                                  aacaacTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTAtt                                 >  3:70591/1‑129 (MQ=255)
                                                                                                                                  aacaacTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTAtt                                 <  4:70591/129‑1 (MQ=255)
                                                                                                                                          gAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTATTAAGCAGCACAACCGTGTGCGCCAGCgg      <  2:133840/148‑1 (MQ=255)
                                                                                                                                             gCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTATTAAGCAGCACAACCGTGTGCGCCAGCGGCTGa  <  2:24956/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CCTGCACATTTTGCGGCGATAATTCGTCACCAAAAGAGAGAATGTCCGGCACGCCCTGACCAACCCAACGCCAGTCTTCGCCCCCTCCAGTTAACTTATGTTTAGCCCATGCGCCCGCAGGCGCAGGCAAAACAACTGGAGGCATTGTCGCCATATATTTCGCCAGCCAGTTCATGGCCAGCGCGCGCTCAATAGTGGATTGCTCAGCGGTCAAACCACTTAAAATAGCGGCCGTCTGGGTAGCGTCGCCACCGCTATTAAGCAGCACAACCGTGTGCGCCAGCGGCTGA  >  NC_000913/2330307‑2330596

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.