Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,110,680 25702 bp→GA [yghD]–pitB 20 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3110680 3136381 25702 34 [0] [0] 33 [yghD]–pitB [yghD],yghE,yghF,yghG,pppA,yghJ,glcA,glcB,glcG,glcF,glcE,glcD,glcC,yghO,insH‑9,yghQ,yghR,yghS,yghT,pitB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 31106790 (0.000)32 (0.890)
+GA		 28/94 0.1 100% coding (448/537 nt) yghD putative type II secretion system M‑type protein
?NC_000913 3136382 = 0 (0.000)intergenic (‑11/+281) pitB/gss metal phosphate:H(+) symporter PitB/fused glutathionylspermidine amidase/glutathionylspermidine synthetase

AGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3110632‑3110679
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGG  >  NC_000913/3136382‑3136429
                                                ||                                                
AGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGA                                                  <  1:1795356/50‑1
AGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAG                                                   >  1:659632/1‑49
 GGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAA                                                 >  1:1141070/1‑50
 GGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAA                                                 <  1:52752/50‑1
     TCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAACGC                                              >  1:2280728/1‑49
       TCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCAT                                           >  1:290061/1‑50
         AGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTT                                         >  1:2124164/1‑50
         AGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCAT                                           <  1:975614/48‑1
          GCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTT                                        <  1:216896/50‑1
          GCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTT                                        <  1:2137586/50‑1
            AGCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGA                                      <  1:1020641/50‑1
             GCACTGACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAA                                     <  1:2233771/50‑1
              CACTGACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAAT                                    >  1:986911/1‑50
                  GACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAATTTAC                                >  1:1786332/1‑50
                   ACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAATTTACG                               <  1:1212059/50‑1
                   ACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAATTTAC                                <  1:2399943/49‑1
                   ACATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAATTTAC                                >  1:511842/1‑49
                    CATCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAATTTACGA                              <  1:1145217/50‑1
                      TCAATTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAATTTACGAAC                            <  1:601976/50‑1
                          TTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAATTTACGAACGGAC                        >  1:1828422/1‑50
                          TTTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAATTTACGAACGGAC                        >  1:199679/1‑50
                           TTGTGTCACCCGCAGCGCATAGAACGCCATTTTGAATTTACGAACGGACG                       >  1:2496603/1‑50
                                 CACCCGCAGCGCATAGAACGCCATTTTGAATTTACGAACGGACGCATTAT                 >  1:1331414/1‑50
                                  ACCCGCAGCGCATAGAACGCCATTTTGAATTTACGAACGGACGCATTATC                <  1:2412580/50‑1
                                  ACCCGCAGCGCATAGAACGCCATTTTGAATTTACGAACGGACGCATTAT                 >  1:2613087/1‑49
                                   CCCGCAGCGCATAGAACGCCATTTTGAATTTACGAACGGACGCATTAT                 >  1:501572/1‑48
                                     CGCAGCGCATAGAACGCCATTTTGAATTTACGAACGGACGCATTATCAGT             >  1:474687/1‑50
                                       CAGCGCATAGAACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGA           >  1:1370347/1‑50
                                       CAGCGCATAGAACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGA           >  1:1851915/1‑50
                                           GCATAGAACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTT        >  1:2765407/1‑49
                                              TAGAACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAC    >  1:2697861/1‑50
                                              TAGAACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAC    <  1:983613/50‑1
                                               AGAACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACG   >  1:926978/1‑50
                                               AGAACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACG   >  1:969013/1‑50
                                                GAACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGG  <  1:1722554/50‑1
                                                ||                                                
AGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3110632‑3110679
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGG  >  NC_000913/3136382‑3136429

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