Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,543,608 Δ1 bp :: IS186 (–) +6 bp :: Δ1 bp coding (494‑499/549 nt) fimA → type 1 fimbriae major subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2514274 =NA (NA)28 (0.660) 27/480 0.6 87.2% noncoding (2/1345 nt) IS186 repeat region
?NC_000913 4543608 = 4 (0.100)coding (494/549 nt) fimA type 1 fimbriae major subunit
* ? NC_000913 = 2515616NA (NA)41 (0.960) 37/480 0.2 90.9% noncoding (1344/1345 nt) IS186 repeat region
?NC_000913 = 4543613 4 (0.100)coding (499/549 nt) fimA type 1 fimbriae major subunit

TGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGCCCGGCGGAAATATGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATT  >  NC_000913/4543364‑4543862
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
tGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGGACGGCGATTGATGCGGGTCATACCAACGTTCAGGCTCTGCAGAGTTCCGCTGCGGGTAGCGCAACAACCGTTGGGGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACACCCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                                                                                                                                                                                                          >  1:43684/1‑251 (MQ=255)
ggCGATACCAATGGTGCCTATAAAGCCGCGGTGTTCGTTTTAGGGACGGCGATTGATGCGGTGCATACCATAGTTCGTGCTCTGGAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTCGGCAGAGCCTGGACAGAACGGGTGTTGTGCTGACGCTGGATGGTGCGGTATTGAGTTCAGAAGCAACCCTGAATAGCGGACCTAATACCATGCCTTTCCAGGCGCGTTATTGTGCACCCGGGGCCGAAAc                                                                                                                                                                                                                                                          <  2:207082/250‑1 (MQ=255)
  gcagtgacttcacgtaatgacatccccccggggccgtaagccagccccagacgtagcagagttgcagcatcacgaatttcgcggcggcgggttagagccccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAAcc                                                                                                                                                                                                                                                         <  1:195295‑M2/8‑1 (MQ=255)
   gATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccc                                                                                                                                                                                                                                                       >  2:164470/1‑251 (MQ=255)
               gCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                                                                                                                                                                                                            >  1:256585‑M2/1‑235 (MQ=255)
                       aGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATGCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                                                                                                                                                                                                   >  1:214851‑M2/1‑227 (MQ=255)
                         ccgcttgtgccttttTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATTCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttga                                                                                                                                                                                                                                 >  1:142238‑M2/8‑225 (MQ=255)
                            cccggggccgtaagccagccccagacgtagcagagttgcagcatcacgaatttcgcggcggcgggttagagccccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                                                                                                                                                                                                              >  2:378895‑M2/217‑251 (MQ=255)
                            cccggggccgtaagccagccccagacgtagcagagttgcagcatcacgaatttcgcggcggcggggtagagccccggcattacgtgccgaagtatccagatcttcgggcttaccaatatgggccagaattgctgacctgttatcgtgagagtaattcataggcacgttaaatcatatcaggcgtaataacacaccccttaagttagcgcttatgggACGCAACCCCGGGTGCTGCTAAAGCGGAAGCGAcc                                                                                                                                                                                                                              >  2:95988‑M2/217‑251 (MQ=255)
                               ttGCCTTTTTAGGGACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAACCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTTCAACCGGGGCCGCAAcccataagcgcacacttaagggttgaaccatc                                                                                                                                                                                                                           <  2:460866‑M2/251‑33 (MQ=255)
                                   cTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGATCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctga                                                                                                                                                                                                                        <  1:427269‑M2/250‑36 (MQ=255)
                                     gtaagccagccccagacgtagcagagttgcagcatcacgaatttcgcggcggcgggttagagccccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                                                                                                                                                                                                     >  1:127005‑M2/208‑251 (MQ=255)
                                         gccagccccagacgtagcagagttgcagcatcacgaatttcgcggcggcgggttagagccccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                                                                                                                                                                                                 <  1:45310‑M2/48‑1 (MQ=255)
                                                    acggagcagagtggcagcatcacgaattgcgcggcggcgggttagcgccccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGTGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                                                                                                                                                                                                      <  2:245888‑M2/59‑1 (MQ=255)
                                                     ttaggcgggTCATACCACCGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcttacttaagggttgaaccagctgaagaattcgacgcctcggt                                                                                                                                                                                                      >  2:451387‑M2/5‑197 (MQ=255)
                                                         gCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctc                                                                                                                                                                                                 >  1:460866‑M2/1‑193 (MQ=255)
                                                               cATACCAACGTTCTGTCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGCTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccagaagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaag                                                                                                                                                                                           >  1:288487‑M2/1‑187 (MQ=255)
                                                                agaccAACGTTCTGGCTCGGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgatgcctcggggcctcgttaag                                                                                                                                                                                           >  2:50344‑M2/3‑186 (MQ=255)
                                                                    cAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgat                                                                                                                                                                                      >  1:428306‑M2/1‑182 (MQ=255)
                                                                          tCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGAGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcg                                                                                                                                                                                <  1:445863‑M2/251‑76 (MQ=255)
                                                                                    cAGAGGTCATCTTCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGCGTGCTGCGCTGACGCGGGATGGTGCGGCAGTTAGTTCAAAAACACCCCTGAATAACGGATCCAATTCCATTCCGTTCCAGGCGCGTTATTTTGAAACCGGGGCCGCAAcccataagcgctaacttaagggttgcaccagctgaagaatgcgaagcctcggtgccgcgttaagacgatgccgcgcgttcttcta                                                                                                                                                                      >  2:19747‑M2/1‑166 (MQ=255)
                                                                                     agagTTCAGCTGCGGGAAGCGCAACAAACGGTGGTGTGCAGATCCTGGAAAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGCAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcatcgcgttcttcaat                                                                                                                                                                     <  1:37398‑M2/251‑87 (MQ=255)
                                                                                            aGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttt                                                                                                                                                              >  1:433659‑M2/1‑158 (MQ=255)
                                                                                             gcagcggGTAGCGCAACATACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgtttt                                                                                                                                                             >  2:100622‑M2/4‑157 (MQ=255)
                                                                                                gCGGGTATCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtggctcgttaagacgatggctcgcgttcttcaattgcgttttgta                                                                                                                                                          >  1:175188‑M2/1‑154 (MQ=255)
                                                                                                 cGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtag                                                                                                                                                         >  1:397528‑M2/1‑153 (MQ=255)
                                                                                                   ggTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggc                                                                                                                                                       <  1:62851‑M2/251‑101 (MQ=255)
                                                                                                    gTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggct                                                                                                                                                      >  1:267522‑M2/1‑150 (MQ=255)
                                                                                                     gagCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAAGAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctg                                                                                                                                                     >  2:384035‑M2/2‑149 (MQ=255)
                                                                                                             caaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggat                                                                                                                                             >  1:218929‑M2/1‑141 (MQ=255)
                                                                                                                         atccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGaaaa                                                                                                                                  <  2:127005‑M2/127‑1 (MQ=255)
                                                                                                                              gttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaactcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCTCCCTTGTGCGATAAAAATAACGATGAAAAGGAaga                                                                                                                            <  1:142004‑M2/133‑1 (MQ=255)
                                                                                                                                   tcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGat                                                                                                                         <  1:13009‑M2/136‑1 (MQ=255)
                                                                                                                                   ggACATAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctctcattcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacc                                                                                                                        >  1:184674‑M2/1‑119 (MQ=255)
                                                                                                                                       gcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTc                                                                                                                   <  1:224139‑M2/142‑1 (MQ=255)
                                                                                                                                        aaacGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTATTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaag                                                                                                                  >  1:64299‑M2/2‑114 (MQ=255)
                                                                                                                                          taccaatatgggccagaattgctgaccatttatcgggagagtaattcatcgtcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttaggggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCt                                                                                                                  <  2:249658‑M2/143‑1 (MQ=255)
                                                                                                                                           cGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACACCCCTGAATAACGGAACCAATACCATTCCGTGCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttaaaccatctgaagaatgcgacgcctcggtgcctcgtgaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgttccacgaatggccacttgtaagctc                                                                                                               >  2:115912‑M2/1‑111 (MQ=255)
                                                                                                                                            ccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCCCCCTTGTGCGAt                                                                                                                                                   >  2:128888‑M2/105‑214 (MQ=255)
                                                                                                                                            ccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGAt                                                                                                                                                   <  1:128888‑M2/110‑1 (MQ=255)
                                                                                                                                            ccaatatgggccagaattgctgaccagttatcgtgagagtaactcatcggcactttaaatcatatcaggcgtaataccacatcccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTa                                                                                                              <  1:86367‑M2/147‑1 (MQ=255)
                                                                                                                                              gTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAACCAACCCTGAATAACGCGACCAATCCCATTCCCGTCCCGTCGCGTTACTTTGCAACCGGGGCCGCAAcccataagcgcgaacttaagggttgaaccacctgaagaacgcgacgcctcggtgcctcgttaagacgatgcctcgcgtgcttcaattgcgtgttgtaggctgtccgggatacggtcccacgaatggccccctgtaagctctag                                                                                                            >  2:37398‑M2/1‑108 (MQ=255)
                                                                                                                                                   gcgcTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatga                                                                                                        >  2:62851‑M2/1‑103 (MQ=255)
                                                                                                                                                           attgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggaCCGCAACCCCGGGTGCTGCTAATGCGGATGCGCACTTCAAGGTTCCGTATCAATAACCTCCCCAGGTTCAGGGACGTCATTACTGGCAGGGATGCCCACCCTTGTGCGATAAAAATATCGATGCAACGGAGGAGTTTATTTCTATTAGCTTCGTTGCTgccc                                                                                               >  2:336819‑M2/90‑250 (MQ=255)
                                                                                                                                                             tGGCTGGAGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccctctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagctgcccatttttgtt                                                                                             <  2:374068‑M2/251‑159 (MQ=255)
                                                                                                                                                              ggATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttt                                                                                            <  1:178067‑M2/251‑160 (MQ=255)
                                                                                                                                                                 gaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGttt                                                                                         >  1:419019‑M2/84‑251 (MQ=255)
                                                                                                                                                                  accagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTg                                                                                        <  1:43569‑M2/169‑1 (MQ=255)
                                                                                                                                                                     gCGACATTTAGTTCAGAAACAACCCGGAATAACGGAACCAATACCAGTCCGTTCCAGGCGCGTTATTTTGCAATCGGGTCCGCAActcataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcgtcgttaagacgatgcctcgcgttcttcaattgcgctttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattcgcca                                                                                     >  2:177496‑M2/1‑85 (MQ=255)
                                                                                                                                                                            gtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGAGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCgg                                                                              >  1:336264‑M2/73‑251 (MQ=255)
                                                                                                                                                                               gTTCAGAAACATCCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATGTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagtta                                                                           >  1:133537‑M2/1‑75 (MQ=255)
                                                                                                                                                                               agagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGACACGGACGAAATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGaa                                                                            >  2:420173‑M2/70‑250 (MQ=255)
                                                                                                                                                                                    gAAACTACCCTGATTAACGGAACCAATGACATTTCGTTCCAGGCGCTCTATTTTGCAACCGGGGCCGCAAcccataagtgctaacttaagggttgaactagctgaagaatgcgacgcctcggtgccttgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagtgatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccaccacgagttagttct                                                                      >  2:317340‑M2/1‑70 (MQ=255)
                                                                                                                                                                                       tcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCCTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATgg                                                                    >  2:344198‑M2/62‑250 (MQ=255)
                                                                                                                                                                                                 aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaccttaagggttgcaccatctgcagaatgcgacgcctcggtgcatcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggttactgtcccacgaatggccacctgtaagctccagatgaccat                                                                                                    <  2:284043‑M2/208‑152 (MQ=255)
                                                                                                                                                                                                 aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccat                                                                                                    >  1:284043‑M2/1‑57 (MQ=255)
                                                                                                                                                                                                      atcatatcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCCCGTTGCCcg                                                    >  2:20604‑M2/47‑251 (MQ=255)
                                                                                                                                                                                                         aCCAATACCATTACGCTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagagttgaaccatctgcagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttggtattctccacaacgagttagttcttcttttcggctccggcacttc                                                 <  2:426930‑M2/251‑203 (MQ=255)
                                                                                                                                                                                                           cAATACCATTCCGTTCCAGGCGCGTTATGTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttctacaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggagccggcacttctg                                               >  2:43684‑M2/1‑47 (MQ=255)
                                                                                                                                                                                                           atcaggcgtaataccacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGCCcggcgg                                                >  1:177635‑M2/42‑250 (MQ=255)
                                                                                                                                                                                                              aggcgtaataccacaacccttacgttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGACCGGCGGAAat                                            >  2:99325‑M2/39‑251 (MQ=255)
                                                                                                                                                                                                                        ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagcatgagacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatccgtcacttctgggggggaaatcc                                   <  2:272338‑M2/250‑217 (MQ=255)
                                                                                                                                                                                                                         cacaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTCATGCGGATGCGCCCTTCAAGGTTCAGTATCACTAACCTACCCAGGTTCAGGGACGTCATTACGGGCCGGGAAGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGCCCGGCGGAAATATGCAATTTc                                  >  2:177342‑M2/28‑250 (MQ=255)
                                                                                                                                                                                                                          ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccagaagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatccggcacttctgggggggaaatccagc                                >  2:93672‑M2/1‑32 (MQ=255)
                                                                                                                                                                                                                               ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGCCCGGCGGAAATATGCAATTTCAGGGCGt                           <  1:20604‑M2/230‑1 (MQ=255)
                                                                                                                                                                                                                                gcgatATTTTGCAACCGGGGCCGCAACCCCGGTTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGGATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGACCGGAAATAAATGGAATACCACGTTGCCCGGCGGAAATATGCAATTTAAGGGCGTc                          >  2:285125/5‑251 (MQ=255)
                                                                                                                                                                                                                                   aagtgtgcgagtaatggCCGCCACCCCGGGGGCTGAGAATGCGGAGGCGACATTCAGGGTTACGTGTCAACACCCTACCCAGGTTCAGGGACGTCATGACGGGCGGGGATGCCCACTCTCGGGCGATAAAAATAACGAAGAAAAGGAAGAGATGATTTCTATTAAGGTCGTTGAGGCCAATGTTTGCTCTGGGCGGAAATAAAAGTAATACCACGTTGCCCGGCGGAAATATGCAAGTTCAGGGCGTCatt                       <  2:445318/236‑1 (MQ=255)
                                                                                                                                                                                                                                     aTTTTGCAACCGTGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatccggcacttctgggggggaaatccagcgatggctgga                      >  1:453768‑M2/1‑21 (MQ=255)
                                                                                                                                                                                                                                        agcgcttatggtCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGCCCGGCGGAAATATGCAATTTCAGGGCGTCATTATTg                   <  2:257629‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                                         gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGCCCGGCGGAAATATGCAATTTCAGGGCGTCATTATTGCg                 >  1:429502‑M2/12‑251 (MQ=255)
                                                                                                                                                                                                                                           gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGCCCGGCGGAAATATGCAATTTCAGGGCGTCATTATTGCGGa               <  1:91662‑M2/242‑1 (MQ=255)
                                                                                                                                                                                                                                               cGGGGCCGCAAcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatccggcacttctgggggggaaatccagcgatggctggattatgtcgtc            <  1:389858‑M2/250‑240 (MQ=255)
                                                                                                                                                                                                                                                 gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGCCCGGCGGAAATATGCAATTTCAGGGCGTCATTATTGCGGAAACtt          <  1:398606/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        aacccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatccggcacttctgggggggaaatccagcgttggctggattatgtcgtcaagctagga   <  2:58277‑M2/250‑249 (MQ=255)
                                                                                                                                                                                                                                                        aacccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatccggcacttctgggggggaaatccagcgatggctggattatgtcgtca           >  1:58277‑M2/1‑2 (MQ=255)
                                                                                                                                                                                                                                                        aaCCCCGGGTGCTGCTAATGAGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAACGGAGGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAACTGGAATACCACGTTGCCCGGCGGCAATATGCCATTTCAGGGCGTCATTCTTGCGGAAACTTGCCGGAtt  >  2:462673/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
TGCGATACCAATGTTGCATCTAAAGCCGCTGTTGCCTTTTTAGGTACGGCGATTGATGCGGGTCATACCAACGTTCTGGCTCTGCAGAGTTCAGCTGCGGGTAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCGTCGTTGCTGCCAATGTTTGCTCTGGCCGGAAATAAATGGAATACCACGTTGCCCGGCGGAAATATGCAATTTCAGGGCGTCATTATTGCGGAAACTTGCCGGATT  >  NC_000913/4543364‑4543862

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