Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 4,163,094 Δ3 bp coding (2239‑2241/4029 nt) rpoB → DNA‑directed RNA polymerase subunit beta

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 4163094 4163096 3 54 [0] [0] 55 rpoB DNA‑directed RNA polymerase subunit beta

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 41630930 (0.000)54 (1.010) 42/150 0.2 100% coding (2238/4029 nt) rpoB DNA‑directed RNA polymerase subunit beta
?REL606 4163097 = 0 (0.000)coding (2242/4029 nt) rpoB DNA‑directed RNA polymerase subunit beta

TGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTG  >  REL606/4163021‑4163168
                                                                        |                                                                           
tgtCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatc                                                                          <  2:950730‑M1/76‑4 (MQ=255)
 gtCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcg                                                                         >  2:1030267‑M1/1‑72 (MQ=255)
   cGTTCAGTACGTGGATGCTTCCCGTATCGCTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgac                                                                       <  2:292156‑M1/76‑7 (MQ=255)
     ttCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacat                                                                     >  1:89233‑M1/1‑68 (MQ=255)
      tCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatc                                                                    >  1:899270‑M1/1‑67 (MQ=255)
       cAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatct                                                                   <  1:429222‑M1/76‑11 (MQ=255)
       cAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGAAatcgacatct                                                                   >  2:190190‑M1/1‑66 (MQ=255)
        aGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatcta                                                                  >  1:404625‑M1/1‑65 (MQ=255)
              tGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacct                                                            >  2:1580900‑M1/1‑59 (MQ=255)
                 aTGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgac                                                         >  1:620815‑M1/1‑56 (MQ=255)
                 aTGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgac                                                         <  1:1060440‑M1/76‑21 (MQ=255)
                  tGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgacc                                                        <  1:1331593‑M1/76‑22 (MQ=255)
                         cGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaataca                                                 <  2:1548577‑M1/76‑29 (MQ=255)
                           tATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacc                                               >  2:792543‑M1/1‑46 (MQ=255)
                           tATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacc                                               >  2:616868‑M1/1‑46 (MQ=255)
                            aTCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacaccc                                              >  1:342797‑M1/1‑45 (MQ=255)
                            aTCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacaccc                                              >  2:1337698‑M1/1‑45 (MQ=255)
                              cGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgt                                            >  1:67174‑M1/1‑43 (MQ=255)
                               gTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgtt                                           >  1:1233402‑M1/1‑42 (MQ=255)
                               gTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgtt                                           <  2:899088‑M1/76‑35 (MQ=255)
                                ttATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttc                                          >  2:695701‑M1/1‑41 (MQ=255)
                                 tATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttct                                         >  1:1291712‑M1/1‑40 (MQ=255)
                                  aTCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttcta                                        >  2:1049343‑M1/1‑39 (MQ=255)
                                  aTCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttcta                                        <  1:12292‑M1/76‑38 (MQ=255)
                                     aaaGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaacc                                     >  2:1074510‑M1/1‑36 (MQ=255)
                                      aaGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaacca                                    >  2:1203931‑M1/1‑35 (MQ=255)
                                       aGTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccag                                   <  1:696159‑M1/76‑43 (MQ=255)
                                        gTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctggccaaatacacccgttctaaccaga                                  <  2:699627‑M1/76‑44 (MQ=255)
                                        gTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccaga                                  <  2:699616‑M1/76‑44 (MQ=255)
                                        gTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccaga                                  <  1:179742‑M1/76‑44 (MQ=255)
                                        gTTAACGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccaga                                  <  2:1650460‑M1/76‑44 (MQ=255)
                                          tAACGAAGACGAGATGTATCCGGGTGACGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaac                                >  2:1512761‑M1/1‑31 (MQ=255)
                                           aaCGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaaca                               <  2:1206142‑M1/76‑47 (MQ=255)
                                           aaCGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaaca                               >  1:1349772‑M1/1‑30 (MQ=255)
                                            aCGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacac                              <  2:1380216‑M1/76‑48 (MQ=255)
                                            aCGAAGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacac                              <  2:375990‑M1/76‑48 (MQ=255)
                                               aaGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctg                           <  2:801193‑M1/76‑51 (MQ=255)
                                               aaGACGAGATGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctg                           <  1:1545712‑M1/76‑51 (MQ=255)
                                                    gagaTGTATCCGGGTGAAGCGatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatca                      <  1:557088‑M1/76‑56 (MQ=255)
                                                    gagaTGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatca                      <  1:151737‑M1/76‑56 (MQ=255)
                                                     agaTGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaa                     <  1:1497081‑M1/76‑57 (MQ=255)
                                                       aTGTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaacc                   >  1:43480‑M1/1‑18 (MQ=255)
                                                         gTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccag                 >  1:495274‑M1/1‑16 (MQ=255)
                                                         gTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccag                 <  1:430011‑M1/76‑61 (MQ=255)
                                                         gTATCCGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccag                 <  2:1299743‑M1/76‑61 (MQ=255)
                                                             ccGGGTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccagatgc             <  2:1129927‑M1/76‑65 (MQ=255)
                                                               gggTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccagatgccg           >  2:1607841‑M1/1‑10 (MQ=255)
                                                               gggTGAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccagatgccg           >  2:763664‑M1/1‑10 (MQ=255)
                                                                   gAAGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccagatgccgtgtg       <  2:686506‑M1/76‑71 (MQ=255)
                                                                    aaGCAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccagatgccgtgtgt      >  2:1349999‑M1/1‑5 (MQ=255)
                                                                       cAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccagatgccgtgtgtgtc   >  1:1269408‑M1/1‑2 (MQ=255)
                                                                       cAatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccagatgccgtgtgtgtc   >  1:1233467‑M1/1‑2 (MQ=255)
                                                                        aatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccagatgccgtgtgtgtct  >  1:1621195‑M1/1‑1 (MQ=255)
                                                                        aatcgacatctacaacctgaccaaatacacccgttctaaccagaacacctgtatcaaccagatgccgtgtgtgtct  >  1:814791‑M1/1‑1 (MQ=255)
                                                                        |                                                                           
TGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTG  >  REL606/4163021‑4163168

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