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

GCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTCAGTTTATTGTCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGCCTTTTTGAATGAACTTGCTCGTCTGGTGGGTTCTTCACACCTGCTCAC  >  CP009273/2215490‑2215743
                                                                                                                                 |                                                                                                                            
gCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGCGAGAGAAAAAGAATGCTCAGTTTATTGTCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTCTAGTTTGTGa                                                                                                               >  1:180449/1‑145 (MQ=255)
       cattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcggcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCTATTCTAGTTTGTGATAttttt                                                                                                        >  2:317136‑M2/123‑145 (MQ=255)
              cGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTCAGTTTATTGTCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGACGCTATTCTAGTTTGTGATATTTTTTCGccac                                                                                                 >  2:336786/1‑145 (MQ=255)
                   tggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCTATTCTAGTTTGTGATATTTTTTCGCCACCACaa                                                                                            <  2:279041‑M2/35‑1 (MQ=255)
                                     tttGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTCAGTTTATTGTCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTtagactggccccctgaatctccagacaaccaatatcactaaaataagt                                                                          <  1:122688‑M2/145‑49 (MQ=255)
                                     ctgggttatcgctcgccacggcaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTc                                                                          <  1:115727‑M2/53‑1 (MQ=255)
                                                        agagagaAAAAGAATGCTCAGTTTATTGTCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGataat                                                       >  2:266398/1‑145 (MQ=255)
                                                             gaAAAAGAATGCTCAGTTTATTGTCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTCTAGCTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGc                                                  >  2:226986/1‑145 (MQ=255)
                                                                ctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCCTGACAACAACAGAAAATAAAACctg                                               >  1:241802‑M2/66‑144 (MQ=255)
                                                                      cagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGCCTTTTtgaa                                         <  2:261563‑M2/86‑1 (MQ=255)
                                                                              ttATTGTCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTtagactggctccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaagactagtttttagactagtcattggagaac                                 >  2:335892‑M2/1‑56 (MQ=255)
                                                                               tATTGTCTCAATTTACAAAATAATCACTCGCTGAATTGTTATACAAGGCGCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGCCTTTTTGAATGAACttac                                <  2:13781/145‑3 (MQ=255)
                                                                                    tCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCTCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGCCTTTTTGAATGAACTTGCTCGTc                           <  1:257586/145‑1 (MQ=255)
                                                                                    tCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTCTAGTTTGTGATATTTTTTCGTCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGCCTTTTTGAATGACCTTGCTCGTc                           <  2:87275/145‑1 (MQ=255)
                                                                                       gAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtctttagactagacagtggagcacggatgattg                        <  1:120558‑M2/145‑99 (MQ=255)
                                                                                       gAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGCCTTTTTGACTGAACTTGATCGTCtgg                        >  2:54523/1‑145 (MQ=255)
                                                                                           tttCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGCCTTTTTGAATGCACTTGCTCGTCTGGTggg                    <  2:167804/145‑1 (MQ=255)
                                                                                                             gaaatataggggcaaatccaCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGCCTTTTTGAATGAACTTGCTCGTCTGGTGGGTTCTTCACACCTGCTCAc  <  2:14129‑M2/125‑1 (MQ=255)
                                                                                                                                 |                                                                                                                            
GCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTCAGTTTATTGTCTGAATTTTCAAAATATTCACTCGCTGAATTGTTATACAAGGCGCTATTCTAGTTTGTGATATTTTTTCGCCACCACAAGGAGTGGAAAATGTCTTCCATGACAACAACTGATAATAAAGCCTTTTTGAATGAACTTGCTCGTCTGGTGGGTTCTTCACACCTGCTCAC  >  CP009273/2215490‑2215743

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