New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE = 80092815 (0.600)3 (0.130) 3/266 3.0 15.7% intergenic (+62/‑289) hisC_1/USA300HOU_RS03915 Histidinol‑phosphate aminotransferase/Putative 5'(3')‑deoxyribonucleotidase
?USA300TCH1516_ALE 801948 = 18 (0.760)intergenic (+189/+435) USA300HOU_RS03915/USA300HOU_RS03920 Putative 5'(3')‑deoxyribonucleotidase/Putative lipid kinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TAATTGTGACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTC  >  USA300TCH1516_ALE/801813‑802079
                                                                                                                                       |                                                                                                                                   
tAATTGTGACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAg                                                                                                                                <  1:15644/141‑1 (MQ=255)
gaaTTGTGGCTGAGATGAACTGTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAg                                                                                                                                <  2:198370/140‑1 (MQ=255)
  atggcgACTGAGATGAACGTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGaa                                                                                                                              <  1:37013/136‑1 (MQ=255)
  aTTGTGACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAACTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGaa                                                                                                                              <  1:246151/141‑1 (MQ=255)
             gaagAACTTTTCTGCCTTAGCCACTACAACAATATATTGGCAGTAGTTGACTGCGGGGCCCCAACATATAGAAATCGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCt                                                                                                                   <  1:102808/138‑1 (MQ=21)
                            cTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTAtt                                                                                                    >  1:31825/1‑141 (MQ=255)
                            cTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTAtt                                                                                                    >  2:236581/1‑141 (MQ=255)
                                         acacTATATTGGCACTAGTTGACTGCGGGGCGCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGAGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAg                                                                                       >  2:250979/1‑141 (MQ=255)
                                          cacTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAgt                                                                                      <  2:194790/141‑1 (MQ=255)
                                            cTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAgtgt                                                                                    >  2:225967/1‑141 (MQ=255)
                                                        tagtTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATCTCTTATGCATGAGTGTAACTCATGCAtt                                                                        >  2:270212/1‑141 (MQ=255)
                                                             gACGGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGTGGTGGGACCCAATATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCatat                                                                   <  2:215561/141‑1 (MQ=255)
                                                                        cccAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTaca                                                        <  1:247386/141‑1 (MQ=255)
                                                                                      actGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTATATCATAATGATACCTCACGCTCTTTTGTGACTTATTCATTCATATTTTTAAGTACACATTAGCTGTGACt                                          <  1:135789/139‑1 (MQ=9)
                                                                                         ggAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAAt                                       >  1:4057/1‑141 (MQ=255)
                                                                                                                    aGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCAAAATCTTCATTTCGCATTAGCTGTGACTAATGATAAAGAATCGCTACATAatcaatca            >  2:31931/1‑141 (MQ=18)
                                                                                                                       tGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCAtca         <  2:219599/141‑3 (MQ=255)
                                                                                                                              ggCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTACTCATACATAATCGCTACATAATAAATCATTAGTCgacc  >  2:201810/1‑138 (MQ=17)
                                                                                                                                       |                                                                                                                                   
TAATTGTGACTGAGATGAACTTTTATGTCTTAGACACTACAACACTATATTGGCAGTAGTTGACTGCGGGGCCCCAACATAGAGAAATTGGAACGCCAATTTCTACAGACAATGCAAGTTGGGGTGGGCCCCAACATAAAGAAATACTTTTTCTTTAGAAATTAGTATTTCTTATGCATGAGTGTAACTCATGCATTCATATTTTTAAGTACACATTAGCTGTGACTAATGATAAAGAATCGCTACATAATCAATCATTAGTCGTTC  >  USA300TCH1516_ALE/801813‑802079

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