Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 190,553 Δ14 bp coding (680‑693/726 nt) rpsB → 30S ribosomal subunit protein S2

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1905521 (0.020)42 (0.790) 34/196 0.8 97.7% coding (679/726 nt) rpsB 30S ribosomal subunit protein S2
?NC_000913 190567 = 1 (0.020)coding (694/726 nt) rpsB 30S ribosomal subunit protein S2

GTTATCCCGGGTAACGACGACGCAATCCGTGCTGTGACCCTGTACCTGGGCGCTGTTGCTGCAACCGTACGTGAAGGCCGTTCTCAGGATCTGGCTTCCCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAGGCGGCAAGCTCGAGAATTCCCGGG  >  NC_000913/190468‑190666
                                                                                                   |                                                                                                   
acggtgttgacttcgttatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctcggcgctgttgctgcaaccgtacgttaaggccgttctcc                                                                                                     <  2:806426‑M2/1‑1 (MQ=255)
    tgttgacttcgttatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGc                                                                                                 <  2:1269352‑M2/5‑1 (MQ=255)
       tgacttcgttatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGga                                                                                              >  1:709133‑M2/93‑100 (MQ=255)
       cgtcgtcgttacccgggataacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGga                                                                                              >  2:148615‑M2/93‑100 (MQ=255)
          cttcgttatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGgaaga                                                                                           >  2:945169‑M2/90‑100 (MQ=255)
          cttcgttatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGgaaga                                                                                           >  2:501084‑M2/90‑100 (MQ=255)
          cttcgttatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGgaaga                                                                                           >  1:236265‑M2/90‑100 (MQ=255)
            aaCGACGACGCAATCCGTGCTGTGACCCTGTACCTGGGCGCTGTTGCTGCAACCGTACGTGAAGGCCGTTCTCAGGATCTGGCTTCCCAGGCGGAAGaaa                                                                                         <  1:1175105/100‑1 (MQ=255)
             cgttatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAg                                                                                        <  1:27113‑M2/14‑1 (MQ=255)
             cgttatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAg                                                                                        >  2:807860‑M2/87‑100 (MQ=255)
              gttatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGc                                                                                       >  1:594306‑M2/86‑100 (MQ=255)
                tatcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGTCGGAAGAAagcct                                                                                     <  1:884536‑M2/17‑3 (MQ=255)
                 atcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTc                                                                                    <  1:443021‑M2/18‑1 (MQ=255)
                 atcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTc                                                                                    >  2:492946‑M2/83‑100 (MQ=255)
                 atcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTc                                                                                    >  2:626191‑M2/83‑100 (MQ=255)
                  tcccgggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCg                                                                                   >  1:334702‑M2/82‑100 (MQ=255)
                      gggtaacgacgacgcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGa                                                                               <  2:96525‑M2/23‑1 (MQ=255)
                                   gcaatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGtaataa                                                                  >  2:1236550‑M2/65‑100 (MQ=255)
                                     aatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAgg                                                                <  1:54271‑M2/38‑1 (MQ=255)
                                     aatccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAgg                                                                >  1:278147‑M2/63‑100 (MQ=255)
                                      atccgtgctgtgaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGc                                                               >  1:1148718‑M2/62‑100 (MQ=255)
                                                 gaccctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTcc                                                    >  2:257546‑M2/51‑100 (MQ=255)
                                                     ctgtacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCa                                                >  1:146987‑M2/47‑100 (MQ=255)
                                                        tacctgggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCaaaa                                             <  1:346190‑M2/57‑1 (MQ=255)
                                                             gggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGATTAATAAGGCTTGATAACTCCCCCAAAATAGtt                                        <  1:1236561‑M2/62‑1 (MQ=255)
                                                              ggcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTc                                       <  2:146985‑M2/63‑1 (MQ=255)
                                                               gcgctgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCg                                      <  2:824942‑M2/64‑1 (MQ=255)
                                                                   tgttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGtt                                  >  2:101100‑M2/33‑100 (MQ=255)
                                                                    gttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTg                                 <  2:973104‑M2/69‑1 (MQ=255)
                                                                    gttgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTg                                 <  2:334699‑M2/69‑1 (MQ=255)
                                                                      tgctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCa                               >  1:953559‑M2/30‑100 (MQ=255)
                                                                        ctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGa                             >  2:63404‑M2/28‑100 (MQ=255)
                                                                        ctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGa                             >  1:14068‑M2/28‑100 (MQ=255)
                                                                        ctgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGa                             >  2:1005573‑M2/28‑100 (MQ=255)
                                                                         tgcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGaa                            >  1:1158366‑M2/27‑100 (MQ=255)
                                                                          gcaaccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGaaa                           >  2:599426‑M2/26‑100 (MQ=255)
                                                                             accgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTCCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAggc                        >  1:718042‑M2/23‑100 (MQ=255)
                                                                              ccgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAggcg                       >  2:1254640‑M2/22‑100 (MQ=255)
                                                                               cgtacgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAggcgg                      >  1:361547‑M2/21‑100 (MQ=255)
                                                                                  acgtgaaggccgttctcCAGGCGTAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAGGCGGCa                    <  2:1296195‑M2/82‑1 (MQ=255)
                                                                                   cgtgaaggccgttctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAGGCGGCAAg                  >  1:37459‑M2/17‑100 (MQ=255)
                                                                                               tctcCAGGCGGAAGAAAGCTTCGTATAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAGGCGGCAAGCTCGAGAATTcc      <  2:1061148‑M2/96‑1 (MQ=255)
                                                                                               tctcCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAGGCGGCAAGCTCGAGAATTcc      >  1:1146072‑M2/5‑100 (MQ=255)
                                                                                                  cCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAGGCGGCAAGCTCGAGAATTCCCgg   >  2:958443‑M2/2‑100 (MQ=255)
                                                                                                   cAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAGCTCCCCCAAAATAGTTCGAGTTGCAGAAAGGCGGCAAGCTCGAGAATTCCCggg  <  1:449349/100‑1 (MQ=255)
                                                                                                   cAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAGGCGGCAAGCTCGAGAATTCCCggg  <  2:963709/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GTTATCCCGGGTAACGACGACGCAATCCGTGCTGTGACCCTGTACCTGGGCGCTGTTGCTGCAACCGTACGTGAAGGCCGTTCTCAGGATCTGGCTTCCCAGGCGGAAGAAAGCTTCGTAGAAGCTGAGTAATAAGGCTTGATAACTCCCCCAAAATAGTTCGAGTTGCAGAAAGGCGGCAAGCTCGAGAATTCCCGGG  >  NC_000913/190468‑190666

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