Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,542,930	A→G	intergenic (+297/‑185)	fimE → / → fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	4,542,930	0	A→G	100.0%	212.5	55	intergenic (+297/‑185)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
Reads supporting (aligned to +/- strand):  new base (51/4):  ref base (0/0):  total (51/4)

ACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAA  >  NC_000913/4542766‑4543121                                                                                                                                                                     |                                                                                                                                                                                                aCTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGAt                                                                                                                                                              >  4:252394/1‑200 (MQ=255)                                           caTCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAATAATAGAATATTGAGCAAACAAATAAACTGa                                                                                                                                                                               >  1:535865/1‑141 (MQ=255)                                                ccTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCcgc                                                                                                                                                                 >  3:713331/1‑150 (MQ=255)                                                   ccGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGttt                                                                                                                                                                       >  1:493096/1‑141 (MQ=255)                                                         tatGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtg                                                                                                      >  4:109733/1‑200 (MQ=255)                                                          atGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGt                                                                                                                                                                         >  2:561704/1‑132 (MQ=255)                                                          atGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTcc                                                                                                                                                       >  3:321259/1‑150 (MQ=255)                                                                        tttCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCAGCGATGCTTTCCTCTATGAGTCaaaa                                                                                                                                         >  3:614815/1‑150 (MQ=255)                                                                                        cATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaa                                                                       >  4:601854/1‑200 (MQ=255)                                                                                             tCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCTCGATGCTTTCCTCTATGAGTCAAAATGGCACCAAATGTTTCATCtt                                                                                                                    >  3:88452/1‑150 (MQ=255)                                                                                             tCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATgg                                                                                                                                      >  2:437580/1‑132 (MQ=255)                                                                                                tCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                                                                                          >  1:542312/1‑141 (MQ=255)                                                                                                tCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTg                                                                                                                 >  3:685561/1‑150 (MQ=255)                                                                                                 cATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGcccc                                                                                                                                  >  2:81705/1‑132 (MQ=255)                                                                                                          ctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACt                                                                                                       >  3:189978/1‑150 (MQ=255)                                                                                                          ctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAAc                                                       <  4:189983/198‑1 (MQ=255)                                                                                                                     accTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                                                                                            >  3:63707/1‑150 (MQ=255)                                                                                                                           tGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                                                                                        <  2:143834/132‑1 (MQ=255)                                                                                                                               atatAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                                                                    >  2:326258/1‑132 (MQ=255)                                                                                                                               atatAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                                                                           >  1:22636/1‑141 (MQ=255)                                                                                                                               atatAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                                                                           >  1:693696/1‑141 (MQ=255)                                                                                                                               atatAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc                                                                                  >  3:333267/1‑150 (MQ=255)                                                                                                                               atatAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc                                                                                  >  3:21931/1‑150 (MQ=255)                                                                                                                               atatAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc                                                                                  >  3:15351/1‑150 (MQ=255)                                                                                                                               atatAAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGAAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGaaa                                >  4:380494/1‑200 (MQ=255)                                                                                                                                 ataAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATa                              >  4:396285/1‑200 (MQ=255)                                                                                                                                 ataAACATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCACAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCa                                                                                         >  1:593625/1‑141 (MQ=255)                                                                                                                                      cATCTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaa                                                                                    >  1:411844/1‑141 (MQ=255)                                                                                                                                         cTATAAATAAAGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCATACTCGTTTacaa                                                                        >  3:193961/1‑150 (MQ=255)                                                                                                                                                 aaaGATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCTGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCt                                                                                                                     >  1:327830/1‑97 (MQ=255)                                                                                                                                                    gATAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                                                                      >  1:354107/1‑141 (MQ=255)                                                                                                                                                     aTAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc                                                            <  3:317925/150‑1 (MQ=255)                                                                                                                                                      tAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGTAAAACTGTTCAGTGTTTGCAGTCAAACTCGTTTacaaaaca                                                                    >  1:458268/1‑141 (MQ=255)                                                                                                                                                      tAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacc                                                                    >  1:695646/1‑140 (MQ=255)                                                                                                                                                      tAACAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAagcagc         >  4:498277/1‑200 (MQ=255)                                                                                                                                                         naaTAGAATANNNNNNNAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACANNNNGTACAGaa                                                        >  3:135018/2‑150 (MQ=255)                                                                                                                                                         cAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTac                                                                          >  2:660520/1‑132 (MQ=255)                                                                                                                                                         cAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                                                                 >  1:621767/1‑141 (MQ=255)                                                                                                                                                         cAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                        >  3:514092/1‑150 (MQ=255)                                                                                                                                                         cAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATg      >  4:285400/1‑200 (MQ=255)                                                                                                                                                         cAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATg      >  4:665405/1‑200 (MQ=255)                                                                                                                                                         cAATAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTGTACAAAACAAAGTGTACAGaa                                                        >  3:92983/1‑150 (MQ=255)                                                                                                                                                           aTAGAATATTGAGCCAACAAATAAACTGCAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacta                                                                        >  2:171742/1‑130 (MQ=255)                                                                                                                                                           aTAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGCTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAATTGTGCAGCGTTGGCAGTCAATCTCGTTTACAAAACAAAGTGTACAGACCg                                                      >  3:176189/1‑150 (MQ=255)                                                                                                                                                           aTAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaa                                                                        >  2:239098/1‑132 (MQ=255)                                                                                                                                                           aTAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaa                                                                        >  2:185470/1‑132 (MQ=255)                                                                                                                                                           aTAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACg                                                      >  3:317735/1‑150 (MQ=255)                                                                                                                                                           aTAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACg                                                      >  3:223611/1‑150 (MQ=255)                                                                                                                                                           aTAGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCGTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCTAACTCGTTTacaa                                                                        >  2:258409/1‑132 (MQ=255)                                                                                                                                                             aGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTa                                                             >  1:17460/1‑141 (MQ=255)                                                                                                                                                             aGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGAc                                                    <  3:252384/150‑1 (MQ=255)                                                                                                                                                             aGAATATTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGaaaa  >  4:541705/1‑200 (MQ=255)                                                                                                                                                                atatTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAACTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa                                                                   >  2:119754/1‑132 (MQ=255)                                                                                                                                                                atatTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacacaacaa                                                                   >  2:655004/1‑132 (MQ=255)                                                                                                                                                                atatTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAg                                                          >  1:543389/1‑141 (MQ=255)                                                                                                                                                                atatTGAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGc                                                 >  3:855717/1‑150 (MQ=255)                                                                                                                                                                     |                                                                                                                                                                                                ACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAA  >  NC_000913/4542766‑4543121

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