Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 3,720,136 IS3 (–) +3 bp :: +TCA coding (88‑90/213 nt) cspA → cold shock protein CspA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 315229 =NA (NA)20 (0.380)
+TCA
19/480 1.1 100% noncoding (1/1255 nt) IS3 repeat region
?NC_000913 3720136 = 0 (0.000)coding (88/213 nt) cspA cold shock protein CspA
* ? NC_000913 = 2171428NA (NA)51 (0.970) 50/484 0.1 100% noncoding (1258/1258 nt) IS3 repeat region
?NC_000913 = 3720138 0 (0.000)coding (90/213 nt) cspA cold shock protein CspA

TATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGT  >  NC_000913/314999‑315465
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
aatagattttatcacgcgatcgattatttatttcctgaaaacaaataacaccatccccgccaaatggcagggatcttagattctgtacttttaagaagagattacaggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttatggatagcagagaagtgtacgcacactgaTCTTACCCAGCAATAGTGGa                                                                                                                                                                                                                           <  2:310430‑M1/20‑1 (MQ=255)
          attacgcgatcgattatttatttcctgaaaacaaatcaaaaaatccccgccaaatggcagggatcttagattctgtgcttttaagcagcgattacaggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcttcccgagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTa                                                                                                                                                                                                                 <  1:40666‑M1/30‑1 (MQ=255)
            tacgcgatcgattatttatttcctgaaaacaaataaaaaaatccccgccaaatggcagggatcttagattctgtgcttttaagcagagattacaggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAg                                                                                                                                                                                                               >  1:272189‑M1/219‑250 (MQ=255)
            tacgcgatcgattatttatttcctgaaaacaaataaaaaaatccccgccaaatggcagggatcttagattctgtgcttttaagcagagattacaggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGt                                                                                                                                                                                                              <  1:196329‑M1/33‑1 (MQ=255)
                  atcgattatttatttcctgaaaacaaataaaaaaatccccgccaaatggcagggatcttagattctgtgcttttaagcagagattacaggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTa                                                                                                                                                                                                         <  1:108317‑M1/38‑1 (MQ=255)
                                  ctgaaaacaaataaaaaaatccccgccaaatggcagggatcttagattctgtgcttttaagcagagattacaggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGt                                                                                                                                                                                        <  2:125743‑M1/55‑1 (MQ=255)
                                                  aaatccccgccaaatggcagggatcttagattctgtgcttttaagcagagattacaggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaa                                                                                                                                                                        >  1:125743‑M1/181‑251 (MQ=255)
                                                            caaatggcagggatcttagattctgtgcttttaagcagagattacaggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACCCTCAGTCAGAGGTGAATCACATGACCAAAACAGTATCa                                                                                                                                                               <  2:397454‑M1/80‑1 (MQ=255)
                                                                 ggcagggatcttagagtctgtgcgtgtaagcagagattacaggctggttacgttaccagctgccgggccgttagcgccgctgtcgatggtgaaggacactttctgaccttcgtccagagagttgtaaccatcgttctggatagcagagaagtggacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCGAAGGGAGTAAACACTCAGTCAGAGGTGAATCACATGACAAAAACAGTATCAACCAGt                                                                                                                                                         >  2:40666‑M1/166‑251 (MQ=255)
                                                                            tagattctgtgcttttaagcagagattacaggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacacttgctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAGCCAGTAAAAAAc                                                                                                                                                  >  1:397454‑M1/155‑247 (MQ=255)
                                                                                                        caggctggttacgttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGt                                                                                                                  <  1:356665‑M1/125‑1 (MQ=255)
                                                                                                                    gttaccagctgccgggcctttagcgccgctttcgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGa                                                                                                        <  1:277324‑M1/135‑1 (MQ=255)
                                                                                                                                  ggcctttagcgccgctttcgctggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCAt                                                                                        <  1:404636‑M1/151‑1 (MQ=255)
                                                                                                                                                    cgatggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACATTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGc                                                                      <  2:396378‑M1/169‑1 (MQ=255)
                                                                                                                                                    cgatggtgaaggacactttctgaccttcgtccagagatttgtaaacaccgttctggatagcagacaagtgtacgaacactgaTCTTACCCAGCAATAGTGGCCACACGGCTAAGTGAGTGAACTCTCAGTCAGCGGTGACTCACATGACAACAACAGTATCAACCAGAAAAAAACCCCGTACACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGc                                                                      <  1:287866‑M1/169‑1 (MQ=255)
                                                                                                                                                      atggtgaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGccc                                                                    <  1:251287‑M1/171‑1 (MQ=255)
                                                                                                                                                           gaaggacactttctgaccttcgtccagagatttgtaaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGa                                                                <  1:318544‑M1/175‑1 (MQ=255)
                                                                                                                                                                                              aaccatcgttctggatagcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGgcg                            <  1:335172‑M1/211‑1 (MQ=255)
                                                                                                                                                                                                              agcagagaagtgtacgaacactgaTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAAt            <  2:349575‑M1/227‑1 (MQ=255)
                                                                                                                                                                                                                         gtacgaacactgaTCTTACCCAGCAATAGTTGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTAGCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAAGCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGt  >  1:485681‑M1/14‑250 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
TATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGT  >  NC_000913/314999‑315465

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