Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,486,010	IS2 (–) +5 bp :: Δ32 bp +GTCTTAATACTAGTTTTTAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTA	intergenic (‑192/‑106)	yhfA ← / → crp	OsmC family protein YhfA/DNA‑binding transcriptional dual regulator CRP

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3186128 =	NA (NA)	11 (0.660) +GTCTTAATACTAGTTTTTAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTA	10/348	0.7	100%	noncoding (33/1331 nt)	IS2	repeat region
	?	NC_000913	3486010 =	0 (0.000)					intergenic (‑192/‑110)	yhfA/crp	OsmC family protein YhfA/DNA‑binding transcriptional dual regulator CRP
*	?	NC_000913	= 3187426	NA (NA)	27 (1.140)	24/496	0.1	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	= 3486014	0 (0.000)					intergenic (‑196/‑106)	yhfA/crp	OsmC family protein YhfA/DNA‑binding transcriptional dual regulator CRP

TTCGCCCAGAAAAGTTAACCCTTCGACCCACTTCACTCGCGCTTGCATTTTTGCTACTCCACTGCGTCAATTTTCCTGACAGAGTACGCGTACTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGT  >  NC_000913/3485771‑3486262                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         caacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTg                                                                                                                                                                                                                                                   >  1:142252‑M2/240‑251 (MQ=255)    gCCCAGAAAAGTTAACCCTTCGACCCACTTCACTCGCGCTTGCATTTTTGCTACTCCACTGCGTCAATTTTCCTGACAGAGTACGCGTACTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcc                                                                                                                                                                                                                                                >  1:210710‑M2/1‑241 (MQ=255)     gagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGAt                                                                                                                                                                                                                                                 >  1:119790‑M2/236‑249 (MQ=255)                     cTTCGACCCACTTCACTCGCGCTTGCATTTTTGCTACTCCACTGCGTCAATTTTCCTGACAGAGTACGCGTACTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacat                                                                                                                                                                                                                               <  2:52793‑M2/251‑28 (MQ=255)                       tCGACCCACTTCACTCGCGCTTGCATTTTTGCTACTCCACTGCGTCAATTTTCCTGACAGAGTACGCGTACTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatct                                                                                                                                                                                                                             <  1:149208‑M2/251‑30 (MQ=255)                          ttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGa                                                                                                                                                                                                                           >  1:116483‑M2/215‑250 (MQ=255)                          ttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGa                                                                                                                                                                                                                           >  1:180991‑M2/215‑250 (MQ=255)                             cACTTCACTCGCGCTTGCATTTTTGCTACTCCACTGCGTCAATTTTCCTGACAGAGTACGCGTACTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatc                                                                                                                                                                                                                       >  2:79420‑M2/1‑216 (MQ=255)                                                         cgccacTGCGTCAATTTTCCTGACAGAGTACGCGTTCTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGCAGCGAGACCCCAGGAGACACAAAGATAAAGCTATGCTAAATCAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCATATTACCGTGCAGTtggattggcccctagatttccagacatctgtgatcacttaacccattccaagcccgctgcctc                                                                                                                                                                                           >  2:913‑M2/3‑188 (MQ=255)                                                             gtgtggtacgccgtctgcgtttctccggccctcagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatagtcttaatactagtttttagactatcacttatttcagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGACGCATCTTTCGACAATCCAGAGACAGCGGCGTTATc                                                                                                                                                                                       >  1:235321‑M2/180‑251 (MQ=255)                                                                tggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTgg                                                                                                                                                                                    >  2:46623‑M2/177‑251 (MQ=255)                                                                    cAATTTTCCTGACAGAGTACGCGTACTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacagcggttagcacttaacccattacaagcccgctgccgcagatattccc                                                                                                                                                                                 >  1:146443‑M2/1‑177 (MQ=255)                                                                    acgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGctct                                                                                                                                                                                >  1:102982‑M2/173‑251 (MQ=255)                                                                                      aCGCGTACTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataaccca                                                                                                                                                               >  2:153361‑M2/1‑159 (MQ=255)                                                                                            aCTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcact                                                                                                                                                        <  2:235314‑M2/251‑99 (MQ=255)                                                                                              tAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactat                                                                                                                                                      <  1:85563‑M2/251‑101 (MQ=255)                                                                                                    caatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCa                                                                                                                                                <  2:11136‑M2/111‑1 (MQ=255)                                                                                                                   ggCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataat                                                                                                                                 >  2:164656‑M2/1‑130 (MQ=255)                                                                                                                       aCCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctc                                                                                                                             <  1:156910‑M2/251‑126 (MQ=255)                                                                                                                         cTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcg                                                                                                                            <  1:134622‑M2/250‑127 (MQ=255)                                                                                                                         cTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgtgataatgctcg                                                                                                                            <  1:158064‑M2/250‑127 (MQ=255)                                                                                                                           gggTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaac                                                                                                                         >  1:235314‑M2/1‑122 (MQ=255)                                                                                                                                     gAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgc                                                                                                                 >  2:170660‑M2/1‑112 (MQ=255)                                                                                                                                         tattaagactatcacttatttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTAAAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCAGACTCTCg                                                                                                            >  1:216632‑M2/104‑250 (MQ=255)                                                                                                                                                 cAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctg                                                                                                   <  2:84448‑M2/251‑152 (MQ=255)                                                                                                                                                          atttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCAtt                                                                                          >  1:87338‑M2/87‑251 (MQ=255)                                                                                                                                                            ttaagtgatagtcttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGc                                                                                        >  1:224170‑M2/85‑251 (MQ=255)                                                                                                                                                                         ttaatactagtttttagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGt                                                                           >  2:187914‑M2/72‑251 (MQ=255)                                                                                                                                                                          cTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaccgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcat                                                                          <  2:209656‑M2/251‑177 (MQ=255)                                                                                                                                                                                cAGTCAGGATGCTACAGTAATACATTGATTTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgcccttcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatact                                                                    >  2:184168‑M2/1‑69 (MQ=255)                                                                                                                                                                                 aGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactg                                                                   <  1:69443‑M2/251‑184 (MQ=255)                                                                                                                                                                                  gtttttagactatcacttatttaagtgatattggttgtctggtgattcagtgggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGCAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATcc                                                                   >  2:68589‑M2/63‑250 (MQ=255)                                                                                                                                                                                       gATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgt                                                               <  1:164518‑M2/249‑188 (MQ=255)                                                                                                                                                                                                   aTACATTTATGTACTGCATGTATGCAAAGGACGTCACATTACCGTTCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttctc                                                 <  2:65530‑M2/251‑202 (MQ=255)                                                                                                                                                                                                   aTACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGAAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatc                                                 <  1:214910‑M2/251‑202 (MQ=255)                                                                                                                                                                                                     aCATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataaggctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgcttgatcgt                                               >  2:34865‑M2/1‑48 (MQ=255)                                                                                                                                                                                                         agtgatattggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGt                                           <  1:14356‑M2/212‑1 (MQ=255)                                                                                                                                                                                                            tGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacg                                        <  2:239631‑M2/251‑211 (MQ=255)                                                                                                                                                                                                               attggttgtctggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGaaaaa                                     <  2:180991‑M2/218‑1 (MQ=255)                                                                                                                                                                                                                  gCATGTATGCAAAGGACGTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctc                                  <  1:172862‑M2/251‑217 (MQ=255)                                                                                                                                                                                                                         tggagattcagggggccagtctaGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGc                           >  2:128488‑M2/24‑251 (MQ=255)                                                                                                                                                                                                                               ggACGTCACATTCACGTTCAGGtggatttgcccctatctttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgcacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgt                     <  2:20343‑M2/251‑230 (MQ=255)                                                                                                                                                                                                                                   gTCACATTACCGTGCAGTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttact                 <  1:217583‑M2/251‑234 (MQ=255)                                                                                                                                                                                                                                                  aGTACAGTTGATATCCCCTTCCCAGGTAGCGGGAAGCATAGTTCGGCAATCCAGAGACAGCGGCTTTATCTGGCTCTGGATAAAGCTTATAACAGAGGATAACCGCGCATTGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGATCACGCTTATTCACCAGGGTGAAAAAGCGTAAACGCTTTACTACATCGTTAAAGGCTCtgt  <  1:122090/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         TTCGCCCAGAAAAGTTAACCCTTCGACCCACTTCACTCGCGCTTGCATTTTTGCTACTCCACTGCGTCAATTTTCCTGACAGAGTACGCGTACTAACCAAATCGCGCAACGGAAGGCGACCTGGGTCATGCTGAAGCGAGACACCAGGAGACACAAAGCGAAAGCTATGCTAAAACAGTCAGGATGCTACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGT  >  NC_000913/3485771‑3486262

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