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 41 [0] [0] 41 [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)41 (0.530) 37/494 0.7 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

GCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGCATCC  >  NC_000913/404571‑405064
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
tgacatgtattatacgtacgcttcaaaaacaaaattatcagcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttcaccataaattatatatgtccacattcggacttaggggaaataataattgcaccattcgtctgtcacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgatacttgatggatcaacttattacgtccctgaggagggatgacaaTGACGt                                                                                                                                                                                                                                                     >  2:68689‑M2/246‑251 (MQ=255)
                                aattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGtttt                                                                                                                                                                                                                     <  1:540018‑M2/38‑1 (MQ=255)
                                    atatgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaccttattacgtccctgaggagggatgacaaTGACGACCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTg                                                                                                                                                                                                                 >  2:612430‑M2/210‑251 (MQ=255)
                                      ctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGAt                                                                                                                                                                                                               >  2:208328‑M2/208‑251 (MQ=255)
                                         caacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGACCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGAtttt                                                                                                                                                                                                            >  1:11822‑M2/205‑251 (MQ=255)
                                                          gtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGATTTTTTCTGATTTTTCCATACACTTCGTc                                                                                                                                                                                             <  2:622265‑M2/62‑1 (MQ=255)
                                                            gcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATc                                                                                                                                                                                          <  2:109341‑M2/65‑1 (MQ=255)
                                                                           atgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTgc                                                                                                                                                                            >  1:462011‑M2/171‑249 (MQ=255)
                                                                                   ttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGa                                                                                                                                                                    >  1:105533‑M2/163‑249 (MQ=255)
                                                                                         ccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCttt                                                                                                                                                             <  1:162478‑M2/94‑1 (MQ=255)
                                                                                               attatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTcc                                                                                                                                                      >  2:36177‑M2/151‑251 (MQ=255)
                                                                                                 tatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTcc                                                                                                                                                      >  2:37977‑M2/149‑249 (MQ=255)
                                                                                                  atatatgtccacattcgggcttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTTAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGt                                                                                                                                                   <  1:249502‑M2/104‑1 (MQ=255)
                                                                                                      atgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTAcgc                                                                                                                                               <  1:596598‑M2/108‑1 (MQ=255)
                                                                                                              attcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaaaataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCTTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACCCGTACCGCt                                                                                                                                       >  2:339065‑M2/136‑251 (MQ=255)
                                                                                                                                                       aacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCa                                                                                               <  1:118112‑M2/156‑1 (MQ=255)
                                                                                                                                                       aacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCa                                                                                               <  1:208328‑M2/156‑1 (MQ=255)
                                                                                                                                                                 aatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCAGGTGAGCAGACCATATCTTTCAGTGCCCCCgg                                                                                     >  1:791‑M2/85‑250 (MQ=255)
                                                                                                                                                                    cgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGt                                                                                 >  2:598096‑M2/82‑251 (MQ=255)
                                                                                                                                                                         gcgcgaaacataatctgttttcaatgcattgataattgatggctcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCCGGCTTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGACTCCCGGTGAGCAGAACATATCTTTCAGTGCCCCCGGATgttctac                                                                            >  2:180483‑M2/77‑249 (MQ=255)
                                                                                                                                                                          cgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTccc                                                                           >  1:166770‑M2/76‑251 (MQ=255)
                                                                                                                                                                            cgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGt                                                                         >  1:503045‑M2/74‑251 (MQ=255)
                                                                                                                                                                               aacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTGTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCGCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTc                                                                      <  1:174265‑M2/181‑1 (MQ=255)
                                                                                                                                                                                         gttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGCGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCACTGGTGGTGCCTCCCGGTGAGCAGACCATCTCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCAtt                                                            >  2:614460‑M2/61‑251 (MQ=255)
                                                                                                                                                                                            ttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCAttttt                                                         >  1:177040‑M2/58‑251 (MQ=255)
                                                                                                                                                                                                   cattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGACGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAAc                                                  >  1:540995‑M2/51‑251 (MQ=255)
                                                                                                                                                                                                      tgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAAccc                                                <  1:37977‑M2/203‑1 (MQ=255)
                                                                                                                                                                                                       gataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCa                                               <  1:515497‑M2/204‑1 (MQ=255)
                                                                                                                                                                                                           attgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTAcc                                          <  1:295077‑M2/209‑1 (MQ=255)
                                                                                                                                                                                                           attgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGCTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGAACACTGCAGCACGGAAGCCTTTCTCTTCCAGTCCGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGCTTCCAGCACCATTTTTGCGGAAACCATTAcc                                          >  2:600399‑M2/43‑251 (MQ=255)
                                                                                                                                                                                                             tgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGCGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAATCCTTTCTCTTCCAGAACGCGTACCCCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATATTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGc                                        >  1:240663‑M2/41‑251 (MQ=255)
                                                                                                                                                                                                              gatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGcc                                       <  1:138565‑M2/212‑1 (MQ=255)
                                                                                                                                                                                                               atggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCt                                      >  1:379853‑M2/39‑251 (MQ=255)
                                                                                                                                                                                                                  gatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTg                                     >  1:554008‑M2/36‑249 (MQ=255)
                                                                                                                                                                                                                   atcaacttatttcgtacctgaggagggatgacaaTGACGTCCGGAGAAAGTCATCAGGGTTTGCTGAGTTTTTATGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGAACCGCTTCAATGGTGGGGTCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTTCCGCCTGAGc                                  >  2:544093‑M2/35‑251 (MQ=255)
                                                                                                                                                                                                                     caacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGCTTTGCTGCGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGAAAGCCTTTCTCTTCCAGTACGCATACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTGCAGTGCCACCGGAGGTTCTCCCGTTACCAGCACCATTTTTGCGGAACCccttaccgctggagagg                                >  2:292697‑M2/33‑243 (MQ=255)
                                                                                                                                                                                                                                gtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACCCTTCGTCATCGCTTCGATCACTGCATCACGGAAGCCGTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGACCCCATTACCGCCTGAGCGGCAAATTTATAc                     <  1:264911‑M2/230‑1 (MQ=255)
                                                                                                                                                                                                                                      gaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCAccatttttgcggaacccattaccgccgggccgccaactttatacccttgg               >  2:204918‑M2/16‑227 (MQ=255)
                                                                                                                                                                                                                                                gacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGc      >  2:248281‑M2/6‑250 (MQ=255)
                                                                                                                                                                                                                                                  caaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGCATc   >  2:412183‑M2/4‑251 (MQ=255)
                                                                                                                                                                                                                                                   aaTGACGTACGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTATTCCAGTACGCGTACCGCTTCAATGGCGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGCTGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGATCGGCAAATTTATACGCCTGGGCCCGGGGCATcc  >  2:194789‑M2/3‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
GCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGCATCC  >  NC_000913/404571‑405064

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