Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 3,367,990 IS2 (–) +5 bp coding (164‑168/1128 nt) yhcG → DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2996361 =NA (NA)30 (0.810) 22/260 0.9 100% noncoding (1331/1331 nt) IS2 repeat region
?NC_000913 = 3367994 0 (0.000)coding (168/1128 nt) yhcG DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
* ? NC_000913 3186096 =NA (NA)50 (1.350) 40/260 0.1 100% noncoding (1/1331 nt) IS2 repeat region
?NC_000913 3367990 = 0 (0.000)coding (164/1128 nt) yhcG DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily

TTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCT  >  NC_000913/2996219‑2996495
                                                                                                                                              |                                                                                                                                      
ttatttCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGAttt                                                                                                                                  >  3:150991/1‑149 (MQ=255)
        aGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCta                                                                                                                           >  1:28535/1‑148 (MQ=255)
            gTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATAttt                                                                                                                      <  8:16156/149‑1 (MQ=255)
                  cTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGa                                                                                                                 >  1:126797/1‑148 (MQ=255)
                  cTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGa                                                                                                                 >  3:75232/1‑148 (MQ=255)
                     ataaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATc                                                                                                             <  4:75232/149‑1 (MQ=255)
                        aataatCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGt                                                                                                          <  4:1142/149‑1 (MQ=255)
                                cTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATc                                                                                                      <  2:90893/145‑1 (MQ=255)
                                cTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATc                                                                                                      >  1:90893/1‑145 (MQ=255)
                                       cATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCAt                                                                                           <  2:28535/149‑1 (MQ=255)
                                              aaCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGc                                                                                    <  4:150991/149‑1 (MQ=255)
                                                                        cAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTg                                                             <  2:126797/146‑1 (MQ=255)
                                                                           aCCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                                                <  4:79696/140‑1 (MQ=255)
                                                                           aCCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                                                >  3:79696/1‑140 (MQ=255)
                                                                                                  ccACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCt                                                                               >  4:9360/1‑102 (MQ=255)
                                                                                                  ccACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCt                                                                               <  3:9360/102‑1 (MQ=255)
                                                                                                       gAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCa                           >  1:53358/1‑149 (MQ=255)
                                                                                                        aGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                          >  3:46065/1‑149 (MQ=255)
                                                                                                            aaCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGt                      >  3:123109/1‑149 (MQ=255)
                                                                                                             aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGtt                     <  2:53358/149‑1 (MQ=255)
                                                                                                                    aaTGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGc              >  8:26869/1‑149 (MQ=255)
                                                                                                                      tGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGAt                               <  2:163754/130‑1 (MQ=255)
                                                                                                                      tGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGAt                               >  1:163754/1‑130 (MQ=255)
                                                                                                                         aaGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa         <  4:46065/149‑1 (MQ=255)
                                                                                                                          aGTAAATTCCAATTGATTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAAc        >  5:109577/1‑149 (MQ=255)
                                                                                                                                ttCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCctct  >  2:129980/1‑149 (MQ=39)
                                                                                                                                    aaTTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                          >  2:39210/1‑121 (MQ=37)
                                                                                                                                              |                                                                                                                                      
TTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCT  >  NC_000913/2996219‑2996495

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