New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 4040683107 (1.680)21 (0.340) 12/254 5.6 20.5% intergenic (+47/+223) rrfA/mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
?NC_000913_3_pae_tpiA = 4213078 58 (0.930)intergenic (+3/‑72) rrfE/yjaA 5S ribosomal RNA of rrnE operon/stress‑induced protein
Rejected: Coverage evenness skew score above cutoff.

ACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_tpiA/4040545‑4040683
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTG  <  NC_000913_3_pae_tpiA/4213078‑4212930
                                                                                                                                                                                                                                                                                                
ACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCAAAGACGGAAGACTGGGCCTTTCGTTTTATTTGATGCCTGG                                                                                                                                            <  1:367601/150‑1
                                                                             GGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGTGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGG                                                               >  2:344282/1‑150
                                                                                          ATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATG                                                                 <  4:123583/135‑1
                                                                                          ATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATG                                                                 >  3:123583/1‑135
                                                                                                    ACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAC                                        >  4:113116/1‑150
                                                                                                       AAAGGCTCAGTCGGGAGACTGGGCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGAC                                                   >  1:158954/1‑136
                                                                                                       AAAGGCTCAGTCGGGAGACTGGGCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGAC                                                   <  2:158954/136‑1
                                                                                                         AGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCG                                   >  1:375618/1‑150
                                                                                                         AGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACA                                                                                                           <  1:235514/78‑1
                                                                                                         AGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACA                                                                                                           >  2:235514/1‑78
                                                                                                           GCTCAGTCGGAAGACTGGGCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCC                                 >  3:456835/1‑150
                                                                                                                    GAAGACTGGGCCTTTCGTTTTATCTGATGCCTGGCAGTTCCCTACTCTCGCATGAGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATT                        <  3:113116/150‑1
                                                                                                                       GACTGGGCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTG                     <  1:336726/150‑1
                                                                                                                             GCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTAT               <  1:723622/150‑1
                                                                                                                             GCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTA                                                                                                                                  <  1:141764/35‑1
                                                                                                                             GCCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTA                                                                                                                                  >  2:141764/1‑35
                                                                                                                              CCTTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATC              <  1:498106/150‑1
                                                                                                                                TTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCT                                                                             >  1:681770/1‑85
                                                                                                                                TTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCT                                                                             <  2:681770/85‑1
                                                                                                                                TTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCACACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTG                     <  2:144381/141‑1
                                                                                                                                TTTCGTTTTATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCACACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTG                     >  1:144381/1‑141
                                                                                                                                          TTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTG  <  3:120157/150‑1
                                                                                                                                                                                                                                                                                                
ACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_tpiA/4040545‑4040683
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTG  <  NC_000913_3_pae_tpiA/4213078‑4212930

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