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

GTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACGCATTACGGATTCCGTTTGAATT  >  REL606/2322196‑2322494
                                                                                                                                                     |                                                                                                                                                     
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:283933‑M1/150‑79 (MQ=255)
                                                                              ttACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataac                                                                         <  2:18481‑M1/150‑79 (MQ=255)
                                                                              ttACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataac                                                                         <  2:143218‑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:228220‑M1/150‑124 (MQ=255)
                                                                                                                           gactgaAGCCCGGCAGTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggat                            <  2:96219‑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:636998‑M1/1‑12 (MQ=255)
                                                                                                                                          gTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaat             >  2:451002‑M1/1‑12 (MQ=255)
                                                                                                                                          gTGGGGAGATCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaat             >  2:73962‑M1/1‑12 (MQ=255)
                                                                                                                                              ggAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACGCATTACGGATTCCGt         <  2:126262‑M1/148‑1 (MQ=255)
                                                                                                                                                 gaTCTcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaatgccgggg      >  2:559105‑M1/1‑5 (MQ=255)
                                                                                                                                                    ctcCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCCAACATACAGCCCGTCGCTGGTCTTTAACCGTCTCTCGCGCAGCACCACCGCCCCGGCCCGTGGGTACTTCCACCGGTGTACCTTCGCTCCACCCATTACGGAtttctctttgat   >  1:555834/1‑140 (MQ=255)
                                                                                                                                                     tcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgataactggtcagcaattctggcccatattggtaagcccgaagaactggatacttcggcacgtaatgccggggctct  >  2:671965‑M1/1‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACGCATTACGGATTCCGTTTGAATT  >  REL606/2322196‑2322494

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.