Predicted mutation
evidence seq id position mutation freq annotation gene description
JC USA300TCH1516_ALE 2,373,573 Δ9 bp 16.3% coding (195‑203/1335 nt) pbuG ← Guanine/hypoxanthine permease PbuG

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE = 237357216 (0.680)3 (0.130) 3/256 NT 16.3% coding (204/1335 nt) pbuG Guanine/hypoxanthine permease PbuG
?USA300TCH1516_ALE 2373582 = 15 (0.650)coding (194/1335 nt) pbuG Guanine/hypoxanthine permease PbuG

CTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGC  >  USA300TCH1516_ALE/2373445‑2373717
                                                                                                                                         |                                                                                                                                       
cTGTTAAACCAACGTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCACACCTGCTACTAACGGGATTGAATATTTAGCTATGAGCCCCATGAATAGCGAGCCTACAAATGCTGCTa                                                                                                                                      <  1:64340/141‑1 (MQ=255)
 tgatAAACCAACGTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTACTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTaa                                                                                                                                     <  1:106663/138‑1 (MQ=255)
       aCCAACTTGCCAAGGAATACCCATGGTTGACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAAATAGATATTTAGCTATTAGCCCCATGAATAGCGAGCACACACATTATGCTAATGCAGt                                                                                                                               >  2:193506/1‑141 (MQ=255)
        ccAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTACCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCg                                                                                                                             <  2:244225/141‑1 (MQ=255)
                ccAAGGTATACCCCTGGTTAACACAACAGTTAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTACCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaaaa                                                                                                                     <  1:175491/141‑1 (MQ=255)
                         aCCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCAcc                                                                                                             >  1:274173/1‑141 (MQ=255)
                                          aggTAAATGAAAAGAATGCATTTAATCCCATACCTGGTGCGAACGAAATTAGATATTTATCTAATAGCCCGATGAATAGCGAGCCTACAAATTTTGCTAATAAAAAAAAAAAAAAAATGGCACCTTGGTCCATTTTCatat                                                                                            <  2:77295/139‑1 (MQ=255)
                                                                 aaTCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTCCATATCTTCCGATCACCCCTTAacacgc                                                                     >  2:46033/1‑139 (MQ=255)
                                                                   tgcatttaatcccatacctggtgctaacgcaattggatatttagctattagccccatgaatagcgagcctGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTCACACCTGc                                                                   >  2:77706‑M2/71‑141 (MQ=255)
                                                                              gTGCTAACGCAAATTGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTaaaa                                                        >  2:107899/1‑141 (MQ=255)
                                                                                            ggATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACt                                          <  1:241790/141‑1 (MQ=255)
                                                                                             gATATTTAGCTATTACCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGc                                        >  2:101175/1‑141 (MQ=255)
                                                                                                caatttgatatttagctattagccccatgaatagcgagcctGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTa                                      <  2:265666‑M2/100‑1 (MQ=255)
                                                                                                atTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTa                                      >  2:53618/1‑141 (MQ=255)
                                                                                                 ttgaGCAATTAGACACATGAATAGCGAGTCTACAAATGCTGATAATGCAGTCGCTACAAAAATGGCACCTTGGTGCATTTCCATATCTTCTGATACGTCTCTATCGAGAAATAGTCTTCAAACGTGCGGGTTAACTGCTaa                                     <  2:134999/138‑1 (MQ=255)
                                                                                                       tactAGCCCCATGAAGAGCGAGCGTACAAATGCTGCTAATGCAGTAGTCACAAAAAGGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAAtata                               <  2:285747/138‑1 (MQ=255)
                                                                                                              cccATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCCCAGGCGATGGCTCGCCTTTCCCGGCTGCGAAACTGAGAACGTGCGGGTTAACTGCTAAAATATAGGCCATa                        <  1:216930/141‑1 (MQ=255)
                                                                                                                        ccatgaatatcgagcctGCTAATTCAGGCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAg              <  2:219913‑M2/124‑1 (MQ=255)
                                                                                                                         cGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGt             <  1:213127/141‑1 (MQ=255)
                                                                                                                               tACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTATTTAGTGCGAAACACAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGAt       >  1:138705/1‑141 (MQ=255)
                                                                                                                                    aTGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTCCGGGTTAACTGCTACAATATAGGCCCTCGATACAAAACTTCTCATAccac  >  1:274047/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGC  >  USA300TCH1516_ALE/2373445‑2373717

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