Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 3,968,900 IS5 (–) +4 bp 100% coding (985‑988/1104 nt) wecA → UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2288919 =NA (NA)28 (1.530) 25/380 NT 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 3968903 0 (0.000)coding (988/1104 nt) wecA UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
* ? NC_000913 = 3653230NA (NA)29 (1.580) 22/380 NT 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 3968900 = 0 (0.000)coding (985/1104 nt) wecA UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

GCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATC  >  NC_000913/2288715‑2289134
                                                                                                                                                                                                            |                                                                                                                                                                                                                        
gcCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                                                                                                                                                                                         <  2:86821/238‑1 (MQ=255)
            aaTGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATCTTTTTTCAGGACAAAACCATCAACTGTCTCACCATACTCAACAAGTGATGCCAATGCTATGTGTGAATATaagg                                                                                                                                                                                                                         >  2:180575/1‑194 (MQ=255)
                                       aGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATCTTTTTTCAGGACAACACCATCAACTGTCTCCCCATACTCCACAGTTGATGCCAATGCTATGGGGGAATATaaggaag                                                                                                                                                                                                                      >  2:124282/1‑170 (MQ=255)
                                                   tacACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTTCCATAAACCAGCTGGGTATCTTTTTTCAGGACAACACCATCAACTTTCTCACCATAATCAACAGTTGATGCCAATGCTATGGGGGAATATAAGGAAGGTGCGAATAAGCGGGGGAATTCTTCTCAGCTGACTCCGTCATTTCATTTCTTCATGTTTTAGCCGATTTTTTTTCCCGTaa                                                                                                                                     >  2:153648/1‑239 (MQ=255)
                                                                  tataTGCTGTTCCCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATACACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                                                                                                     <  2:149220/239‑1 (MQ=255)
                                                                                                         ggCGGTCCCCTAAACCAGCTGGATATCTTTTTTCAGGACAACACCCTCAACTGTCTCACCATACTCCACAGTTGATGCCCATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGGAATTCTTTTCGGCTGAATTAGTCATTTTATTTTTTCCTGTTTGAGCCGCTTTTTTTTCCCGTAAATGCCTTTAATCAGGCCTATTTGAACGCTTCTTTGCCCTTTatagggctgttc                                                                               >  2:107590/1‑227 (MQ=255)
                                                                                                                                   cTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACt                                                    <  1:153648/239‑1 (MQ=255)
                                                                                                                                                                        cTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCa                >  1:54845/1‑238 (MQ=255)
                                                                                                                                                                                  gATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAAc     >  2:212379/1‑239 (MQ=40)
                                                                                                                                                                                     gCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGCAACAGGTTGGCCAGCGTGAATAACATc  >  1:84843/1‑239 (MQ=39)
                                                                                                                                                                                     gCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACa    >  1:136373/1‑237 (MQ=39)
                                                                                                                                                                                     gCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACa    <  2:136373/237‑1 (MQ=39)
                                                                                                                                                                                                            |                                                                                                                                                                                                                        
GCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATC  >  NC_000913/2288715‑2289134

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