Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 2,866,695 IS1 (+) +9 bp 50.0% coding (849‑857/993 nt) rpoS ← RNA polymerase, sigma S (sigma 38) factor

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1050545NA (NA)4 (0.300) 4/262 NT 42.9% noncoding (768/768 nt) IS1 repeat region
?NC_000913 2866695 = 5 (0.400)coding (857/993 nt) rpoS RNA polymerase, sigma S (sigma 38) factor
* ? NC_000913 = 28667035 (0.400)6 (0.450) 5/264 NT 52.8% coding (849/993 nt) rpoS RNA polymerase, sigma S (sigma 38) factor
?NC_000913 3583428 = NA (NA)noncoding (1/768 nt) IS1 repeat region

ACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCT  >  NC_000913/2866534‑2866886
                                                                                                                                                                         |                                                                                                                                                                                       
aCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATGTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTgg                                                                                  <  1:242969/273‑1 (MQ=255)
                                              gATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGggtgatgctgccaacttactgatttagtgtatgatggtgttttt                                                                                                                                             >  2:155408‑M1/1‑124 (MQ=255)
                                                 atTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgag                                                                                                                                          <  1:155409‑M1/168‑48 (MQ=255)
                                                   cacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGCCGCTTCGTACCCCAGCAAACCGAATCGACGt                                                                                                                                                                 >  1:240703‑M2/111‑143 (MQ=255)
                                                      ggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGcc                                                                                                                                                              <  2:240702‑M2/36‑1 (MQ=255)
                                                                                                                              acataaaacactatcaataagttggagtcattaccGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCAt                                                                                                                       >  2:100170‑M2/36‑110 (MQ=255)
                                                                                                                                 taaaacactatcaataagttggagtcattaccGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTg                                                                                                                    <  1:100171‑M2/78‑1 (MQ=255)
                                                                                                                                 ggCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACtt                                                                                                                                                       >  2:217602/1‑75 (MQ=255)
                                                                                                                                    cAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCAc                                                                                                                                                    <  1:217603/75‑1 (MQ=255)
                                                                                                                                             gACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGa                                                                                                                >  2:66631/1‑102 (MQ=255)
                                                                                                                                                cTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGc                                                                                                             <  1:66632/102‑1 (MQ=255)
                                                                                                                                                          cAGTGTTGCCGCTTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagc                                                             >  2:40518‑M1/1‑16 (MQ=255)
                                                                                                                                                          cAGTGTTGCCGCTTCGggtgatgctgccaacttactgatttagtgtatgatggtgtt                                                                                                                                                >  2:205677‑M1/1‑16 (MQ=255)
                                                                                                                                                             tgtTGCCGCTTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcacc                                                          <  1:40518‑M1/140‑128 (MQ=255)
                                                                                                                                                             tgtTGCCGCTTCGggtgatgctgccaacttactgatttagtgtatgatggtgttttt                                                                                                                                             <  1:205678‑M1/57‑45 (MQ=255)
                                                                                                                                                              gtTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCg                                                                                     >  1:266173/1‑112 (MQ=255)
                                                                                                                                                              gtTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCtttt                             >  2:169538/1‑168 (MQ=255)
                                                                                                                                                                 GCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTgg                                                                                  <  2:266172‑M2/112‑1 (MQ=255)
                                                                                                                                                                 GCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCg                          <  1:169539‑M2/168‑1 (MQ=255)
                                                                                                                                                                  ccGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCt  >  2:170074/1‑191 (MQ=255)
                                                                                                                                                                         |                                                                                                                                                                                       
ACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTCAGCTCGAACAGCCATTTGACGATGCTCTGCTTCATATCGTCATCTTGCGTGGTATCTTCCGGACCGTTCTCTTTTTCATCGGCCAGGATGTCCAGCAACGCTTTTTCGGAATCACCACCCAGCGGGGTGTCT  >  NC_000913/2866534‑2866886

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