Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 2,322,340 IS186 (+) +6 bp 66.0% coding (131‑136/963 nt) menC ← O‑succinylbenzoate synthase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2322340 =9 (0.360)7 (0.270) 5/272 NT 42.6% coding (136/963 nt) menC O‑succinylbenzoate synthase
?REL606 = 2774196 NA (NA)noncoding (1343/1343 nt) IS186 repeat region
* ? REL606 = 23223459 (0.360)28 (1.070) 19/272 NT 74.8% coding (131/963 nt) menC O‑succinylbenzoate synthase
?REL606 2772854 = NA (NA)noncoding (1/1343 nt) IS186 repeat region

TCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACGCATTACGGATTCCGT  >  REL606/2322193‑2322487
                                                                                                                                                   |                                                                                                                                                   
tCAGTTTGGCCAATGAACAGCTTACGCAAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGgaga                                                                                                                                                   <  1:169420/150‑1 (MQ=255)
tCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGgaga                                                                                                                                                   <  2:511310/150‑1 (MQ=255)
tCAGTTCTGCCAATGAACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCCGCCGGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCGCTTCCCAGGGTTCCTGACTGGAGCCCGGCAGTGGGgaga                                                                                                                                                   <  1:237804/150‑1 (MQ=255)
  aGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGAtc                                                                                                                                                 >  1:446785/1‑150 (MQ=255)
   tctggagcttacaggtggccattcgtgggacagtatccctgacagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaaggaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgtttatgggAGAtct                                                                                                                                                <  2:593519‑M1/6‑1 (MQ=255)
        gCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccat                                                                                                                                           >  2:132041‑M1/1‑145 (MQ=255)
             acaggtggctattcgtgggacagtatccctgacagcctacaaaacgcaattgaagaacgtgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggAGATCTCCCCCCAccc                                                                                                                                      >  1:591779‑M1/135‑150 (MQ=255)
             acaggtggccattcgtgggacagtatccctgacagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggAGATCTCCCCCCAccc                                                                                                                                      >  2:652303‑M1/135‑150 (MQ=255)
             acaggtggccattcgtgggacagtagccctgacagcctacaaaacgcaattgaagaacgcgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggAGATCTCCCCCCAccc                                                                                                                                      >  2:680919‑M1/135‑150 (MQ=255)
                   gCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCGCACTTTGCGCCTCTTCCAAGGTTTCCTGACTGAAGCCCCGCAGTGAGGAGATCTtccataagcgctaact                                                                                                                                >  1:74811‑M1/1‑134 (MQ=255)
                   gCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTcgcg                                                                                                                                >  2:299995/1‑150 (MQ=255)
                        cGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagg                                                                                                                           >  2:73707‑M1/1‑129 (MQ=255)
                          ccAAAGGCGACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGACGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggt                                                                                                                         <  1:463202‑M1/150‑24 (MQ=255)
                          ccAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggt                                                                                                                         <  2:5503‑M1/150‑24 (MQ=255)
                               ggCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtgg                                                                                                                    >  2:603254‑M1/1‑122 (MQ=255)
                                 ccACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggta                                                                                                                  <  2:239126‑M1/150‑31 (MQ=255)
                                    cGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtatta                                                                                                               >  1:227572‑M1/1‑117 (MQ=255)
                                        aGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcc                                                                                                           <  2:154745‑M1/150‑38 (MQ=255)
                                        aGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCAAGTTATTTACCCAGGCAAGCAGCACACCTGGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgccaacttaagggttgtggtattacgcc                                                                                                           <  1:419764‑M1/150‑38 (MQ=255)
                                          gCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACAt                                                                                                         <  2:448054/150‑1 (MQ=255)
                                            aTCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATAc                                                                                                       <  1:119436/150‑1 (MQ=255)
                                                      gCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                                                                                             <  1:502458‑M1/150‑52 (MQ=255)
                                                       cTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacg                                                                                            <  2:437243‑M1/150‑53 (MQ=255)
                                                                 ccTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGACGCGGGTTTTTAAc                                                                                  <  1:126962/150‑1 (MQ=255)
                                                                            gTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTcgcg                                                                       >  2:63354/1‑150 (MQ=255)
                                                                              tATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCa                                                                     >  2:19509/1‑150 (MQ=255)
                                                                                 ttACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataac                                                                  <  2:143218‑M1/150‑79 (MQ=255)
                                                                                 ttACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataac                                                                  <  2:283933‑M1/150‑79 (MQ=255)
                                                                                 ttACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataac                                                                  <  2:18481‑M1/150‑79 (MQ=255)
                                                                                  tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAac                                                                 >  2:514545‑M1/66‑150 (MQ=255)
                                                                                       aGGCAAGCAGCACACTTTTCCCCTCTCCCCATGTTTCCTGACTGAAGCCTGGCCGGTGGGATATCTCCCCCCCCCCTTCGCCCTCGCCTTCACGCAGGCAAACATACACCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCAccc                                                            <  1:489782/150‑1 (MQ=255)
                                                                                           ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCcgc                                                        <  2:243916‑M1/94‑1 (MQ=255)
                                                                                             gcagcaCACTTTGCGCCTCTTCCAAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTaccataagcgctcacttaagggttgtggtattacgcctgctatgatttaacgtgccgatgaattactctcacgataactggtcagcaatt                                                      >  1:138351‑M1/1‑60 (MQ=255)
                                                                                                       ttGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccata                                            >  2:44107‑M1/1‑50 (MQ=255)
                                                                                                             ctctTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccatacgcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggta                                      <  1:282404‑M1/150‑107 (MQ=255)
                                                                                                             ctctTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctcacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggta                                      <  1:73491‑M1/150‑107 (MQ=255)
                                                                                                             ctctTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggta                                      <  2:218419‑M1/150‑107 (MQ=255)
                                                                                                                 tCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcc                                  >  2:202431‑M1/1‑40 (MQ=255)
                                                                                                                       gTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATAcc                            <  1:299146/150‑1 (MQ=255)
                                                                                                                         ttCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTg                          <  2:159894/150‑1 (MQ=255)
                                                                                                                            acccttaagttagcgcttatgggAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTgcgc                       <  2:364082‑M1/127‑1 (MQ=255)
                                                                                                                              gactgaAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggat                     <  2:96219‑M1/150‑124 (MQ=255)
                                                                                                                              gactgaAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggat                     <  2:228220‑M1/150‑124 (MQ=255)
                                                                                                                                   aaGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttc                >  2:362637‑M1/1‑22 (MQ=255)
                                                                                                                                     gCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcgg              >  2:56586‑M1/1‑20 (MQ=255)
                                                                                                                                            aGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgctttaacgtgccgatgaattactcgcactataactggtcagcaattctggcccatactggtaagcccgaagaactggatactgcggcacgtaa       >  1:462827‑M1/1‑13 (MQ=255)
                                                                                                                                             gTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaat      >  2:451002‑M1/1‑12 (MQ=255)
                                                                                                                                             gTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaat      >  2:636998‑M1/1‑12 (MQ=255)
                                                                                                                                             gTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaat      >  2:73962‑M1/1‑12 (MQ=255)
                                                                                                                                                 ggAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACGCATTACGGATTCCGt  <  2:126262‑M1/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACGCATTACGGATTCCGT  >  REL606/2322193‑2322487

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.