Predicted mutation
evidence seq id position mutation annotation gene description
JC REL606 3,145,826 Δ7 bp intergenic (+157/‑32) dnaG → / → rpoD DNA primase/RNA polymerase sigma factor

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 31458251 (0.020)47 (0.820) 33/146 0.6 97.9% intergenic (+156/‑39) dnaG/rpoD DNA primase/RNA polymerase sigma factor
?REL606 3145833 = 1 (0.020)intergenic (+164/‑31) dnaG/rpoD DNA primase/RNA polymerase sigma factor

ATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGT  >  REL606/3145758‑3145908
                                                                           |                                                                           
cctcgtaattatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgcca                                                                             <  1:1713472‑M2/1‑1 (MQ=255)
 ctcgtaattatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccaa                                                                            >  1:1558450‑M2/75‑76 (MQ=255)
  tcgtaattatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAAc                                                                           <  1:1792157‑M2/3‑1 (MQ=255)
   cgtaattatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAAcc                                                                          >  2:482475‑M2/73‑76 (MQ=255)
   cgtaattatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAAcc                                                                          >  1:1365234‑M2/73‑76 (MQ=255)
     taattatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTc                                                                        >  2:979081‑M2/71‑76 (MQ=255)
     taattatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTc                                                                        <  2:1625892‑M2/6‑1 (MQ=255)
       attatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCAt                                                                      <  1:1176424‑M2/8‑1 (MQ=255)
       attatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCAt                                                                      <  1:789626‑M2/8‑1 (MQ=255)
         tatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGa                                                                    <  1:1062386‑M2/10‑1 (MQ=255)
         tatcgttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGa                                                                    <  2:79610‑M2/10‑1 (MQ=255)
             gttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATa                                                                >  1:688554‑M2/63‑76 (MQ=255)
              ttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATaa                                                               >  1:1789767‑M2/62‑76 (MQ=255)
              ttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATaa                                                               >  1:1355056‑M2/62‑76 (MQ=255)
              ttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATaa                                                               >  1:1402415‑M2/62‑76 (MQ=255)
              ttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATaa                                                               >  1:431907‑M2/62‑76 (MQ=255)
              ttggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATaa                                                               >  1:231341‑M2/62‑76 (MQ=255)
               tggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAg                                                              >  1:1047948‑M2/61‑76 (MQ=255)
                ggcggtaaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAgt                                                             >  2:1860364‑M2/60‑76 (MQ=255)
                     taaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGa                                                        >  2:1136264‑M2/55‑76 (MQ=255)
                      aaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGAt                                                       <  1:1550254‑M2/23‑1 (MQ=255)
                      aaacaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGAt                                                       <  1:716908‑M2/23‑1 (MQ=255)
                        acaatcgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATAc                                                     <  1:1250920‑M2/25‑1 (MQ=255)
                        acaaccgttggatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATAc                                                     <  1:826267‑M2/25‑1 (MQ=255)
                                cAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATAAGTGTGGATACCGTCTTAt                                             >  1:246586/1‑76 (MQ=255)
                                  gatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATgg                                           <  2:293657‑M2/35‑1 (MQ=255)
                                  gatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATgg                                           <  2:1722482‑M2/35‑1 (MQ=255)
                                  gatttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATgg                                           <  2:934466‑M2/35‑1 (MQ=255)
                                     ttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGc                                        >  2:355126‑M2/39‑76 (MQ=255)
                                     ttcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGc                                        >  1:1209251‑M2/39‑76 (MQ=255)
                                      tcagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCa                                       <  1:444128‑M2/39‑1 (MQ=255)
                                       cagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCaa                                      <  1:1804670‑M2/40‑1 (MQ=255)
                                       cagcgttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCaa                                      >  1:157688‑M2/37‑76 (MQ=255)
                                           gttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAAcc                                  <  2:404592‑M2/44‑1 (MQ=255)
                                            ttaacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAAccc                                 <  2:1463458‑M2/45‑1 (MQ=255)
                                             taacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCg                                <  2:654219‑M2/46‑1 (MQ=255)
                                             taacggctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCg                                >  2:1553003‑M2/31‑76 (MQ=255)
                                              aacggctgaaggacatcgggtcaatcgccGACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGc                               <  1:968768‑M2/47‑1 (MQ=255)
                                                   ctgaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTca                          >  1:435766‑M2/25‑76 (MQ=255)
                                                     gaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTcaca                        >  2:1290783‑M2/23‑76 (MQ=255)
                                                     gaaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTcaca                        <  1:486354‑M2/54‑1 (MQ=255)
                                                      aaggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAg                       >  2:1439141‑M2/22‑76 (MQ=255)
                                                        ggacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCt                     <  2:1077118‑M2/57‑1 (MQ=255)
                                                         gacatcgggtcaatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTg                    >  1:579954‑M2/19‑76 (MQ=255)
                                                                   caatcgccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTg          >  1:1805150‑M2/9‑76 (MQ=255)
                                                                        gccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCAcc     >  2:1666430‑M2/4‑76 (MQ=255)
                                                                        gccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCAcc     >  2:1353058‑M2/4‑76 (MQ=255)
                                                                        gccAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCAcc     >  2:724389‑M2/4‑76 (MQ=255)
                                                                           aaCCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGt  <  1:1001828/76‑1 (MQ=255)
                                                                           |                                                                           
ATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATGAAATAAGTGTGGATACCGTCTTATGGAGCAAAACCCGCAGTCACAGCTGAAACTTCTTGTCACCCGT  >  REL606/3145758‑3145908

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