New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4195804 =63 (0.970)7 (0.110) 6/246 NT 10.6% coding (296/1896 nt) thiC phosphomethylpyrimidine synthase
?NC_000913 4195826 = 59 (0.960)coding (274/1896 nt) thiC phosphomethylpyrimidine synthase

AGGGCTGGCAAAACTACGCCAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4195827‑4195804
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagccCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTTTCTTCGTACTGCGGCTGTTCTTTGCTACCGCCAATTAGCGTCGGGCTAAGCTGGATCTCACGCATCGGCAC  >  NC_000913/4195826‑4195970
                                                                                                                                                                         
AGGGCTGGCAAAACTACGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTT                                                                         <  1:62883/98‑1
AGGGCTGGCAAAACTACGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTT                                                                         >  2:62883/1‑98
  GGCTGGCAAAACTACGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTTTCTTCGTACTGCGGCTGTTCTTTGCTACCGCCAATTAGCGTCGGGCTAAGCTG                    >  4:35510/1‑149
   GCTGGCAAAACTACGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGA                                                                                                     >  1:64043/1‑67
   GCTGGCAAAACTACGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGA                                                                                                     <  2:64043/67‑1
   GCTGGCAAAACTACGCCAGCCCTTGCTGCACGTTAATAGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTTTCTTCGTACTGCGGCTGTTCTTTGCTACCGCCAATTAGCGTCGGGCTAAGCTGG                   >  3:551540/1‑149
       GCAAAACTACGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTTTCTTCGTACTGCGGCTGTTCTTTGCTACCGCCAATTAGCGTCGGGCTAAGCTGGATCT               >  1:304739/1‑149
                    AGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTTTCTTCGTACTGCGGCTGTTCTTTGCTACCGCCAATTAGCGTCGGGCTAAGCTGGATCTCACGCATCGGCAC  >  1:419836/1‑149
                                                                                                                                                                         
AGGGCTGGCAAAACTACGCCAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4195827‑4195804
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccagccCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTTTCTTCGTACTGCGGCTGTTCTTTGCTACCGCCAATTAGCGTCGGGCTAAGCTGGATCTCACGCATCGGCAC  >  NC_000913/4195826‑4195970

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