Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 2,016,437 Δ2 bp :: IS186 (–) +8 bp :: Δ1 bp coding (570‑577/687 nt) fliH → flagellar biosynthesis protein FliH

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2016437 =1 (0.020)75 (1.460) 69/482 0.0 98.6% coding (570/687 nt) fliH flagellar biosynthesis protein FliH
?NC_000913 2514274 = NA (NA)noncoding (2/1345 nt) IS186 repeat region
* ? NC_000913 = 20164442 (0.040)61 (1.190) 59/484 0.1 96.7% coding (577/687 nt) fliH flagellar biosynthesis protein FliH
?NC_000913 = 2515615 NA (NA)noncoding (1343/1345 nt) IS186 repeat region

GCAGAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGGC  >  NC_000913/2514025‑2514521
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
gCAGAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTc                                                                                                                                                                                                                                                          <  2:218961/250‑1 (MQ=255)
                 cGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGCAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTa                                                                                                                                                                                                                                         >  2:243084/1‑250 (MQ=255)
                  gCTCTCTCCGCGTGCTGCAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTa                                                                                                                                                                                                                                         >  1:172674/1‑249 (MQ=255)
                       ctcCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAggg                                                                                                                                                                                                                                     <  1:422401/248‑1 (MQ=255)
                                  gAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                                                                                                                                                                                                         <  1:274/249‑1 (MQ=255)
                                     ggCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTACTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                                                                                                                                    >  1:188779/1‑251 (MQ=255)
                                                     tCCTGGTTTAATTCGCTACCTTCTATTCAATCGGGATTATCGCAGGCGCAGTTAAAGGACTGAATGAAGAGCAAGGAACCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCCCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGGTGTGGTATTACGCCTGATATGATTTACCGTGc                                                                                                                                                                                                    >  2:502809/1‑251 (MQ=255)
                                                      ccTGGTTTCCGTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGAGATGATTTAACGTGcc                                                                                                                                                                                                   <  2:277086/251‑1 (MQ=255)
                                                              ccTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAACGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt                                                                                                                                                                                           >  2:100615/1‑251 (MQ=255)
                                                                  cGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGGATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                                                                                                                                                        <  2:442587/250‑1 (MQ=255)
                                                                         ttCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAt                                                                                                                                                                                >  2:9118/1‑251 (MQ=255)
                                                                                caaacGGGATTATCGCAGGGGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGTGGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCTCTAACTTAAGGGTTGTGGTTTTACGCCTTATATGATTTAACGTGCCTATGAGTTACTCTCACGATAACTGGt                                                                                                                                                                         <  2:535322/247‑1 (MQ=255)
                                                                                  aTCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa                                                                                                                                                                       <  1:9118/251‑1 (MQ=255)
                                                                                     gggATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGc                                                                                                                                                                     >  2:173370/1‑250 (MQ=255)
                                                                                          tATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTACCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAAc                                                                                                                                                                                                                                            >  1:514641/1‑174 (MQ=255)
                                                                                          tATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGTTCTGAAGCGGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTTCGTCAATCGCAGCACTTTTGCTGTAAGACCATACATAAAAAAGCCGGGTATAATTCCCATTAGCGCTAAc                                                                                                                                                                                                                                            <  2:514641/174‑1 (MQ=255)
                                                                                                      aGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTCAGGGTTGTGGTATTACGCCTGCTATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTgg                                                                                                                                                  >  2:357169/1‑251 (MQ=255)
                                                                                                            aGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGcc                                                                                                                                            >  1:196474/1‑251 (MQ=255)
                                                                                                                cTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCga                                                                                                                                         >  2:488711/1‑250 (MQ=255)
                                                                                                                          agCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGAt                                                                                                                               >  1:533350/1‑250 (MQ=255)
                                                                                                                           gCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATAc                                                                                                                             <  1:357169/251‑1 (MQ=255)
                                                                                                                                gTAACGTGGTTTCTCGCGTCGGTCTGAAGCTGGTTTACTGCTCTACCCTGGTGATGGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGGATTACTCCTGATATGATTTAACGTGCCTATGAATTACTCTCACGATAACTGTGCAGGAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCgg                                                                                                                        <  1:51025/251‑1 (MQ=25)
                                                                                                                                       ggTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa                                                                                                                  <  1:488711/250‑1 (MQ=255)
                                                                                                                                        gTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt                                                                                                                <  1:167789/251‑1 (MQ=255)
                                                                                                                                                  tCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGGTTTGGCTGTCTGCGTCAATCGCAGCACTGGTTCTGGAAGTCCATACATAAAAAAGCCGGGTATAATTCCCATAAGCGCTACCTTAAGTGTTTTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCTGCAATTCTGGCGCATATCGGTAAGCCCGAAGACCTGGATACTTCGGCACGTAGTGCCGGGGct                                                                                                       <  2:488722/250‑1 (MQ=18)
                                                                                                                                                         gAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATAGGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCTCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAAcc                                                                                                 <  1:373088/249‑1 (MQ=255)
                                                                                                                                                           aGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACccgcc                                                                                             <  2:172674/251‑1 (MQ=255)
                                                                                                                                                                    aCGGCTCTACCCTGGTGAGTGTGCTGGCTGCGTCAATCGCATCACTGGTGCTGTAAGACCATACATAAAAAAGGCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTTTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGTTCAGCAATCCGTGGCCATATTGGTAAGCCCGCAGAAGTGGATACTTCGGCACGTAATGCCGGGATTCTGACCCGCCGCCGCGaa                                                                                     <  2:65588/250‑1 (MQ=18)
                                                                                                                                                                       gCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTc                                                                                 >  1:22918/1‑251 (MQ=255)
                                                                                                                                                                                 ggTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGTTACTTCGGCACGTAAGGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGAtgct                                                                         <  2:42183/249‑1 (MQ=255)
                                                                                                                                                                                      gtgtGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAActc                                                                  >  2:341792/1‑251 (MQ=255)
                                                                                                                                                                                         tgCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATCCGTGCTGCTGCACCTCTGc                                                               >  2:101116/1‑251 (MQ=255)
                                                                                                                                                                                                  gCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTc                                                         <  1:178442/249‑1 (MQ=255)
                                                                                                                                                                                                                      gCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCgg                                   >  1:99059/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTg                  >  1:116188/1‑250 (MQ=37)
                                                                                                                                                                                                                                          aaaGCCGGGGATAATGCCCATAAGCGCTCTCTTAAGGGTTGTGGTATTACGCCTGTTATGTTTTAACGTGCTGTGGAATTGCTCGCACGATAACTGGTAAGCAATGCTGTCTCAGATTGGTAAGCCCGAAGACCTGGATTCTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGGCGCGAAATTCGTGAGGCTGCAACTCTGCTACGTTGGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAg               <  2:394035/251‑1 (MQ=11)
                                                                                                                                                                                                                                                       aTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGGc  >  1:4157/1‑251 (MQ=11)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
GCAGAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGGC  >  NC_000913/2514025‑2514521

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