Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 362,204 IS1 (–) +9 bp coding (968‑976/1254 nt) lacY ← lactose permease

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 19796 =NA (NA)46 (0.970) 43/486 0.2 100% noncoding (768/768 nt) IS1 repeat region
?NC_000913 = 362212 0 (0.000)coding (968/1254 nt) lacY lactose permease
* ? NC_000913 = 258675NA (NA)49 (1.030) 47/486 0.1 100% noncoding (1/768 nt) IS1 repeat region
?NC_000913 362204 = 0 (0.000)coding (976/1254 nt) lacY lactose permease

AAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTG  >  NC_000913/19641‑20043
                                                                                                                                                           |                                                                                                                                                                                                                                                       
aaGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGTAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACtt                                                                                                                                                          >  1:56949/1‑251 (MQ=255)
aaGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACtt                                                                                                                                                          >  1:237063/1‑251 (MQ=255)
       gggAGCTAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGAtttt                                                                                                                                                                           >  2:333464/1‑227 (MQ=255)
       gggAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGAtttt                                                                                                                                                                           <  1:333464/227‑1 (MQ=255)
                                                                                           ggCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGTAAAGAGTAGTAATTGCTGGTAATGTCTCCAACTGATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTccc                                                               <  2:56949/251‑1 (MQ=255)
                                                                                           ggCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTccc                                                               <  2:237063/251‑1 (MQ=255)
                                                                                                                                   aagcgaatgtctttctgaagagatGGTAATGCCTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTGTGTCATGCAGCGCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATGCATACAGCGGCCAGCCATCCGTCATcc                       <  1:288861/228‑1 (MQ=11)
                                                                                                                                            tctttctggagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCACGATGACTTTGTCATGCAGCTCCACCGTTTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATccacc              >  2:81066/15‑247 (MQ=11)
                                                                                                                                                tctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCa                                                                       >  1:12893/11‑190 (MQ=2)
                                                                                                                                                tctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAACGCCCGATGCCTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCa                                                                       <  2:12893/180‑1 (MQ=255)
                                                                                                                                                   gaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGCGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCCCGTCaa        >  1:254075/8‑250 (MQ=2)
                                                                                                                                                    tctcgttGGTAATGGCTCCACCTTATTGGTAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCCTGCAGCTCCACCGATTTTGCGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGCTTATGCCGCTCAATTCGCTGCGACTATCGCTTGCTGCTTACGTGCAGCTTTCCCTTCAGGCGGGCTTGATTCAGCTGCCCGCCATCCGTCATCCATCTCCACACGTCaccg      >  2:258731/7‑248 (MQ=11)
                                                                                                                                                       cgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAggg  >  1:115528/4‑250 (MQ=11)
                                                                                                                                                        gCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCACTTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTg  >  2:343635/1‑251 (MQ=14)
                                                                                                                                                           |                                                                                                                                                                                                                                                       
AAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTG  >  NC_000913/19641‑20043

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