Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC CP009273 2,215,619 IS2 (+) +5 bp 33.3% intergenic (‑151/‑41) bglX ← / → dld beta‑D‑glucoside glucohydrolase, periplasmic/D‑lactate dehydrogenase, FAD‑binding, NADH independent

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 2062422 =NA (NA)6 (0.260) 6/268 NT 39.1% noncoding (1331/1331 nt) IS2 repeat region
?CP009273 2215619 = 9 (0.400)intergenic (‑151/‑45) bglX/dld beta‑D‑glucoside glucohydrolase, periplasmic/D‑lactate dehydrogenase, FAD‑binding, NADH independent
* ? CP009273 = 2063752NA (NA)3 (0.130) 3/268 NT 24.3% noncoding (1/1331 nt) IS2 repeat region
?CP009273 = 2215623 9 (0.400)intergenic (‑155/‑41) bglX/dld beta‑D‑glucoside glucohydrolase, periplasmic/D‑lactate dehydrogenase, FAD‑binding, NADH independent

GACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATT  >  CP009273/2062286‑2062533
                                                                                                                                        |                                                                                                               
gACTATACAGGTTATTGACCAACAAACATTGCAGAACTCGACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACGCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGc                                                                                                         >  1:258291/1‑145 (MQ=255)
                gACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGCGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTcacgac                                                                                         >  1:105114/1‑140 (MQ=255)
                                     tCCACCAATGCTTGCATGGCCTATAATTGGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCAt                                                                    <  1:136411/145‑1 (MQ=38)
                                     tCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCAt                                                                    >  1:369935/1‑145 (MQ=255)
                                        accaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGAGTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTAc                                                                 >  2:171723/1‑145 (MQ=255)
                                                             aaTTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATAtt                                            <  2:406334/145‑1 (MQ=255)
                                                                ttGCAGGACATCCCCGGAGGAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                         <  1:291497/145‑1 (MQ=255)
                                                                 tGCAGGACATCCCCGGAGTAGAAATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCg                                        <  2:157434/145‑1 (MQ=255)
                                                                   cAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTg                                      <  2:369935/145‑1 (MQ=255)
                                                                             ccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                            >  1:381383/1‑145 (MQ=255)
                                                                             ccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                            >  2:182764/1‑145 (MQ=255)
                                                                                                    tCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCGCGCGCACTATGCGGATggc     >  1:96539/1‑143 (MQ=255)
                                                                                                      cTTGGCAGGCCGGAGACGAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCATCGCACTATGCGGATGCCAt   <  1:171723/145‑1 (MQ=255)
                                                                                                       ttGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCACCGCACTATGCGGATGCCCtt  >  1:347868/1‑145 (MQ=255)
                                                                                                                                        |                                                                                                               
GACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATT  >  CP009273/2062286‑2062533

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ 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.