New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4400624112 (0.870)14 (0.120) 7/166 NT 11.6% intergenic (+28/‑48) hfq/hflX global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication/GTPase, stimulated by 50S subunit binding; Mn(2+) homeostasis regulator
?NC_000913 = 4400640 109 (0.910)intergenic (+44/‑32) hfq/hflX global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication/GTPase, stimulated by 50S subunit binding; Mn(2+) homeostasis regulator
Rejected: Frequency below cutoff threshold.

ACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4400522‑4400624
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggggaACGCAGGATCGCTGGCTCCCCGTGTAAAAAAACAGCCCGAAACCTTATTCGGTTTCTTCGCTGTCCTGTTGCGCGGAAGTATTCTGCGCGCTGCTAC  <  NC_000913/4400640‑4400544
                                                                                                                                                                                                        
ACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGG                                                                                                      <  1:1518507/100‑1
  CAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGG                                                                                                    <  2:3593802/100‑1
   gcgcgctgCTACCATGATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGA                                                                                                   >  1:645483/9‑100
     CAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGCATAAGGTTTCGGGCTGTTTTTTTAAACAGG                                                                                                     >  2:214433/1‑96
                ATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAGCGCAGGATCGC                                                                                       <  2:222144/99‑1
                       TAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCC                                                                               <  1:640435/100‑1
                       TAGCAGCGCGCAGAATACTTCCGCCCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTGTTTTTACACGGGGAACGCAGTATCGCTGGCTCCC                                                                               <  2:1597550/100‑1
                          CAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCCCGT                                                                            <  2:717792/100‑1
                           AGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGATCGCAGGATCGCTGGCTCCCCGTG                                                                           <  2:1297299/100‑1
                           AGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCCCGTG                                                                           <  1:3023546/100‑1
                           AGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCCCGTG                                                                           <  1:1800502/100‑1
                            GCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCCCGTGT                                                                          <  1:3351542/100‑1
                            GCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCCCGTGT                                                                          <  1:447960/100‑1
                            GCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAAGAAGGTTCCGGGATGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCCCGTGT                                                                          <  2:202712/100‑1
                             CGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCCCGTGTA                                                                         <  2:3446925/100‑1
                                GCAGAATACTTCAGCGCAACAGGACAGCGAAGAAGCCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGAAACGCAGGATCGCTGGGTCCCCGTG                                                                           <  2:3284355/95‑1
                                  AGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCCCGTGTAAAAAA                                                                    <  1:1786940/100‑1
                                  AGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGAACGCAGGATCGCTGGCTCCCCGTGTAAAAAA                                                                    <  1:3608554/100‑1
                                                                                                  GGGGAACGCAGGATCGCTGGCTCCCCGTGTAAAAAAACAGCCCGAAACCTTATTCGGTTTCTTCGCTGTCCTGTTGCGCGGAAGTATTCTGCGCGCTGCT    >  2:3093480/1‑100
                                                                                                  GGGGAACGCAGGATCGCTGGCTCCCCGTGTAAAAAAACAGCCCGAAACCTTATTCGGTTTCTTCGCTGTCCTGTTGCGCGGAAGTATTCTGCGCGCTGCT    <  2:1342673/100‑1
                                                                                                  GGGGAACGCAGGATCGCTGGCTCCCCGTGTAAAAAAACAGCCCGAAACCTTATTCGGTTTCTTCGCTGTCCTGTTGCGCGGAAGTATTCTGCGCGCTGCT    >  1:2868744/1‑100
                                                                                                   GGGAACGCAGGATCGCTGGCTCCCCGTGTAAAAAAACAGCCCGAAACCTTATTCGGTTTCTTCGCTGTCCTGTTGCGCGGAAGTATTCTGCGCGCTGCTA   >  1:2918333/1‑100
                                                                                                    GGAACGCAGGATCGCTGGCTCCCCGTGTAAAAAAACAGCCCGAAACCTTATTCGGTTTCTTCGCTGTCCTGTTGCGCGGAAGTATTCTGCGCGCTGCTAC  >  2:2028522/1‑100
                                                                                                    GGAACGCAGGATCGCTGGCTCCCCGTGTAAAAAAACAGCCCGAAACCTTATTCGGTTTCTTCGCTGTCCTGTTGCGCGGAAGTATTCTGCGCGCTGCTAC  >  2:2820461/1‑100
                                                                                                                                                                                                        
ACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGGACAGCGAAGAAACCGAATAAGGTTTCGGGCTGTTTTTTTACACGGGGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4400522‑4400624
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggggaACGCAGGATCGCTGGCTCCCCGTGTAAAAAAACAGCCCGAAACCTTATTCGGTTTCTTCGCTGTCCTGTTGCGCGGAAGTATTCTGCGCGCTGCTAC  <  NC_000913/4400640‑4400544

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