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)63 (1.040) 56/482 0.1 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)65 (1.080) 62/482 0.1 100% noncoding (1/768 nt) IS1 repeat region
?NC_000913 362204 = 0 (0.000)coding (976/1254 nt) lacY lactose permease

GAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTG  >  NC_000913/19571‑20043
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
gAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAgt                                                                                                                                                                                                                                 <  1:124123/251‑1 (MQ=255)
                                                                             gggAGATAGATTAGTTGTACATTACCACGATTTTGACTCAGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGCATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCTCCGATTTTGAGAACGACAGCGACTTCCGTccc                                                                                                                                                    >  2:347380/1‑251 (MQ=255)
                                                                                                                                         agaCATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGc                                                                                        >  1:170236/1‑251 (MQ=255)
                                                                                                                                                                                                       ggATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCCACTTATTGATAGTGTTTTCTGTTCAGATAATGCCCGATGAATTTGTCATGCAACTCCACCGATTTTGAGAACGACAGCGCCTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATCCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGctg                                                                       >  1:376245/1‑204 (MQ=14)
                                                                                                                                                                                                       ggATGCTGAAAAGAGAGGTAATTGCTGGTAATGACTCCAACTTATTGAGAGTGTTTTATGTTCAGATAATGCCCGATGTCTTTGTCATGCAACCCCACCGATGTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCGGCCTCAGATCCAGGTCATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACTTGCat                                                                       <  2:376245/206‑2 (MQ=14)
                                                                                                                                                                                                       acaagcgaatggcgttcggaagagatGGTAAGGACTCCAACTTAGTAATAGTGTTTTATGTTCAGATAAGGCCCGATTACTTTGTCATGCAGCGCCACCGATGTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATGCGCTGCGTATATCGCGTGCTGATTACGTGCAGCTTTCCCGTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCAt                         <  2:327428/226‑1 (MQ=11)
                                                                                                                                                                                                        caagcgaatgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATc                        >  1:314102/25‑251 (MQ=2)
                                                                                                                                                                                                           gcgaatgtctttctgaagagatGGTACTGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCat                     >  1:590497/22‑251 (MQ=255)
                                                                                                                                                                                                             gaatgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGTTTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGAGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCTGGAGTCCT‑ACAGCGGCCAGCCATCCGTCATCCatat                   >  2:193308/20‑251 (MQ=11)
                                                                                                                                                                                                             gaatgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTAATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCACTTCAGGCGGGATTCAT‑ACAGCGGCCAGACATCCGTCATCCatct                   >  1:255358/20‑249 (MQ=11)
                                                                                                                                                                                                              atatatagtccctctcgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATc                  >  1:163712/19‑251 (MQ=2)
                                                                                                                                                                                                              aatgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATc                  >  2:390538/19‑251 (MQ=2)
                                                                                                                                                                                                               aaaaGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCTCCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATca                 <  1:347380/251‑1 (MQ=38)
                                                                                                                                                                                                                tgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATcac                >  1:126632/17‑251 (MQ=2)
                                                                                                                                                                                                                    agtccctctcgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGTTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACg            <  1:27839/239‑1 (MQ=2)
                                                                                                                                                                                                                        cctctcgttGGTAATGACTCCAACTTCTTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCCCAATTCGATGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTAAGGCGGGATTCATTAC‑GCGGCCAGCACTCCGTCCTCCATATCCCCCCGTCa         >  2:544541/9‑250 (MQ=11)
                                                                                                                                                                                                                           ctcgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAgg     >  1:455342/6‑251 (MQ=11)
                                                                                                                                                                                                                             cgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGcc                                                                                                                                                >  1:401685/4‑111 (MQ=11)
                                                                                                                                                                                                                             cgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCGCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGcc                                                                                                                                                <  2:401685/108‑1 (MQ=2)
                                                                                                                                                                                                                              gttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGgcca                                     <  2:281558/214‑1 (MQ=11)
                                                                                                                                                                                                                              gttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGgcca                                     >  1:281558/3‑216 (MQ=11)
                                                                                                                                                                                                                               atGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTg  >  1:135542/2‑250 (MQ=11)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        
GAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTG  >  NC_000913/19571‑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.