Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 3,801,411 IS5 (–) +4 bp coding (586‑589/1017 nt) waaR ← lipopolysaccharide 1,2‑glucosyltransferase; UDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2288919 =NA (NA)26 (0.700) 13/264 1.4 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 3801414 0 (0.000)coding (586/1017 nt) waaR lipopolysaccharide 1,2‑glucosyltransferase; UDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase
* ? NC_000913 = 2290113NA (NA)38 (1.030) 26/264 0.5 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 3801411 = 0 (0.000)coding (589/1017 nt) waaR lipopolysaccharide 1,2‑glucosyltransferase; UDP‑glucose:(glucosyl)LPS alpha‑1,2‑glucosyltransferase

GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGA  >  NC_000913/2289969‑2290255
                                                                                                                                                |                                                                                                                                              
gACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATaa                                                                                                                                            >  3:46882/1‑149 (MQ=32)
       aaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTa                                                                                                                                     <  4:46882/149‑1 (MQ=38)
                                          aaTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAagag                                                                                                   >  4:27326/1‑148 (MQ=255)
                                            ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaa                                                                                                 <  4:152142/148‑1 (MQ=255)
                                            ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaa                                                                                                 >  3:152142/1‑148 (MQ=255)
                                                 tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGaa                                                                                           >  1:177046/1‑149 (MQ=255)
                                                 tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGaa                                                                                           >  1:179844/1‑149 (MQ=255)
                                                       cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATCATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGAt                                                                                     >  8:59901/1‑149 (MQ=255)
                                                        gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGAc                                                                                                         <  1:102159/128‑1 (MQ=255)
                                                        gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGAc                                                                                                         >  2:102159/1‑128 (MQ=255)
                                                                 ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTa                                                                                                                                     <  2:178409/91‑1 (MQ=255)
                                                                 ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTa                                                                                                                                     >  1:178409/1‑91 (MQ=255)
                                                                           atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATc                                                                 >  2:157623/1‑149 (MQ=255)
                                                                              atgaACCCTGTTCTATGGCTCCAGATGACAAACATTATCTCATATCAGGGACTTGTTCTCACCTTCCATAACTCTGTATCCACAAGAACAGATATTGCGGAACGACAAAGAGAAACATAAACAGATTGATGCAATAAGCTTACATCCTa                                                              >  4:77953/1‑149 (MQ=25)
                                                                                        ttCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGaa                                                                                           <  6:141988/110‑1 (MQ=255)
                                                                                        ttCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGaa                                                                                           >  5:141988/1‑110 (MQ=255)
                                                                                             tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGa                                                                          >  2:175702/1‑122 (MQ=255)
                                                                                             tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGa                                                                          <  1:175702/122‑1 (MQ=255)
                                                                                             tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCt                                                <  2:177046/148‑1 (MQ=255)
                                                                                             tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCt                                                <  2:179844/148‑1 (MQ=255)
                                                                                                     gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGc                                       <  1:157623/149‑1 (MQ=255)
                                                                                                        gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAtt                                    >  4:152387/1‑149 (MQ=255)
                                                                                                                       atatCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGTGTACAGaaaa                        <  7:59901/146‑1 (MQ=255)
                                                                                                                          tCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGt                                              <  4:24016/121‑1 (MQ=255)
                                                                                                                          tCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGt                                              >  3:24016/1‑121 (MQ=255)
                                                                                                                             gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAg               <  3:152387/149‑1 (MQ=255)
                                                                                                                              ggACTTGTTCGCTCCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGATGTTAAAAGAATTGGGTACAGAAAATAc                     >  7:35526/1‑142 (MQ=255)
                                                                                                                              ggACTTGTTCGCACCTTCCATAACGCTTTAGCCACCAGAACAGATATTGAGGAACGACAAAGAGAAACAGAACCAGATTGATGCATGGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATAc                     <  8:35526/142‑1 (MQ=255)
                                                                                                                                       cGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATAcac     >  7:1477/1‑147 (MQ=255)
                                                                                                                                          acattCCATCACGCTTTAGCCACCAGACCAGATATTGCGGAACTACAAAGAGAAACAGAAACAGATTTATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGa  <  3:77953/146‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGA  >  NC_000913/2289969‑2290255

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