Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 896,463 IS150 (+) +3 bp 100% coding (159‑161/432 nt) ECB_00838 → putative phage tail protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 666130NA (NA)17 (0.640) 15/278 NT 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 896463 = 0 (0.000)coding (159/432 nt) ECB_00838 putative phage tail protein
* ? REL606 = 8964650 (0.000)29 (1.100) 23/274 NT 100% coding (161/432 nt) ECB_00838 putative phage tail protein
?REL606 = 2775877 NA (NA)noncoding (1/1443 nt) IS150 repeat region

ACAATCCCTGCGCCACGCCCTCAATAAAGCGGTGCCTTATGTCCGCAATAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAACCCGGCGTTACGGGAAAAACTATTTACCTTTGAGGTGGATATTCTGCGCAACGATGTCTGTGA  >  REL606/896316‑896612
                                                                                                                                                   |                                                                                                                                                     
acaATCCCTGCGCCACGCCCTCAATAAAGCGGTGCCTTATGTCCGCAATAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGAt                                                                                                                                                     >  2:605968/1‑150 (MQ=255)
 gaactgaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGatc                                                                                                                                                    <  1:179711‑M2/4‑2 (MQ=255)
 gaactgaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGAtt                                                                                                                                                    <  2:68830‑M2/4‑1 (MQ=255)
 gaactgaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGAtt                                                                                                                                                    <  1:232235‑M2/4‑1 (MQ=255)
     cccTGTGCCACGCCCTCAATAAAGCGGTGCCTTATGTCCGCAATAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtac                                                                                                                                                <  1:218787‑M1/150‑6 (MQ=255)
     cccTGCGCCACGCCCTCAATAAAGCGGTGCCTTATGTCCGCAATAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtac                                                                                                                                                >  2:551613‑M1/1‑145 (MQ=255)
      gaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaagggctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaacttttttgggtcagtacaGATTGAGGa                                                                                                                                               >  1:311412‑M2/142‑150 (MQ=255)
                          atattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATc                                                                                                                           <  1:64109‑M2/29‑1 (MQ=255)
                          atattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATc                                                                                                                           >  2:64300‑M2/122‑150 (MQ=255)
                          atattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaataagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATc                                                                                                                           >  2:109550‑M2/122‑150 (MQ=255)
                               gTGCCTTATGTCCGCAATAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaa                                                                                                                      >  2:50668‑M1/1‑119 (MQ=255)
                                 gCCTTATGTCCGCAATAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatg                                                                                                                    >  2:223218‑M1/1‑117 (MQ=255)
                                  ccTTATGTCCGCAATAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatgg                                                                                                                   <  2:40710‑M1/150‑35 (MQ=255)
                                           agcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGtt                                                                                                          >  2:416387‑M2/105‑150 (MQ=255)
                                               aTAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaat                                                                                                      <  2:117829‑M1/150‑48 (MQ=255)
                                                tAACCCGGACAAACTGCTTCTGGTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaata                                                                                                     <  1:604194‑M1/150‑49 (MQ=255)
                                                tAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTTCACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaata                                                                                                     <  2:432947‑M1/150‑49 (MQ=255)
                                                           aaCTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagcca                                                                                          >  2:501888‑M1/1‑91 (MQ=255)
                                                           aaCTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGATACACCATGAACGCGGGGATagtactgcacccaatttgttggacgaggaagtggactcgcccatactatatcaatgcca                                                                                          >  1:282926‑M1/1‑91 (MQ=255)
                                                             cTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaa                                                                                        >  1:457420‑M1/1‑89 (MQ=255)
                                                                 taaaaggtctgactccaattgaatatcggaatcagacttatatgcatcgtgtttaactgtccaacgttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAAcc                                                                                    <  1:76538‑M2/68‑1 (MQ=255)
                                                                 aTCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatac                                                                                    <  2:591606‑M1/150‑66 (MQ=255)
                                                                      ggtctgactccaattgaatatcgggctcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCGGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGccc                                                                               >  1:20884‑M2/78‑150 (MQ=255)
                                                                         gtgGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttga                                                                            >  1:27336‑M1/1‑77 (MQ=255)
                                                                           gtctccaatggaatatcggaatcatacctatatacctcgtgtttaagtgtccaactttttgggggcagtacaGATTGAGGATTTCAGCGGCGACCAGATTCTGCTGATGGACCCGGTTTTGCTGGGGCTGCGGGATAACCAGCCCGATGc                                                                          <  1:656450‑M2/78‑1 (MQ=255)
                                                                           ggATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcacctaatatggcaaagccaaaataccattttgaac                                                                          >  1:601222‑M1/1‑75 (MQ=255)
                                                                                  attgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAAt                                                                   <  2:561379‑M2/85‑1 (MQ=255)
                                                                                      aatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggtgtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAAcc                                                               <  1:19834‑M2/89‑1 (MQ=255)
                                                                                       atatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAAccc                                                              <  2:529251‑M2/90‑1 (MQ=255)
                                                                                        cTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagt                                                             <  2:458722‑M1/150‑89 (MQ=255)
                                                                                               ccACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGACCGCGGTGATtgtactgcatccattttgttggacgatgaaatgcaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaat                                                      <  1:550020‑M1/150‑96 (MQ=255)
                                                                                               ccACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagagcctttaagtcgtgaat                                                      >  1:515449‑M1/1‑55 (MQ=255)
                                                                                                     gggCCCGTTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgtagaacgatgaaatggaatagtcccgaagatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaaacactac                                                <  1:222581‑M1/150‑102 (MQ=255)
                                                                                                      ggCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactact                                               >  1:40318‑M1/1‑48 (MQ=255)
                                                                                                          aGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcac                                           >  2:351528‑M1/1‑44 (MQ=255)
                                                                                                                 gtgtttaactgtccaactttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAACCCGGCGTTACGGGAAAAACTATTTACCt                                    >  2:395101‑M2/35‑150 (MQ=255)
                                                                                                                 gtgtttaactgtccaactttttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAACCCGGCGTTACGGGAAAAACTATTTACCt                                    >  2:5026‑M2/35‑150 (MQ=255)
                                                                                                                  gTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatg                                   <  2:602991‑M1/150‑115 (MQ=255)
                                                                                                                      tGGGAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatgg                               >  2:411917‑M1/1‑32 (MQ=255)
                                                                                                                       gggAGTACCGTTACTCCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgcctggaagtcgtgaatcactacttcacaactgatgatggt                              <  1:500458‑M1/150‑120 (MQ=255)
                                                                                                                         gAGTACCGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggtta                            >  2:184300‑M1/1‑29 (MQ=255)
                                                                                                                           gtccaactctttggggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATCACCAGCCCGATGCCATCAATAACCCGGCCTTATGGGATAAACTATTTACCTTTGAGGTGGa                          >  1:297435‑M2/25‑150 (MQ=255)
                                                                                                                              ccGTTCCGCCCTGAACGCGGTGATtgtactgaacccattttgttggacgatgaaatggaatagcccctaatatgtcaaatctaaaatacccttttgaaaagcgcctttaagtcgtgaatcactacttcacaactgatgatggttacagga                       >  1:517318‑M1/1‑24 (MQ=255)
                                                                                                                              ccGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacagga                       >  2:200025‑M1/1‑24 (MQ=255)
                                                                                                                              ccGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacagga                       >  2:610726‑M1/1‑24 (MQ=255)
                                                                                                                              ccGTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacagga                       >  2:126502‑M1/1‑24 (MQ=255)
                                                                                                                                gTTACACCCTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatc                     >  2:27373‑M1/1‑22 (MQ=255)
                                                                                                                                      cccTGAACGCGGTGATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcg               >  2:276557‑M1/1‑16 (MQ=255)
                                                                                                                                      acctGAACGCGGTGATtgtacttcacccatattgttggacgatgaaatggaatagccactaatatatcaaagccaaaatacccttttgaaaagcgccttgaattcgagaatcactacttcacaacagacgatggttacaggatcatctcg               >  1:506408‑M1/2‑16 (MQ=255)
                                                                                                                                        gggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAACCCGGCGTTACGGGAAAAACTATTTACCTTTGAGGTGGATATTCTGCGCAAc             >  2:625765‑M2/12‑150 (MQ=255)
                                                                                                                                         ggtcagtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAACCCGGCGTTACGGGAAAAACTATTTACCTTTGAGGTGGATATTCGGCGCAACg            >  1:508138‑M2/11‑150 (MQ=255)
                                                                                                                                             agtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAACCCGGCGTTACGGGAAAAACTATTTACCTTTGAGGTGGATATTCTGCGCAACGATGt        >  2:630376‑M2/7‑150 (MQ=255)
                                                                                                                                              gtacaGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAACCCGGCGGTACGGGAAAAACTATTTACCTTTGAGGTGGATATTCTGCGCAACGATGTc       >  1:257246‑M2/6‑150 (MQ=255)
                                                                                                                                                   gATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtg  <  2:268296‑M1/150‑148 (MQ=255)
                                                                                                                                                   gATtgtactgcacccattttgttggacgatgaaatggaatagcccctaatatgtcaaagccaaaatacccttttgaaaagcgccttgaagtcgtgaatcactacttcacaactgatgatggttacaggatcatctcggcacgttttggtg  <  2:40441‑M1/150‑148 (MQ=255)
                                                                                                                                                   gATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAACCCGGCGTTACGGGAAAAACTATTTACCTTTGAGGTGGATATTCTGCGCAACGATGTCTGTGa  <  1:623934/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAATCCCTGCGCCACGCCCTCAATAAAGCGGTGCCTTATGTCCGCAATAACCCGGACAAACTGCATCTGTTTGTGGATAACGGTTCGCTGGTTGCCACGGGGGCCAGCTCCATGTCGTGGGAGTACCGTTACACCCTGAACGCGGTGATTGAGGATTTCAGCGGCGACCAGAATCTGCTGATGGCCCCGGTTTTGCTGTGGCTGCGGGATAACCAGCCCGATGCCATCAATAACCCGGCGTTACGGGAAAAACTATTTACCTTTGAGGTGGATATTCTGCGCAACGATGTCTGTGA  >  REL606/896316‑896612

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