Predicted mutation
evidence seq id position mutation freq annotation gene description
JC NC_007779 2,865,278 +TTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGC 20.5% coding (813/876 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_007779 2865229 =157 (0.790)20 (0.190) 16/96 NT 20.5% coding (862/876 nt) rpoS RNA polymerase, sigma S (sigma 38) factor
?NC_007779 = 2865278 133 (0.670)coding (813/876 nt) rpoS RNA polymerase, sigma S (sigma 38) factor

CGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_007779/2865266‑2865229
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCG  <  NC_007779/2865278‑2865235
                                      ||||||||||||||||||||||||||||||||||||||||||                                            
CGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCT                         >  1:8629409/1‑101
 GAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTG                        >  1:6744074/1‑101
     TCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAA                    <  1:5655829/101‑1
      CCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAAC                   >  1:6920789/1‑101
         GCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCA                <  1:865756/101‑1
          CAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAG               >  1:2757506/1‑101
           AAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGACTATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGG              >  1:7866588/1‑101
           AAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGG              >  1:7008192/1‑101
           AAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGG              >  1:6389817/1‑101
              CGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGC           >  1:7183509/1‑101
              CGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGC           <  1:1513358/101‑1
               GCAGGGGCTGAATATCGAAGCGCCGTTCGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCT          <  1:5321605/101‑1
                CAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTG         <  1:1978623/101‑1
                 AGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGA        >  1:2605231/1‑101
                  GGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAA       <  1:6750926/101‑1
                    GGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATA     <  1:6732631/101‑1
                    GGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATA     <  1:4895362/101‑1
                    GGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATA     <  1:8515183/101‑1
                      CTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATC   >  1:9307499/1‑101
                       TGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGTAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCG  >  1:6263745/1‑101
                                      ||||||||||||||||||||||||||||||||||||||||||                                            
CGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_007779/2865266‑2865229
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCG  <  NC_007779/2865278‑2865235

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