Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 502,543 IS150 (–) +3 bp 100% coding (848‑850/1095 nt) ybbB ← tRNA 2‑selenouridine synthase, selenophosphate‑dependent

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 502543 =1 (0.040)12 (0.450) 11/278 NT 92.2% coding (850/1095 nt) ybbB tRNA 2‑selenouridine synthase, selenophosphate‑dependent
?REL606 588495 = NA (NA)noncoding (1/1977 nt) IS150 repeat region
* ? REL606 = 5025451 (0.040)27 (1.010) 21/278 NT 96.4% coding (848/1095 nt) ybbB tRNA 2‑selenouridine synthase, selenophosphate‑dependent
?REL606 = 666130 NA (NA)noncoding (1443/1443 nt) IS150 repeat region

AATATTCTTCAAGTAACGGCACCAGCCAGGCCAGATGACCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTC  >  REL606/502396‑502691
                                                                                                                                                     |                                                                                                                                                  
aaTATTCTTCAAGTAACGGCACCAGCCAGGCCAGATGACCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGAt                                                                                                                                                    <  2:520485/150‑1 (MQ=255)
 accaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATg                                                                                                                                                   >  2:631600‑M1/147‑150 (MQ=255)
 aaccaaaacgtgccgatatgatcctgtaaccatcatcagttgtgaagagtgattcacgactttaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATg                                                                                                                                                   >  1:99890‑M1/147‑150 (MQ=255)
     tcttcAAGTAACGGCACCAGCCAGGCCAGATGACCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtac                                                                                                                                               <  2:138071‑M1/150‑6 (MQ=255)
            gTAACGGCACCAGCCAGGCCAGATGACCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATGCAGGTATTCGc                                                                                                                                        <  2:608094/150‑1 (MQ=255)
            ccgagatgatcctgtaaccataatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggcgttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGc                                                                                                                                        >  1:104843‑M1/136‑150 (MQ=255)
                 atgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGCAATAc                                                                                                                                   >  2:660007‑M1/131‑150 (MQ=255)
                             ccatcatctgttgtgaagtagtgactcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGCAATACTCCTGCCAGccc                                                                                                                       >  1:465578‑M1/119‑150 (MQ=255)
                                 gATGACCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaaca                                                                                                                   <  1:418149‑M1/150‑34 (MQ=255)
                                  ctgACCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCCGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacac                                                                                                                  <  1:251567‑M1/149‑35 (MQ=255)
                                  atcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTtc                                                                                                                  >  2:385122‑M1/114‑150 (MQ=255)
                                  atcagttgtgaagtagtgattcacgacttcaacgcgcttttcaaaagtggattttggctttgacagattaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTtc                                                                                                                  >  1:456320‑M1/114‑150 (MQ=255)
                                  aTGACCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacac                                                                                                                  <  2:107394‑M1/150‑35 (MQ=255)
                                  aTGACCGTCGGTGCTGCCGGTGGGGAGTTGCGTTGTCAGTGCTGCATCCAGTCTTGCAGCCATTTTGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaagagttggacagttaaacac                                                                                                                  <  1:179027‑M1/150‑35 (MQ=255)
                                    gACCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacga                                                                                                                >  2:176924‑M1/1‑114 (MQ=255)
                                     aCCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgag                                                                                                               >  2:485904‑M1/1‑113 (MQ=255)
                                                tgcCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtc                                                                                                    <  2:15447‑M1/150‑49 (MQ=255)
                                                gagacagggtgagttgcgTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtc                                                                                                    <  1:571013‑M1/143‑49 (MQ=255)
                                                     gtggtgAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgatt                                                                                               >  1:50670‑M1/1‑97 (MQ=255)
                                                       ggtgAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtcctgtctc                                                                                             >  2:572679‑M1/1‑95 (MQ=255)
                                                          gAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgat                                                                                          <  2:191013‑M1/150‑59 (MQ=255)
                                                                aaggcgcttttcaaaagggtattttggctttggcatattaggggatattccatttcttcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGAAATGCTCCAGGCAGCCCTGTTCGTCGCCATACTCGTGGGTAAAGTCatgatg                                                                                    >  1:113225‑M1/84‑150 (MQ=255)
                                                                        gTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtc                                                                            >  2:519156‑M1/1‑78 (MQ=255)
                                                                        gTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGGAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagatggacagttaaacacgaggcatataggtctgattccgatattcaattggagtc                                                                            <  1:24433‑M1/150‑73 (MQ=255)
                                                                            tgcATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagac                                                                        <  2:26227‑M1/150‑77 (MQ=255)
                                                                            tgcATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagac                                                                        <  1:484527‑M1/150‑77 (MQ=255)
                                                                               aTCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagacctt                                                                     >  2:437341‑M1/1‑71 (MQ=255)
                                                                                  gtattttggcattgacatattaggggctattccatttcatcgtccaacaaaatgggtgaagtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATActc                                                                  >  1:204‑M1/66‑150 (MQ=255)
                                                                                   tattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATActct                                                                 >  2:642218‑M1/65‑150 (MQ=255)
                                                                                               gacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTCCGCGTGGGTAAAATCATGATGCATACGCAAGCAATACTCTTCATTCAGgcgc                                                     <  1:266904‑M1/98‑1 (MQ=255)
                                                                                                    attaggggctattccatttcatcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAg                                                <  2:673979‑M1/103‑1 (MQ=255)
                                                                                                            ctattccatttcatcgtccaacaaaatgggtgcagtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGAtct                                        <  1:444170‑M1/111‑1 (MQ=255)
                                                                                                             cTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgt                                       >  2:408201‑M1/1‑41 (MQ=255)
                                                                                                             cTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagctggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgt                                       >  2:134033‑M1/1‑41 (MQ=255)
                                                                                                               gTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtag                                     <  2:410407‑M1/150‑112 (MQ=255)
                                                                                                                 aGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagta                                   <  1:436101‑M1/150‑114 (MQ=255)
                                                                                                                  gCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtat                                  <  2:396743‑M1/150‑115 (MQ=255)
                                                                                                                       aGCCGACGCTTAATCGCCGAAAGTCCGTGATtgtactgcccccaaaaagttggacagttaaacccgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtattcaat                             <  1:109811‑M1/150‑120 (MQ=255)
                                                                                                                         ccGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacactttaacacgggccttatcggtctgattcctatattcaattggggtcagaccttttaatttcaggctcagtcttctgctgttgtagtattctatat                           <  1:185014‑M1/150‑122 (MQ=255)
                                                                                                                         ccGACGCTTAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagctggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtattcaatat                           <  1:133645‑M1/150‑122 (MQ=255)
                                                                                                                                ttAATCGCCGAAAGTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtattcaatatattccgt                    <  1:517670‑M1/150‑129 (MQ=255)
                                                                                                                                   agcgccgcaaGTCCGTGATtgtactgacccgaaccagttgagcagttaaacacgactcctataggtctgattccgttattcagttggagtcagaccttttaatttcagcccaattcgtctgctgttgtagtattcaatatattccgtaac                 <  1:325711‑M1/148‑132 (MQ=255)
                                                                                                                                         gtgcagtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATcgcc           <  2:49685‑M1/140‑1 (MQ=255)
                                                                                                                                             gTCCGTGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtattcaatatattccgtaacagcatccttc       <  2:331136‑M1/150‑142 (MQ=255)
                                                                                                                                              gtacagATGCCGGGATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTCCGCGGGGGTAACATCTTCATGCATACGCACGACATACCCTCCGTCCAGGCGCTCAAGCCGGATCTGAAGCGGATCTTCTACCCCCCCACTCGCCGACTg      >  1:525770/6‑150 (MQ=255)
                                                                                                                                              gtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTg      <  2:644447‑M1/145‑1 (MQ=255)
                                                                                                                                              gtacaGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTg      <  2:342740‑M1/145‑1 (MQ=255)
                                                                                                                                                 caGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGt   >  2:445432‑M1/3‑150 (MQ=255)
                                                                                                                                                  tGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtattcaatatattccgtaacagcatccttcctgtc  >  2:458239‑M1/1‑4 (MQ=255)
                                                                                                                                                  tGATtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtattcaatatattccgtaacagcatccttcagttc  >  1:240799‑M1/1‑4 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AATATTCTTCAAGTAACGGCACCAGCCAGGCCAGATGACCGTCGGTGCTGCCGGTGGTGAGTTGCGTTGTCAGTGCTGCATCCAGCCTTGCAGCCAGTTCGTTATAGCGCTGTAGCCCCAGCCGACGCTTAATCGCCGAAAGTCCGTGATGCAGGTATTCGCAATACTCCTGCCAGCCCTGTTCGTCGCCGTACGCGTGGGTAAAATCATGATGCATACGCAAGAAATACTCTTCGTTCAGGCGCTCAAGACGGATCTCAAACGGATCTTCTACCACCGCAATCGCCGCCTGAGTC  >  REL606/502396‑502691

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