Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 3,560,930 IS5 (+) +4 bp coding (520‑523/831 nt) glpG ← rhomboid intramembrane serine protease

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3560930 =0 (0.000)31 (1.120) 28/266 0.2 100% coding (523/831 nt) glpG rhomboid intramembrane serine protease
?NC_000913 3652036 = NA (NA)noncoding (1195/1195 nt) IS5 repeat region
* ? NC_000913 = 35609330 (0.000)24 (0.860) 24/266 0.3 100% coding (520/831 nt) glpG rhomboid intramembrane serine protease
?NC_000913 = 3653230 NA (NA)noncoding (1/1195 nt) IS5 repeat region

AATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAGTGTTGGATCGAATGG  >  NC_000913/3560784‑3561077
                                                                                                                                                  |                                                                                                                                                   
aaTGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAAtt                                                                                                                                                    <  2:731/148‑1 (MQ=255)
 aactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTa                                                                                                                                                   <  1:54861‑M1/3‑1 (MQ=255)
       cTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaagg                                                                                                                                            >  5:119482‑M1/1‑143 (MQ=255)
            acgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTAc                                                                                                                                       <  7:80677‑M1/15‑1 (MQ=255)
            acgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTAc                                                                                                                                       >  7:77633‑M1/135‑149 (MQ=255)
                cttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAg                                                                                                                                   >  3:92335‑M1/131‑149 (MQ=255)
                 ttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAgg                                                                                                                                  <  6:46399‑M1/20‑1 (MQ=255)
                   cgTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtc                                                                                                                                <  6:119482‑M1/149‑19 (MQ=255)
                          cACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgata                                                                                                                         >  5:134114‑M1/1‑124 (MQ=255)
                          cACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGAAAATTAGggaaggtgcgaacaagtccctgata                                                                                                                         <  1:108313‑M1/149‑26 (MQ=255)
                                acggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCcaccgca                                                                                                                   >  5:46399‑M1/115‑149 (MQ=255)
                                  cAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgaga                                                                                                                    <  1:58902‑M1/146‑31 (MQ=255)
                                  cAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgaga                                                                                                                    >  2:58902‑M1/1‑116 (MQ=255)
                                        aataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCAccgcc                                                                                                              <  6:5726‑M1/40‑1 (MQ=255)
                                        aataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCAccgcc                                                                                                              <  4:55820‑M1/40‑1 (MQ=255)
                                        aataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCga                                                                                                            >  4:76662‑M1/107‑148 (MQ=255)
                                             gctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGtttt                                                                                                                            <  1:76370‑M1/26‑1 (MQ=255)
                                             gctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGtttt                                                                                                                            >  2:76370‑M1/102‑127 (MQ=255)
                                                  ttcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATaccaccac                                                                                                 <  1:126092‑M1/53‑1 (MQ=255)
                                                          atttacgggagaaaaaatcggctcaaacaggaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAg                                                                                         >  6:108442‑M1/89‑149 (MQ=255)
                                                           tttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagtcgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAgg                                                                                        <  2:25459‑M1/62‑1 (MQ=255)
                                                              cGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaccagg                                                                                     <  6:134114‑M1/149‑62 (MQ=255)
                                                                ccAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                                                                                   >  7:85950‑M1/1‑86 (MQ=255)
                                                                   agaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAagag                                                                                >  8:80677‑M1/80‑149 (MQ=255)
                                                                    gaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAaga                                                                                 >  6:112717‑M1/79‑147 (MQ=255)
                                                                    gaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAaga                                                                                 <  5:112717‑M1/69‑1 (MQ=255)
                                                                     aaaaaataggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAag                                                                                  >  1:29677‑M1/78‑145 (MQ=255)
                                                                     aaaaaataggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAag                                                                                  <  2:29677‑M1/68‑1 (MQ=255)
                                                                         ccGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                                                                                >  8:49201‑M1/1‑77 (MQ=255)
                                                                              aaaCCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGtgaagttgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatataacagggttcatcatgagtcat                                                                     <  7:33246‑M1/149‑78 (MQ=255)
                                                                                 ccACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaa                                                                  <  8:85950‑M1/149‑81 (MQ=255)
                                                                                  aaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCa                                                                                          >  4:95285‑M1/65‑124 (MQ=255)
                                                                                  aaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCa                                                                                          <  3:95285‑M1/60‑1 (MQ=255)
                                                                                    acatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCgaga                                                               >  2:126092‑M1/63‑149 (MQ=255)
                                                                                         aagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGc                                                          >  3:55820‑M1/58‑149 (MQ=255)
                                                                                          agaaatgaaatgactgagtcagtcgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCa                                                         >  1:25459‑M1/57‑149 (MQ=255)
                                                                                             aatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTa                                                      >  5:5726‑M1/54‑149 (MQ=255)
                                                                                              aaTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatca                                                                             >  7:27518‑M1/1‑56 (MQ=255)
                                                                                              aaTTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatca                                                                             <  8:27518‑M1/125‑70 (MQ=255)
                                                                                                       tgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgag                                                                         <  2:80792‑M1/120‑74 (MQ=255)
                                                                                                       tgcACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgag                                                                         >  1:80792‑M1/1‑47 (MQ=255)
                                                                                                                  gCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagagcagggttcatcatgagtcatcaacataccttcgccgacagtgaattcagcactaag                                 >  7:37488‑M1/1‑36 (MQ=255)
                                                                                                                   agaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGt                                                                                       <  7:110226‑M1/63‑1 (MQ=255)
                                                                                                                   agaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGt                                                                                       >  8:110226‑M1/32‑94 (MQ=255)
                                                                                                                   agaagaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAaga                                                                                 >  2:6492‑M1/32‑100 (MQ=255)
                                                                                                                    ttAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccataga                                                                                          >  3:48238‑M1/1‑34 (MQ=255)
                                                                                                                    ttAACAGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccataga                                                                                          <  4:48238‑M1/90‑57 (MQ=255)
                                                                                                                       gaatttccccgcttattcgcaccttccTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAAc                            <  7:30846‑M1/122‑1 (MQ=255)
                                                                                                                       aCAGGGCGTTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccg                            >  8:107015‑M1/1‑31 (MQ=255)
                                                                                                                         aGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcag                                        >  7:45874‑M1/1‑29 (MQ=255)
                                                                                                                         aGGGCGCTGATAAGCGTAATGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcag                                        <  8:45874‑M1/135‑107 (MQ=255)
                                                                                                                                            tGACAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttc                                                                                 <  4:101027‑M1/75‑66 (MQ=255)
                                                                                                                                              aCAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagca                                      <  1:65092‑M1/116‑109 (MQ=255)
                                                                                                                                              aCAATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagca                                      >  2:65092‑M1/1‑8 (MQ=255)
                                                                                                                                                 cTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAGTGTTGGATCGAATgg  <  3:7475‑M1/148‑1 (MQ=255)
                                                                                                                                                 cTTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAGTGTTGGATCGAATg   >  8:30846‑M1/2‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AATGCCACTTTGCGGATCGCGTTCGCCACGTAGCCAGACGTAGCCCATCAGCGCATACACCACGCCAGAAAGCCCGCCAAACCACGGCCCGCTGAATTTTTGCTGCACATAGCCGCTTAACAGGGCGCTGATAAGCGTAATGACAATTAGCTTACCGCTACCGAGGCGTTTTTCCACCGCACCGCCGAGATACCACCACCAGAGCAGGTTAAAGAGGATATGCATCAGCGAGAAGTGCATTAACGCGTGGGTGAAGTAACGCCAGAACTCAAATTTCAGTGTTGGATCGAATGG  >  NC_000913/3560784‑3561077

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