Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,542,042	IS1 (+) +8 bp	coding (6‑13/597 nt)	fimE →	regulator for fimA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	25 (0.660)	22/486	0.6	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 4542049	0 (0.000)					coding (13/597 nt)	fimE	regulator for fimA
*	?	NC_000913	= 3584195	NA (NA)	53 (1.400)	46/488	0.1	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	4542042 =	0 (0.000)					coding (6/597 nt)	fimE	regulator for fimA

GTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGA  >  NC_000913/3583946‑3584443                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         gTCACACTTTGACGTGGAGCAATAGATGACGGATGGCTGGCCGCTGTATGACTCCCGCCTGAAGGGAAAGCTGCACGTAATAAGCAAGCGATATACGCAGCGACTTGAGCGGCATACCCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCACAATCGGTGGAGCTGCATGACAAAGTCATCGGGCCTTATCTGAACATACAACACTATCAATAAGTTGGAGTCATTACCg                                                                                                                                                                                                                                                         >  2:238566/1‑251 (MQ=11)  tCACCCTTTGACGTGGTGATATGGATGACGGATTGCTTGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGTCAAAGTCATCGGTCATTATCTGAACATAAAACACTATCAATAAGTTTGAGTCATTACCg                                                                                                                                                                                                                                                         <  1:350615/250‑1 (MQ=2)      ccTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAACTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACCTAAAACAATATCAATAAGTTGGGGTCATTCCCGACAtt                                                                                                                                                                                                                                                    >  1:391118/1‑251 (MQ=255)             cGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTa                                                                                                                                                                                                                                             <  1:198836/251‑1 (MQ=25)              gtggtgATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                                                                                                                                                                                                                                            <  1:182232/251‑1 (MQ=25)              gtggtgATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTa                                                                                                                                                                                                                                             >  1:130575/1‑250 (MQ=25)                 gtgATATAGAGGACGGATGGCTGTCCGCTGTATGAATCCCGCCTGAAGTGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCTCCTGGTACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTTGAGCTGCATGAAAAAGTCATCGGGCATTATCTGAACATAAAACTCGAGCAATAAGTTTGAGGCATGACCAACATTAAAAGTAACtt                                                                                                                                                                                                                                         <  2:188432/251‑1 (MQ=11)                  ggATATAGATGACGGATGGCTGGCCGCTGTCTGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                                                                                                                                                                                                                                        <  1:278459/250‑1 (MQ=255)                    atatAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                                                                                                                                                                                                                      <  1:234223/251‑1 (MQ=38)                      atagATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                                                                                                                                                                                                                                     <  1:96976/247‑1 (MQ=39)                          aggACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAt                                                                                                                                                                                                                                <  1:54074/249‑1 (MQ=255)                          aTGACGGATGGCTGGCCGCTTTATGAAGCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCTATATACGCAGCGAATTGAGCGGCATAACCGGAATCTGATGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCTTTCTCAAAATCGGTTGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGGCATTACCGACATTAAAAGTAACTTTTTTGAAAt                                                                                                                                                                                                                                <  1:279984/251‑1 (MQ=14)                                 aTGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGAAATCAGCACGCGATATACGCAGCGCATGGAGCGGCATAATCTGAATTTAAGGCAGCGCCTGGCGCGGCTGGGGCGGAAGGCGCTGGCTCTCTCAAAATCGGAGGAGCTGCAGGACGAAGTCATCGGGCATTATCTAAACTTAAAGCACTATCAATAAGTTGGAGTCAGTACCGACATAACAAGTTACTTTTTGGAAATTAGAGat                                                                                                                                                                                                                         <  2:396309/251‑1 (MQ=2)                                                aTGAATCCCGCCTGAATGGAAAGCTGCACGTAATCAGCAATCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGTTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATGATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTGTGAAATTAGAGATAATAGTGGCACAGa                                                                                                                                                                                                           <  2:338497/250‑1 (MQ=255)                                                           cTGAAGGGAAAGCTACACGAAAACCGCACGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATATGAGGCAGCACCTGGAACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGACATCGGGCATTATCTGAAAATAAAACACTATCAATAAGTTGGCGTCATTACCGACATTCAAAGTAACTTTTTTGCAATTAGAGATAATATTGGCCAAGAACATATATTc                                                                                                                                                                                                 >  2:351302/1‑249 (MQ=14)                                                              aaGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTc                                                                                                                                                                                            <  1:113500/251‑1 (MQ=255)                                                                         gCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCGGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGAt                                                                                                                                                                                  <  1:237808/250‑1 (MQ=255)                                                                                   aGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCCGCACCTGGCACGGCTGGGAAGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCAtt                                                                                                                                                                       >  1:370306/1‑251 (MQ=255)                                                                                   aGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGtt                                                                                                                                                                                         >  1:232695/1‑233 (MQ=255)                                                                                   aGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAGTGGCACAGAAAATATATTCAGTTCGtt                                                                                                                                                                                         <  2:232695/233‑1 (MQ=255)                                                                                       aGCGATATCCGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACTGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAATTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGATATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTa                                                                                                                                                                   <  2:232943/251‑1 (MQ=255)                                                                                          gATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAAGATATTCAGTTCGTGTAGTGATTTCTCTTCATTCTTAAAc                                                                                                                                                                <  2:130575/251‑1 (MQ=255)                                                                                                                gCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCat                                                                                                                                          >  1:93168/1‑251 (MQ=255)                                                                                                                         gAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATc                                                                                                                                   <  1:239953/249‑1 (MQ=255)                                                                                                                                     gcaCCTGGCACGGCTGGGACGGAAGTCGATGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACaaaa                                                                                                                     <  2:328447/251‑1 (MQ=255)                                                                                                                                       aCCTGGCCCGGCTGGGACGGAAGTCGCTGTCTTTCTCAAAAGCGTAGGATCTGCATGACAAAGACATCGGGCACTATCTGAACAGAACACACTATCAATAAATTGGAGCAATTTCCGACATTAAAATTAACTTTTTTGAAAGTAGAGAGAATATAGGCTCACAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATACTTAAACGAAATCACGGACCGTTTTTCATAATAATCAAAAgccaataaaaaa                                                                                                                   <  2:272560/251‑12 (MQ=255)                                                                                                                                            gCACGGATGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCAtt                                                                                                               <  1:236841/250‑1 (MQ=255)                                                                                                                                              aCGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCAtt                                                                                                                                                                       >  2:155868/1‑192 (MQ=255)                                                                                                                                              aCGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTCCGTTTAGTGATTTCTCTTCAtt                                                                                                                                                                       <  1:155868/192‑1 (MQ=255)                                                                                                                                                         ggaatgcggtGTCGTTCTCAAAATCGGTGGAGCTGCAGGACAAAGTCATCGGGGATTATCTGAAGATAAAACACTTTCAATAAGTTGGAGTCATTACCGACATTAAAAGTATCTTTTGTGAAATTAGAGGTAATATTTTCACAGAAAAGATATTCAGTTCTTTTAGTGATTTCTCTTCAGTCTTAAAGGAAGGCATGGACAGTTGTTCATAGTAATCAAAAGGAAACAAAAAGCATTGtggtcgaggccat                                                                                                 <  2:174820/242‑13 (MQ=255)                                                                                                                                                                            aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACcagc                                                                              >  1:332259/1‑251 (MQ=255)                                                                                                                                                                                                aCAAAGTCATCGGGCATTATCTGAACAGAAAGCACTATCAATAAGTTGGAGTCATTACCGACATTAAAATTATCTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATGCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACTAAATCATGGTCAGTTGTTCATAATAATCAAAATGAAACAAAAATCATTGGAAACCAAGCCATCTTGTAATAGATTCACAGCAGCAAATTCCTTGTTCGaaa                                                          <  2:20690/251‑1 (MQ=255)                                                                                                                                                                                                           gggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGc                                               >  1:271396/1‑251 (MQ=255)                                                                                                                                                                                                              cATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACcagcag                                                                            <  2:125726/219‑1 (MQ=255)                                                                                                                                                                                                              cATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACcagcag                                                                            >  1:125726/1‑219 (MQ=255)                                                                                                                                                                                                                         cATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGaaaa                                  >  1:74853/1‑250 (MQ=255)                                                                                                                                                                                                                               acacTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTAc                            >  1:96489/1‑250 (MQ=255)                                                                                                                                                                                                                                           aggtggcgtcgttACCGACATTAAAAGTAACTTTTTTGAAATTAGCGATAATATTGGCCAAGAAAATATATTCTGTTAGTTTAGTAATTTCTATTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGCAACAACAATCATTGTAAACACTGCCATCTTTTTAGAGATTACAAGCAAAAAATTACTTGTTCGAAATTTAGCAGTGAAAGCTATGGAAAAACTTACAACGGAACTacac               >  2:271342/12‑251 (MQ=255)                                                                                                                                                                                                                                                        aCCGACATTAAAAGTAACTTTGTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTgaga  <  1:238566/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         GTCACCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGA  >  NC_000913/3583946‑3584443

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.