New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 335560636 (1.040)5 (0.150) 5/280 2.4 12.5% coding (882/4461 nt) gltB glutamate synthase, large subunit
?NC_000913 = 3355642 36 (1.100)coding (918/4461 nt) gltB glutamate synthase, large subunit
Rejected: Frequency below cutoff threshold.

ACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3355449‑3355606
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caggcgggTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTTATCCATCGAACTGGAGTCAGAGCCGGTTTCGTTGACGAACGGTGCGGCGTCGTGCAGGTCAGGGATAAGCGGTGTCTGGAATT  <  NC_000913/3355642‑3355494
                                                                                                                                                                                                                                                                                                                   
accCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGC                                                                                                                                                              <  2:19664/148‑1
      TCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGG                                                                                                                                                        >  2:154919/1‑151
                                                                         GCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTT                                                                                     >  2:818/1‑151
                                                                                                         GCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGCCAGGCGGGTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTTATCCATCGAACTGGAGTCAGAGCCGGTTTCGT                                                     >  2:186020/1‑151
                                                                                                                CTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTTATCCATCGAACTGGAGTCAGAGCCGGTTTCGTTGACGAA                                              >  2:315787/1‑151
                                                                                                                            GGATAACATGCTGGAACTGCTGCTGCCAGGCGGGTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTTATCCATCGAACTGGAGTCAGAGCCGGTTTCGTTGACGAACGGTGCGGCGT                                   >  2:216953/1‑150
                                                                                                                                 ACATGCTGGAACTGCTGCTGGCAGGCGGGTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTTATCCATCGAACTGGAGTCAGAGCCGGTTTCGTTGACGAACGGTGCGGCGTCGTGCA                             >  2:321414/1‑151
                                                                                                                                                          CGGGTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTTATCCATCGAACTGGAGTCAGAGCCGGTTTCGTTGACGAACGGTGCGGCGTCGTGCAGGTCAGGGATAAGCGGTGTCTGGAA    <  1:113647/151‑1
                                                                                                                                                          CGGGTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTTATCCATCGAACTGGAGTCAGAGCCGGTTTCGTTGACGAACGGTGCGGCGTCGTGCAGGTCAGGGATAAGCGGTGTCTGGAA    <  1:548610/151‑1
                                                                                                                                                            GGTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTTATCCATCGAACTGGAGTCAGAGCCGGTTTCGTTGACGAACGGTGCGGCGTCGTGCAGGTCAGGGATAAGCGGTGTCTGGAATT  >  2:188719/1‑151
                                                                                                                                                                                                                                                                                                                   
ACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3355449‑3355606
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caggcgggTGGTACTAATAGACGCATGGCGCGGATGATATCCATCCCGCCTGCCAGCAGCAGTTCCAGCATGTTATCCATCGAACTGGAGTCAGAGCCGGTTTCGTTGACGAACGGTGCGGCGTCGTGCAGGTCAGGGATAAGCGGTGTCTGGAATT  <  NC_000913/3355642‑3355494

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