Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 177,251 IS2 (–) +5 bp coding (370‑374/624 nt) yadS ← PF03458 family protein YadS

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 177251 =1 (0.040)26 (0.990) 24/486 0.2 96.2% coding (374/624 nt) yadS PF03458 family protein YadS
?NC_000913 381260 = NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 = 1772551 (0.040)28 (1.050) 27/490 0.2 96.5% coding (370/624 nt) yadS PF03458 family protein YadS
?NC_000913 = 382590 NA (NA)noncoding (1331/1331 nt) IS2 repeat region

GAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGT  >  NC_000913/382348‑382822
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
gAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCGTGTAATGGGTTAAGTGATAACAGATGGCTGGAAATATAGGGGCAAATCCAATTACCTa                                                                                                                                                                                                                                  <  1:250023/251‑1 (MQ=17)
     cGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAgg                                                                                                                                                                                                                             >  1:36604/1‑251 (MQ=34)
     cGGAGAGTAAAGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACAGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATAGCATACGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGGAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACATATCAgg                                                                                                                                                                                                                             >  2:215752/1‑251 (MQ=255)
         gagTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGt                                                                                                                                                                                                                         <  2:229338/251‑1 (MQ=35)
                      gCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                                                                                                                                            <  2:36604/251‑1 (MQ=255)
                        agagagCTTCGTGAAAACGATAAAGCGTGACTACATAAGTATCATGCCCAAACCACACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGACTATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGACTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACt                                                                                                                                                                                                           >  1:161507/1‑250 (MQ=255)
                           gagCTTCGTGAAAACGATAAAGCGTGACTACAACAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGTGAATATCTGCGGCATCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTGCCTATCAGGCAGTTTGCGCCAATCACTAcag                                                                                                                                                                                                       <  2:32554/251‑3 (MQ=255)
                                    gAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGa                                                                                                                                                                                               >  2:112172/1‑250 (MQ=255)
                                                        aCATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATa                                                                                                                                                                          <  1:167529/251‑1 (MQ=255)
                                                             aGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGt                                                                                                                                                                     <  2:205983/251‑1 (MQ=255)
                                                                               gACGGGTTACCGGCAGCAAAGAACCTTGCAGAGACGTTAGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGACTAACGGCAACCTGACAAATACGGGTATTCGCGCCGATaacaat                                                                                                                                                   >  2:184601/1‑250 (MQ=255)
                                                                                          ggCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc                                                                                                                                          >  2:27767/1‑249 (MQ=255)
                                                                                                 aaGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCCGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTg                                                                                                                                                                                                                      >  1:166134/1‑166 (MQ=255)
                                                                                                 aaGAACCGTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTg                                                                                                                                                                                                                      <  2:166134/166‑1 (MQ=255)
                                                                                                          gCAGAGGCTTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCTCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAc                                                                                                                        <  1:91870/251‑1 (MQ=255)
                                                                                                              agGCTTTCGTGCATTATATCGAATGGCTTCCGCATTGTGCGCTGTGTTATCGCTCGCCACGGGAATATCTTCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTTGAAATATAGGGGCAAATCCAATGACCTATCAGGCAGTTGGCGCCAATCTCTACCGGACTCTGAATAACGGCACCATGTCAAATACTTGTATTCTCGCCGATAACAACTGCTCCTGCAGTTTCATCAATAATTTCCTgg                                                                                                                    <  2:228383/251‑1 (MQ=18)
                                                                                                                    tCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGCTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGc                                                                                                               >  2:58142/1‑250 (MQ=255)
                                                                                                                                   aaTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGTCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGTACCCTTAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATAGTTTtc                                                                                               <  2:218548/251‑1 (MQ=255)
                                                                                                                                         aTCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTa                                                                                         <  1:112172/251‑1 (MQ=255)
                                                                                                                                               gtagtgctgttggttctcggtggcCAAGGGAATATGTGCGGCAGCGGGCTGGTAAAGGGTTAAGGGAAAACAGATGTCAGGAAATAAAGGGTCAAATCCAATTACCTATCAGGCAGTGTGCGACAGTCAATACCGGACACTGAATAACGGCACCATGACAAATAGGGGTAGTCTCGCCGATAACAACTTCTCCTGCAGTTTCATCAATAATTACCTGGTCTGCACTTTGAATAGTTTGCTCTTTAACagta                                                                                   <  2:300/229‑4 (MQ=255)
                                                                                                                                                     cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTc                                                                             >  1:65307/1‑251 (MQ=255)
                                                                                                                                                          ggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGTCTTGTAAAGTGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGGTTCATCAATAAGTACCTGGTCTGCAATCGGAATATTTTTCTCTTTATCATTAAGCTTCAGATa                                                                        <  2:220291/251‑1 (MQ=38)
                                                                                                                                                                 gcccgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTc                                                                 <  1:2866/248‑1 (MQ=255)
                                                                                                                                                                                                       gggTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACg                           >  1:230054/1‑251 (MQ=255)
                                                                                                                                                                                                              gTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTGCCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACTGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCATTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCCTccac                    <  1:53479/251‑1 (MQ=255)
                                                                                                                                                                                                               tGATAACAGATGTCTGGAAATATAGGGGCAACTCCAATTACCTATCAGGCAGTTTGCGCCAATCAATACCGGACCCTGAATAACGGCACCATGACAACTACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTccacc                   >  2:136715/1‑251 (MQ=255)
                                                                                                                                                                                                                       gATGTCTGGAAATATAGGGGCAAATCCAATTACCTATAAGGCAGTTTGCGCCAATCAATACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGATTCATCAATAATTACCTGGTCTGCAATCTGACTATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCgg            >  1:172186/1‑250 (MQ=255)
                                                                                                                                                                                                                                aaaTATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggt  <  1:136715/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
GAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGT  >  NC_000913/382348‑382822

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