Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 390,251 Δ14,565 bp [yaiT][dgcC] 16 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 390251 404815 14565 44 [0] [0] 44 [yaiT]–[dgcC] [yaiT],ytiB,insF2,insE2,yaiT,iprA,ampH,sbmA,yaiW,yaiY,yaiZ,ddlA,iraP,phoA,psiF,[dgcC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3902500 (0.000)44 (0.560) 41/494 0.5 100% intergenic (‑523/‑1) hemB/yaiT porphobilinogen synthase/putative autotransporter YaiT
?NC_000913 404816 = 0 (0.000)coding (1114/1116 nt) dgcC diguanylate cyclase DgcC

AAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCG  >  NC_000913/390032‑390497
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       
aaaCAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctg                                                                                                                                                                                                                         >  1:293649‑M1/1‑219 (MQ=255)
  aCAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctg                                                                                                                                                                                                                         >  1:198857‑M1/1‑217 (MQ=255)
       aatatCTGCAACGTGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACAGTCGGACTTAGGGGAAAGAATAATTTATCCATTCGTCTGTAACGCAGAATAATCGTTAGCGCGAAACATAATATGGTTTCTATGCATTGATAATTGATGGGTCAACTTATTACGACAGTGAGGAGGTAGAACAAtgacgaccggcgaaagtcatctgggtttgctgagttttg                                                                                                                                                                                                                  <  2:209078‑M1/249‑40 (MQ=255)
                      ttGAAATCCTGAGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATCTGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatac                                                                                                                                                                                                   >  2:445911‑M1/1‑197 (MQ=255)
                                  gCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatc                                                                                                                                                                                        <  2:254869‑M1/250‑66 (MQ=255)
                                   cTTTTATGACTTTGATGCGTGTTTTTCTCCATAAATTATATATGTTTACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAttacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgc                                                                                                                                                                                      >  2:58810‑M1/1‑184 (MQ=255)
                                                 aTGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgc                                                                                                                                                                          >  1:77141‑M1/1‑170 (MQ=255)
                                                   gCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcg                                                                                                                                                                                                                                             >  1:280831‑M1/1‑168 (MQ=255)
                                                   gCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcg                                                                                                                                                                                                                                             <  2:280831‑M1/180‑13 (MQ=255)
                                                             ctcCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaa                                                                                                                                                                                                                                           <  1:26543‑M1/172‑15 (MQ=255)
                                                             ctcCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaa                                                                                                                                                                                                                                           >  2:26543‑M1/1‑158 (MQ=255)
                                                                 aTAAATTATATATGTCCACATTCGGACTTAGGGGAAAGCATAATTGAACAATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctc                                                                                                                                                        >  2:479747‑M1/1‑154 (MQ=255)
                                                                        ctatatgTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggcagcctttctcttccagc                                                                                                                                                 >  2:633424‑M1/2‑147 (MQ=255)
                                                                               tCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgcaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtc                                                                                                                                          >  2:136826‑M1/1‑140 (MQ=255)
                                                                                    aTTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcg                                                                                                                                                                                  >  1:23184‑M1/1‑135 (MQ=255)
                                                                                    aTTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcg                                                                                                                                                                                  <  2:23184‑M1/206‑72 (MQ=255)
                                                                                              aGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtgg                                                                                                                           >  1:366831‑M1/1‑125 (MQ=255)
                                                                                                        ataatTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccgg                                                                                                                 <  2:301640‑M1/251‑137 (MQ=255)
                                                                                                        ataatTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccgg                                                                                                                 >  2:12169‑M1/1‑115 (MQ=255)
                                                                                                              gAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagc                                                                                                            <  2:156534‑M1/250‑142 (MQ=255)
                                                                                                               aaCCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagc                                                                                                            >  2:25998‑M1/1‑108 (MQ=255)
                                                                                                                      cGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgctgcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccata                                                                                                    <  2:381406‑M1/250‑150 (MQ=255)
                                                                                                                       gTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcacttccgcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatc                                                                                                  >  2:619857‑M1/1‑100 (MQ=255)
                                                                                                                       gTCTGTAACGCAGCATAAGCGTTAGCGCGAATCATAATAAGGTTTCTATGCAGTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGTGATGACAAtgacgtccggcgaaagtcagcaggatttgctgagtttttctgatttgtccatacacttcgtcatcgcttttatcacggcagcacggaggcctttctctttcagtccgcgtaccgcgacaatggtggtgcctcccggttaccagaccatatc                                                                                                  <  2:515926‑M1/251‑152 (MQ=255)
                                                                                                                         cTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctt                                                                                                >  1:321554‑M1/1‑98 (MQ=255)
                                                                                                                           gTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatcttt                                                                                               <  2:445691‑M1/250‑155 (MQ=255)
                                                                                                                            tAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttc                                                                                              >  1:129554‑M1/1‑95 (MQ=255)
                                                                                                                                               gcgcGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcc                                                                          >  1:514809‑M1/1‑76 (MQ=255)
                                                                                                                                               gcgcGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCACTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggcggtgcctcccggtgagcagacaatatctttcagtgcccccggatgttctac                                                                          >  2:433519‑M1/1‑76 (MQ=255)
                                                                                                                                                    aaaCATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgttt                                                                     >  1:244339‑M1/1‑71 (MQ=255)
                                                                                                                                                       cATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttcc                                                                   >  1:325519‑M1/1‑68 (MQ=255)
                                                                                                                                                          aataTGTTTTCTCTGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagc                                                                >  2:75267‑M1/1‑65 (MQ=255)
                                                                                                                                                             atGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGTAGGGATGACAAtgacgtccggcgaacgtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccc                                                            >  2:276727‑M1/1‑62 (MQ=255)
                                                                                                                                                             atGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcacca                                                            <  1:108647‑M1/251‑190 (MQ=255)
                                                                                                                                                                      aTGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgc                                                     <  1:447386‑M1/249‑197 (MQ=255)
                                                                                                                                                                              aTAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccatt                                           >  2:50915‑M1/1‑45 (MQ=255)
                                                                                                                                                                                     aTGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcct                                    <  1:17643‑M1/251‑214 (MQ=255)
                                                                                                                                                                                       ggATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctg                                   <  1:161501‑M1/250‑215 (MQ=255)
                                                                                                                                                                                           cAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgagcgg                              >  2:365873‑M1/1‑32 (MQ=255)
                                                                                                                                                                                             aCTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgttccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgagcggca                            <  1:611134‑M1/251‑222 (MQ=255)
                                                                                                                                                                                             aCTTATTACGTCCATGAGGAGGGATGACAAtgacgtccggcgacagtcatcaggaattgatgagtttttctgatttttccatacacttcgtcatcgcttcgataactgcagcacggaagcctttctcttccagtacgcgtaccgcttcactggtgctgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgagcggc                             >  1:489705‑M1/1‑30 (MQ=255)
                                                                                                                                                                                                   taCGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgagcggcaaattt                       >  1:18753‑M1/1‑24 (MQ=255)
                                                                                                                                                                                                                    aTGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggaggttctcccgtttccagcaccatttttgcgtaacccattacagcctgcgcggctaatttatacgactaggcgcgtgg     >  2:566415‑M1/1‑7 (MQ=255)
                                                                                                                                                                                                                       aCAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgagcggcaaatttatacgcctgggcgcgtggcat  <  2:198857‑M1/251‑248 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       
AAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATACCG  >  NC_000913/390032‑390497

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