New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3540274 =26 (0.890)3 (0.110) 3/256 NT 11.1% coding (112/228 nt) feoA ferrous iron transport protein A
?NC_000913 3540310 = 23 (0.810)coding (148/228 nt) feoA ferrous iron transport protein A

ATATGAATGGGGTCGCCGAGTGGAGCGACGCGCACCACATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3540314‑3540274
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cattCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATCCAAATTCTGG  >  NC_000913/3540310‑3540450
                                                                                                                                                                                      
ATATGAATGGGGTCGCCGAGTGGAGCGACGCGCACCACATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATCAA                                         >  1:560906/1‑143
 TATGAATGGGGTCGCCGAGTGGAGCGACGCGCACCACATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTACTACGTTGATAACAACAATGAAA                                        >  2:91203/1‑143
  ATGAATGGGGTCGCCGAGTGGAGCGACGCGCACCACATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATATGGCCGCATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAA                                       >  2:346377/1‑143
                                      ATTCATGTCGTACCCTGTCGTGTGCGGCTGGTATTCCGCCAAAAATATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTACTACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATCCAAATTCTG   <  2:85828/143‑1
                                       TTCGTATCTGAACCCGTCGTGTGAGGCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAAGACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATCCAAATTCTGG  <  2:514463/143‑1
                                                                                                                                                                                      
ATATGAATGGGGTCGCCGAGTGGAGCGACGCGCACCACATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3540314‑3540274
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cattCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATCCAAATTCTGG  >  NC_000913/3540310‑3540450

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