Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 944,912 Δ1 bp :: IS186 (+) +6 bp :: Δ1 bp intergenic (+16/+14) dmsC → / ← ycaC dimethyl sulfoxide reductase, anaerobic, subunit C/putative isochorismatase family hydrolase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 608008 =NA (NA)33 (0.850) 33/292 0.4 100% noncoding (2/1345 nt) IS186 repeat region
?NC_000913 = 944917 0 (0.000)intergenic (+21/+14) dmsC/ycaC dimethyl sulfoxide reductase, anaerobic, subunit C/putative isochorismatase family hydrolase
* ? NC_000913 = 609350NA (NA)40 (1.030) 37/292 0.2 100% noncoding (1344/1345 nt) IS186 repeat region
?NC_000913 944912 = 0 (0.000)intergenic (+16/+19) dmsC/ycaC dimethyl sulfoxide reductase, anaerobic, subunit C/putative isochorismatase family hydrolase

CATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTAT  >  NC_000913/609205‑609500
                                                                                                                                                 |                                                                                                                                                      
cATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTa                                                                                                                                                    >  2:103016/1‑150 (MQ=31)
                   gCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgt                                                                                                                                <  1:602384/151‑1 (MQ=255)
                    ccATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtc                                                                                                                               <  1:449119/151‑1 (MQ=255)
                           tGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggtc                                                                                                                          >  1:432605/1‑149 (MQ=255)
                                          gACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCt                                                                                                         <  2:612111/151‑1 (MQ=255)
                                              gCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCtt                                                                                                     <  1:372464/151‑1 (MQ=255)
                                               ccTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTaa                                                                                                                <  1:270902/139‑1 (MQ=255)
                                               ccTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTaa                                                                                                                >  2:270902/1‑139 (MQ=255)
                                                 tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTac                                                                                                  >  1:200362/1‑151 (MQ=255)
                                                      aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATACCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACttt                                                                                             >  2:579321/1‑151 (MQ=255)
                                                       aCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTc                                                                                            <  1:504099/151‑1 (MQ=255)
                                                           aaTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGccc                                                                                        <  2:32542/151‑1 (MQ=255)
                                                              tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGt                                                                                      <  1:587048/150‑1 (MQ=255)
                                                                aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGt                                                                                   >  2:463568/1‑151 (MQ=255)
                                                                                    tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCg                                                                 <  1:139642/149‑1 (MQ=255)
                                                                                        gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGt                                                           <  1:478080/151‑1 (MQ=255)
                                                                                         aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTa                                                          <  1:410724/151‑1 (MQ=255)
                                                                                             aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGGCAGTGATGGCTGCGGGCGAATCGt                                                           >  2:212106/1‑146 (MQ=255)
                                                                                              ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTAcc                                                        >  2:623084/1‑148 (MQ=255)
                                                                                               cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAt                                                    <  1:70772/151‑1 (MQ=255)
                                                                                                 aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAtgt                                                  <  2:100256/151‑1 (MQ=255)
                                                                                                      tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTc                                              >  2:370640/1‑150 (MQ=255)
                                                                                                          aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAAtt                                          >  1:248423/1‑150 (MQ=255)
                                                                                                           ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTa                                         >  1:294185/1‑150 (MQ=255)
                                                                                                                  gATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGtt                                 <  1:581016/151‑1 (MQ=255)
                                                                                                                  gATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGtt                                 <  2:480984/151‑1 (MQ=255)
                                                                                                                         cAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGcac                            >  2:535911/1‑149 (MQ=255)
                                                                                                                                tAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATg                   >  2:227945/1‑151 (MQ=255)
                                                                                                                                  aGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGAc                 >  1:281754/1‑151 (MQ=255)
                                                                                                                                     tAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTa               <  1:326141/150‑1 (MQ=255)
                                                                                                                                                 gAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTGCTTTTAt  <  1:500219/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
CATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTAT  >  NC_000913/609205‑609500

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