Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,353,837 IS2 (+) +5 bp coding (877‑881/1518 nt) lysU ← lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3186096 =NA (NA)21 (0.810) 18/490 0.5 83.7% noncoding (1/1331 nt) IS2 repeat region
?NC_000913 = 4353841 4 (0.160)coding (877/1518 nt) lysU lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase
* ? NC_000913 4353837 =5 (0.200)21 (0.810) 21/492 0.3 80.5% coding (881/1518 nt) lysU lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase
?NC_000913 = 4499511 NA (NA)noncoding (1331/1331 nt) IS2 repeat region

GTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTC  >  NC_000913/4499264‑4499721
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                  
gTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCt                                                                                                                                                                                                                 <  1:112241/251‑1 (MQ=11)
        tCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTg                                                                                                                                                                                                          >  1:155301/1‑250 (MQ=33)
           ggAGAGAAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGatat                                                                                                                                                                                                      >  2:87205/1‑251 (MQ=25)
              gagTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGaa                                                                                                                                                                                                   <  2:232671/251‑1 (MQ=37)
                                 agCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCAGAGTTCGCTGGGTTATCGCTCGCCACGGGAATATCTGTGGCTGCGGGCGTGTAATTGGTTAAGTGATAACAGATGGGTGGAAATATAGGGGCAAATTGACTCTGTGCTGATATGAATACGGATATTTTATACAGa                                                                                                                                                                                <  1:239564/251‑1 (MQ=14)
                                   cTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGCTAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAAt                                                                                                                                                                              >  1:125836/1‑251 (MQ=255)
                                     tCGTGAAAACGATAAAGCTTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                                                                                                                             <  2:48264/250‑1 (MQ=255)
                                             aCGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGAGATTTTATACAGAATACTTTTCAGt                                                                                                                                                                    <  1:85511/251‑1 (MQ=255)
                                               gATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTAc                                                                                                                                                                  >  2:91666/1‑251 (MQ=255)
                                                  aaagcccgACTACATCAGTAACATGACCAAACCAGACGGGTTAACGGCAGCAAAGAACCTCGCAGACGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGAGGCAGCGGGCTTGTAATGGGTTACGTGCTAACAGATGTCTGGAAATATAGGGGCAACTCCACCTTGTGATGATATGAATACGGATATTTTATACAGAATACTTTTCAGTaaaa                                                                                                                                                                >  2:121672/8‑247 (MQ=17)
                                                      cGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGAc                                                                                                                                                           <  2:122247/251‑1 (MQ=255)
                                                               aTCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTAc                                                                                                                                                   >  2:267474/1‑250 (MQ=255)
                                                                 cAGTATCATGCCCAAACCAGACGGGTTAACGGCCGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                                                                                                                                >  2:252691/1‑251 (MQ=255)
                                                                       cATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                                                                                                                           <  2:155301/250‑1 (MQ=255)
                                                                            ccAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCatatat                                                                                                                                     <  1:116158/251‑1 (MQ=255)
                                                                                                    caAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGTCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGAGCAGGTAGAATa                                                                                                             <  2:207390/251‑1 (MQ=255)
                                                                                                                                        aaTGGCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTc                                                                         <  2:175996/251‑1 (MQ=255)
                                                                                                                                                 cGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTccc                                                                <  1:64971/251‑1 (MQ=255)
                                                                                                                                                                                                  ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGGTATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATg               >  2:132481/1‑251 (MQ=255)
                                                                                                                                                                                                                tAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTc  <  2:125025/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                  
GTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTC  >  NC_000913/4499264‑4499721

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