New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4442010 =80 (1.050)3 (0.040) 3/264 NT 3.8% coding (167/639 nt) msrA methionine sulfoxide reductase A
?NC_000913 4442119 = 75 (1.040)coding (58/639 nt) msrA methionine sulfoxide reductase A
Rejected: Frequency below/above cutoff threshold.

TGATAAAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGGGGTGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4442165‑4442010
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggggtgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAA  >  NC_000913/4442119‑4442265
                                                                                                                                                                                                                                                                                                               
TGATAAAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTG                                                                                                                                                           >  1:919203/1‑150
TGATAAAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTG                                                                                                                                                           <  2:740521/150‑1
 GATAAAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGG                                                                                                                                                          <  1:409979/150‑1
  ATAAAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGGG                                                                                                                                                         <  1:616685/150‑1
  ATAAAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAAAGACCAATGTAACTGACGGAATGGAGATTGCCATGTTTGCGATGGGTTGTTTCAGGG                                                                                                                                                         >  2:1193079/1‑150
   TAAAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGGGG                                                                                                                                                        <  1:900309/150‑1
    AAAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCTTGCGGTCAACGGTCACTCAATGCCCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTGGTTTCTGGGGg                                                                                                                                                       >  2:74687/1‑149
     AAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGGGGTG                                                                                                                                                      >  1:541719/1‑150
     AAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGGGGTG                                                                                                                                                      >  2:1049258/1‑150
                                         GACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGA                                                                                                                  >  2:784570/1‑150
                                             TAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCT                                                                                                              >  1:694888/1‑150
                                                                                                                   ATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGG                                        <  1:784570/150‑1
                                                                                                                                                     GGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTC      <  1:730149/150‑1
                                                                                                                                                     GGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAAAACCACCCTGGAAATGTCT     >  1:600934/1‑150
                                                                                                                                                       GGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGGGTATTTGCACACATTAACCCCCGTGAATAGGTCTA    >  1:546040/1‑150
                                                                                                                                                         TGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAA  <  2:1014456/150‑1
                                                                                                                                                                                                                                                                                                               
TGATAAAAAGCATCTGGTTTCCCCCGCCGATGCCCTGCCTGGACGTAACACCCCGATGCCCGTAGCCACGCTGCATGCGGTCAACGGTCACTCAATGACCAATGTACCTGACGGAATGGAGATTGCCATTTTTGCGATGGGTTGTTTCTGGGGTGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4442165‑4442010
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggggtgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAA  >  NC_000913/4442119‑4442265

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