New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE = 196574366 (0.670)5 (0.060) 4/146 NT 7.4% coding (65/924 nt) hemH Ferrochelatase
?USA300TCH1516_ALE = 1965774 65 (0.720)coding (34/924 nt) hemH Ferrochelatase

GGTGATAAACCACCTATAAATTCATATCTATCTTTCAAATCTTGAAGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATATGGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/1965647‑1965743
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatggctTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGACATGGTAAACGTCCATCTGAAGAAGAACTTCAAG  <  USA300TCH1516_ALE/1965774‑1965687
                                                                                                                                                                                         
GGTGATAAACCACCTATAAATTCATATCTATCTTTCAAATCTTGAAGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAAT                                                                                                >  1:3153145/1‑91
  TGATAAACCACCTATAAATTCATATCTATCTTTCAAATCTTGAAGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATAT                                                                                              >  1:1434973/1‑91
   GATAAACCACCTATAAATTCATATCTATCTTTCAAATCTTGAAGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATATG                                                                                             <  1:362541/91‑1
   GATAAACCACCTATAAATTCATATCTATCTTTCAAATCTTGAAGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATATG                                                                                             >  1:2511040/1‑91
    ATAAACCACCTATAAATTCATATCTATCTTTCAAATCTTGAAGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATATGG                                                                                            >  1:61953/1‑91
                                            AAGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATA                                                    <  1:1524983/91‑1
                                             AGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATAT                                                   <  1:2379398/91‑1
                                             AGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATAT                                                   <  1:2378381/91‑1
                                                 CTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATA                                               <  1:337249/91‑1
                                                             GACGTTTACCATGTCTAATATCTGTATAATATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGAC                                   <  1:1630627/91‑1
                                                                                          TATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGACATGGTAAACGTCCATCTGAAGAAGAACTT      <  1:2824859/91‑1
                                                                                          TATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGACATGGTAAACGTCCATCTGAAGAAGAACTT      <  1:3771871/91‑1
                                                                                           ATGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGACATGGTAAACGTCCATCTGAAGAAGAACTTC     <  1:3989068/91‑1
                                                                                            TGGCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGACATGGTAAACGTCCATCTGAAGAAGAACTTCA    <  1:1557457/91‑1
                                                                                              GCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGACATGGTAAACGTCCATCTGAAGAAGAACTTCAAG  >  1:3027821/1‑91
                                                                                              GCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGACATGGTAAACGTCCATCTGAAGAAGAACTTCAAG  >  1:314953/1‑91
                                                                                              GCTTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGACATGGTAAACGTCCATCTGAAGAAGAACTTCAAG  >  1:316478/1‑91
                                                                                                                                                                                         
GGTGATAAACCACCTATAAATTCATATCTATCTTTCAAATCTTGAAGTTCTTCTTCAGATGGACGTTTACCATGTCTAATATCTGTATAATATGGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  USA300TCH1516_ALE/1965647‑1965743
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatggctTATGGCACACCTTATAAAGAAAGTGACATAGAGCCATATTATACAGATATTAGACATGGTAAACGTCCATCTGAAGAAGAACTTCAAG  <  USA300TCH1516_ALE/1965774‑1965687

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