Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 2,176,844 IS4 (–) +13 bp coding (349‑361/855 nt) gatY ← tagatose‑1,6‑bisphosphate aldolase 2

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2176844 =3 (0.120)23 (0.900) 20/484 0.3 87.9% coding (361/855 nt) gatY tagatose‑1,6‑bisphosphate aldolase 2
?NC_000913 = 4503515 NA (NA)noncoding (1/1426 nt) IS4 repeat region
* ? NC_000913 = 21768563 (0.120)18 (0.700) 18/486 0.4 85.0% coding (349/855 nt) gatY tagatose‑1,6‑bisphosphate aldolase 2
?NC_000913 4502090 = NA (NA)noncoding (1426/1426 nt) IS4 repeat region

GCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTGAGAAACTCAGTTGATTCGGCC  >  NC_000913/4501841‑4502329
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
tcgcTAAAAAATGCTTTATTTTTGAGTACCTCGAACCCTTGCCAATTAAGGGACATAGCCAGAGACGCAACTGGCTGATGCAAGCAGGCTAATGTCTTTAGAGCCGCTGGCCGGCGGGCATGACAATGTCAAACAGTCGAAGATTGTGGCGTTGACGGTCATTCTTGGGGACAAAAAATCTCAGAATTCTCGTTGGCAACCAGATAATTTAACAGTTTATGCGGAAACTAAATAGGCATCTGGGTTCTGt                                                                                                                                                                                                                                                 <  2:250920/249‑1 (MQ=255)
   cTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGAAACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATAGGCAGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGAGGGGAAGTCTTGGGGACAAAAAAGCTCAGTAATATCGTTGACAACCAGGTAATTAAAGAGTTGATGCTGAAAAGAGATAGGCATCTTGGTTATGGAATg                                                                                                                                                                                                                                             >  1:14127/1‑251 (MQ=255)
     aaaaaaTGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGATGGTCAGTCTTGGGGACAAAAAAGTTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGcc                                                                                                                                                                                                                                           <  1:100193/251‑1 (MQ=255)
                agtctagaaaaatcaagcgccggtgcgctggcatacatcccatgagccgtgccgatggcgaccgccagggaatcaattccggttgcctcggcaaattcacgcgcctgagcggggttggtgtacaacgcatcggcttcattgacttgcacatcatcttcctggccgccaagttgccccagctccgcctcgacgctgacatcaaagcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTa                                                                                                                                                                                                                                  >  1:54568‑M2/234‑249 (MQ=255)
                                     ggtgcgctgtcatacatcccatgagccgtgccgatggcgaccgccagggaatcaattccggttgcctcggcaaattcacgcgcctgagcggggttggtgtacaacgcatcggcttccttgacttgcacatcatcttcctggccgccaagttgccccagctccgcctcgacgctgacatcaaagcgctggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCg                                                                                                                                                                                                           >  2:251477‑M2/213‑251 (MQ=255)
                                                        ttatttaccggagcaatacccaTGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTc                                                                                                                                                                                          <  1:18122/228‑1 (MQ=255)
                                                        catgagccgtgccgatggcgaccgccagggaatcaattccggttgcctcggcaaattcacgcgcctgagcggggttggtgtacaacgcatcggcttcattgacttgcacatcatcttcctggccgccaagttgccccagctccgcctcgacgctgacatcaaagcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTcc                                                                                                                                                                                         >  1:115385‑M2/194‑250 (MQ=255)
                                                         atgagccgtgccgatggcgaccgccagggaatcaattccggttgactcggcaaattcaagcgcctgagcggggttggtgtacaacgcatcggcttcattgacttgcacatcatcttcctggacgccaagttgccccagctccgcctcgacgctgacatcaaagcgatggaaaaaatccaccccctctttgacTAATGACAGTCAGTTAAGCAactgaatggctctttttcgtggatttggggagtttcc                                                                                                                                                                                         >  2:285080‑M2/193‑231 (MQ=255)
                                                           cAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGTGCTGTGGGGTATTTCCAggg                                                                                                                                                                                     <  2:173963/251‑1 (MQ=255)
                                                                aCGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCGGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCctc                                                                                                                                                                                 <  1:58380/250‑1 (MQ=255)
                                                                  gccgatggcgaccgccagggaatcaattccggttgcctcggcaaattcacgcgcctgagcggggttggtgtacaacgcatcggcttcattgacttgcacatcatcttcctggccgccaagttgccccagctccgcctcgacgctgacatcaaagcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTcc                                                                                                                                                                              <  1:251477‑M2/68‑1 (MQ=255)
                                                                  gCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTggg                                                                                                                                                                                                 <  1:140025/232‑1 (MQ=255)
                                                                  gCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTggg                                                                                                                                                                                                 >  2:140025/1‑232 (MQ=255)
                                                                           tGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGACAAAGATTGTGGCGTTGACGGTCAGGCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAAAAGTTGATGCTGAAACGAAATAGGAATCTTGGTTCTGTAATGCCAGTCAGATAAGCAACTGACAGGATCTTTTTCGGGGCTGTGGGGAATAGCCAGGGCCACTCCTTTACAAct                                                                                                                                                                     >  2:252679/1‑251 (MQ=255)
                                                                           gaccgccagggaatcaattccggttgcctcggcaaattcacgcgcctgagcggggttggtgtacaacgcatcggcttcattgacttgcacatcatcttcctggccgccaagttgccccagctccgcctcgacgctgacatcaaagcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCACGGGCCTCTCCTTTACCAct                                                                                                                                                                     >  2:142896‑M2/175‑251 (MQ=255)
                                                                               gCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCATCTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACtctc                                                                                                                                                                  <  1:190519/250‑1 (MQ=255)
                                                                                         tAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGACTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGAAACTGACTGGCTCTTTTTCGTGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGCAGGCCCt                                                                                                                                                       >  2:46813/1‑251 (MQ=255)
                                                                                             cccggtttccgcggcaaattcacgcgcctgagcggggttggtgtacaactcatcggcttcattgacttgcacatcatctgcctggccgccaagttgccccagctccgcctcgacgctgacatcaaagcgatggcaaaaagccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTTTGGGGTATTTCCAGGGCCGCTCCTTTACCACTCTCGGGTTGGCCCTTTcc                                                                                                                                                   <  2:182170‑M2/95‑1 (MQ=255)
                                                                                                    gATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTg                                                                                                                                            <  1:183937/251‑1 (MQ=255)
                                                                                                                     gCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGt                                                                                                                           >  1:208022/1‑251 (MQ=255)
                                                                                                                                      aGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGTAGGCCCTTTCCCTTCTTGTCGGTCATTTCTCAAGTTGTCCCATACttaca                                                                                                            >  2:139651/1‑246 (MQ=255)
                                                                                                                                                            ggTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGtat                                                                                     >  1:142303/1‑250 (MQ=255)
                                                                                                                                                                  tCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTccc                                                                              >  1:213967/1‑251 (MQ=255)
                                                                                                                                                                    catcatcttcctggccgccaagttgccccagctccgcctcgacgctgacatcaaagcgatggcaaaaatccaccaccgctttgacTAATGCCATTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCgg                                                                            <  1:218545‑M2/166‑1 (MQ=255)
                                                                                                                                                                                    gccaagttgccccagctccgcctcgacgctgacatcaaagcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAAt                                                            >  1:232193‑M2/70‑251 (MQ=255)
                                                                                                                                                                                       gTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCACTGtc                                                         >  2:270604/1‑251 (MQ=255)
                                                                                                                                                                                       aagttgccccagctccgcctcgacgctgacatcaaagcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGCATACGTCCCGGTGATGCGCCCTGCAATGtc                                                         >  1:252388‑M2/67‑251 (MQ=255)
                                                                                                                                                                                         aTTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTCAACGGCCCGGTGACGCGCCCTACAATGtcat                                                       >  2:270119/1‑251 (MQ=255)
                                                                                                                                                                                           tgccccagctccgcctcgacgctgacatcaaagcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGtcatc                                                      <  1:142896‑M2/188‑1 (MQ=255)
                                                                                                                                                                                            ctcGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGc                                                                                                                                                          >  1:185933/1‑149 (MQ=255)
                                                                                                                                                                                            ctcGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATGTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGc                                                                                                                                                          <  2:185933/149‑1 (MQ=255)
                                                                                                                                                                                                   ctccgcctcgacgctgacatcaaagcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCCCAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATtctc                                             >  1:270774‑M2/55‑251 (MQ=255)
                                                                                                                                                                                                     ccgcctcgacgctgacatcaaagcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATtctc                                             <  1:72550‑M2/197‑1 (MQ=255)
                                                                                                                                                                                                     ccgcctcgacgcggacatcaaagcgatggcataaatccaccacctctttagcTAATGACAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGAGCCTCTCCTTTACCACTCTCGGGCAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGAATTctcac                                           >  2:251623‑M2/53‑250 (MQ=255)
                                                                                                                                                                                                         aGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCCCAAGTTGTCCCATACTTGCAAGATCGCGCATCAACTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCAcc                                       >  1:290767/1‑251 (MQ=255)
                                                                                                                                                                                                         aGATAATTTAACAGTAGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCAAGTAAGTTAAGCAACTGACTGGCTCTTTTTCGAGGCTGTGGGGTATTTACAGGGACTCTCCTTTACCACTCTCGAGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCAGCAGCTACGGTATACGTCACGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCAcc                                       >  1:184545/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaagcgatggcaaaaatccaccacctccttgacGAATTCCAGTCAGTTAAGCAACTGACTGTCTCTTTTTCGGGGCTGTGGGGTACTACCATGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATTCTTGCAAGATCGCGCATCAGCTCCTGTATACGTCCTGGTGAAGCGCCCTGCAATGTCATCAGCATACTCAGCACCATTCCGCATGATTCt                        <  2:206870‑M2/218‑1 (MQ=255)
                                                                                                                                                                                                                         gATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGACCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGAGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTg                       <  1:46813/251‑1 (MQ=255)
                                                                                                                                                                                                                          agcgatggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTga                      <  1:15078‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                                          aagCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCGCTTTTTCGGGGCTGTGGGGTATTTCCAGTGCCGCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGGCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTACACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCTCCATCCCGCATGATTCTga                      <  2:142303/249‑1 (MQ=255)
                                                                                                                                                                                                                                ggcaaaaatccaccacctctttgacTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTGAGAAACt                >  1:56398‑M2/26‑251 (MQ=255)
                                                                                                                                                                                                                                     aaaTAGGCATCTTGGTTCTGTAATGCCAGTCAGTTATGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCTCCATTCCGCATGATTCTGAGAAACTCAGtt           <  1:6620/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ccctcctctttgatCAATGCCATTTAGTTAAGCATCTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCTCTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCTTCAGCAGTCTCATCACCATTCCGCATGATTCTGAGAACCTCATTTGATTCg     <  2:176167‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                                                ttGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTGAGAAACTCAGTTGATTCGGcc  >  1:10310/1‑249 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
GCGCTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTGAGAAACTCAGTTGATTCGGCC  >  NC_000913/4501841‑4502329

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