Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,044,847 IS5 (+) +4 bp coding (649‑652/1431 nt) yihF → DUF945 domain‑containing protein YihF

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3652036 =NA (NA)22 (0.680) 21/488 0.5 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 4044847 = 0 (0.000)coding (649/1431 nt) yihF DUF945 domain‑containing protein YihF
* ? NC_000913 = 3653230NA (NA)44 (1.360) 36/488 0.1 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 = 4044850 0 (0.000)coding (652/1431 nt) yihF DUF945 domain‑containing protein YihF

ATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGACATATTAAAATATCAGCAAGAATTC  >  NC_000913/3652998‑3653455
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                 
aTAAGGACGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCCTGCGGGACAAGAACATCTATTTTCTGGTCTGACGGCGCTTACTGATGAATTCACTGTCGGCGAAGGTACGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACACACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAttttt                                                                                                                                                                                                                 >  1:37447/1‑251 (MQ=11)
  aaGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAttttttt                                                                                                                                                                                                               >  1:287108/1‑251 (MQ=35)
       cGCCGGCCATTACCAGCCTTGGTGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGTCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCATATGTCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATa                                                                                                                                                                                                          <  2:313276/251‑1 (MQ=255)
       cGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGTCTCCAGATTACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACAt                                                                                                                                                                                                           <  1:15808/250‑1 (MQ=255)
           ggCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGa                                                                                                                                                                                                      <  1:262855/251‑1 (MQ=40)
                    cAGCCTTGGGGTAAAACGGCTCGTTGTCTTCCACCATGTTTTGCCATGGCAGAATCGGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCGCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACGCTATTTTTCTACATAAAGAACACCCTAc                                                                                                                                                                                             <  1:106620/251‑1 (MQ=255)
                      gCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAATGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCGCCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccc                                                                                                                                                                                           <  1:131689/251‑1 (MQ=255)
                      gCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATTACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccc                                                                                                                                                                                           <  1:287644/251‑1 (MQ=255)
                        cTTGGGGTAAAAAGGCTCGATGACTTCCCACATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGGGACACTATTTTTTTACATAAAGAACAACCTAccccc                                                                                                                                                                                         >  2:272205/1‑251 (MQ=255)
                                              aCTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAat                                                                                                                                                                   >  1:179538/1‑251 (MQ=255)
                                                   caccaTGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCtt                                                                                                                                                              <  1:283973/251‑1 (MQ=255)
                                                             tGCCATGGCAGAATCTGCCCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGTACAATATCTTGacaaacaa                                                                                                                                                     <  1:268180/250‑1 (MQ=255)
                                                                             gCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGAGACTTGTTCGCCCCTTCCTTAGGTGACACGATTTTTTTACATAAATAACACCCTACCTCGAAATTATATAAAGAGGGTGAATATCTGGATAAACAAATCTAAGGGTTAtttt                                                                                                                                     <  2:62829/250‑1 (MQ=17)
                                                                                         aCAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGCTGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAAc                                                                                                                        >  1:308544/1‑251 (MQ=255)
                                                                                            aGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAt                                                                                                                     >  1:119786/1‑251 (MQ=255)
                                                                                                aaTCTCTTTTCTGGTCTGACGGCGCTTAATGCTGAATTCCCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAAACCTGTTCTATGGCTCAAGATGACAAACATGATCTCATATAAGGGACTTTTTCGAACATTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAACCCCAAATTATATAAAGATGGACAATATCTTGACAAAAAAATCGAAGGTTTATTTTAAACATTATTAACCATTAAt                                                                                                                 >  1:117646/1‑251 (MQ=21)
                                                                                                              tcggACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCGACACTATTTTTTTACATAAAGAacac                                                                                                                                                                                                  >  2:86305/4‑156 (MQ=255)
                                                                                                              tCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAacac                                                                                                                                                                                                  <  1:86305/156‑1 (MQ=255)
                                                                                                                      cgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGCAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCaaa                                                                                             >  2:24376/1‑249 (MQ=255)
                                                                                                                                                                               gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCAt                                  <  2:154523/251‑1 (MQ=255)
                                                                                                                                                                                   tATGGCTCCAGATGACAAACATGGTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACa                               >  2:219714/1‑250 (MQ=255)
                                                                                                                                                                                    atgtcTCCAGATGACAAACATGATCTCATATCAGGGACTGTTTCGCATCTTCCTTATGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTTACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCGTCTACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAAt                             <  1:215997/247‑1 (MQ=255)
                                                                                                                                                                                                               atatCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGACATATTAAAATATCAGCAAGAATTc  <  2:287108/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATtatata                                                                                                                                                                               >  2:124316/1‑61 (MQ=255)
                                                                                                                                                                                                                                gCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATtatata                                                                                                                                                                               <  1:124316/61‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                 
ATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGACATATTAAAATATCAGCAAGAATTC  >  NC_000913/3652998‑3653455

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.