Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	789,812	IS2 (–) +5 bp	coding (164‑168/1149 nt)	galK ←	galactokinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	381260 =	NA (NA)	42 (0.970)	38/486	0.2	93.2%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	789812 =	3 (0.070)					coding (168/1149 nt)	galK	galactokinase
*	?	NC_000913	= 382590	NA (NA)	51 (1.190)	43/482	0.1	96.2%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	= 789816	2 (0.050)					coding (164/1149 nt)	galK	galactokinase

GCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCG  >  NC_000913/382345‑382834                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     gCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTAGCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTAc                                                                                                                                                                                                                                                 <  1:249764/251‑1 (MQ=14)  cGGATTCCGGAGAGGAACGTAATAGCAGATAGCTTCGTGAAAACGATAAAGCGTTACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTAcc                                                                                                                                                                                                                                                <  2:313250/251‑1 (MQ=9)     aGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGATAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTAt                                                                                                                                                                                                                                             >  1:445510/1‑251 (MQ=17)        ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAg                                                                                                                                                                                                                                          <  2:410166/251‑1 (MQ=33)            agagTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGc                                                                                                                                                                                                                                        >  1:169619/1‑249 (MQ=35)               gTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGttt                                                                                                                                                                                                                                   <  2:445510/251‑1 (MQ=37)                      aaTAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCa                                                                                                                                                                                                                             >  1:121501/1‑250 (MQ=40)                       aTAGCAGAGAGCTTCGTGAAAACGATAAAGCGTCACTACATCAGTATCATGCCCACACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCAACCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATACCAGAGGTGTGGAAATATAGGGGCAAATCAAATTACCAATCAGGCAATTTGCGCCAAt                                                                                                                                                                                                                           >  2:29256/1‑251 (MQ=11)                             gagagCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACAAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTAGAGCATTATAACGAATGGCATCCGCGTAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGGGATACCAGATGTCTGTACATGTAGGGGCAAATCCAATTATCTAACGGGCAGTTTGCTACAATCACTAc                                                                                                                                                                                                                     >  2:461450/1‑251 (MQ=11)                              agagCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCCGTTTGCGCCAATCACtccc                                                                                                                                                                                                                    >  2:381630/1‑248 (MQ=255)                                     cgggAAAACGATAAAGCGTGCCTACAGCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTACGAGCATTATAACGAATGGCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCt                                                                                                                                                                                                             <  2:162623/248‑1 (MQ=255)                                                 aaaGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGCACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCAACGGGAATATCTGCGACAGCGGGCTTGTAATGGGTTAAGTGATAACAGCTGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCCGTTTGCGCCACTCACTACCGGACCATGAATAACGGCAc                                                                                                                                                                                                 >  2:49854/1‑251 (MQ=14)                                                            aCATCAGTATCATGCCCAAACCAGACGGTTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCAGTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATATCAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                                                                                                                                                                        <  1:45377/249‑1 (MQ=255)                                                                cAGTAGCATGCCCAAACCAGACGTGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACgg                                                                                                                                                                                   <  1:430445/250‑1 (MQ=255)                                                                           ccAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCCGAGGCGTTCGAGAATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGTACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCg                                                                                                                                                                       >  1:140980/1‑251 (MQ=255)                                                                           ccAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCg                                                                                                                                                                       <  1:181185/251‑1 (MQ=255)                                                                             aaaCCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGAt                                                                                                                                                                     <  1:352765/251‑1 (MQ=255)                                                                               aCCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGAGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGTCAAATACGGGAATTCGCGCCGATa                                                                                                                                                                    >  2:102951/1‑250 (MQ=255)                                                                                 cAGACGGGGTAACGTCAGCAAATAACCTTGTAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGTTTAGCGCTCGCCACGTGAATATCTGCGGCAGCGGTCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTTCGCCAATCACTACCGGACCCTGAATAACGGCACCATTGCAAATACGGGTATTCGCGCCGATaaca                                                                                                                                                                 <  1:418891/251‑1 (MQ=18)                                                                                      gggTTAAGGGCAGAAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGAGGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCAGATAACAACTTc                                                                                                                                                            <  2:275935/251‑1 (MQ=255)                                                                                                              gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGAGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACACTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAc                                                                                                                                    >  1:434863/1‑251 (MQ=255)                                                                                                                aggggtgTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACTGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCCCCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTATCACCAATAATTACCt                                                                                                                                  <  2:261623/245‑1 (MQ=255)                                                                                                                          gtgCATTATAACGAATGGCATCCGCATAGTGCGCTGGGGTATCGCTCGCCACGGGAAAATCTGCGGCAGCGGGCTTTTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAAt                                                                                                                        <  2:427409/249‑1 (MQ=255)                                                                                                                           aGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATc                                                                                                                       >  1:294908/1‑251 (MQ=255)                                                                                                                                    aCGAATGGCCTCCGCATCGTGCTCTGGGTTATCGTGCGCCACGGGAATTTCTGCGGAAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAGTCTCTCCCTGACCCTGAATAACGGCCCCATGCCAAATACGGGTATTCGCGCCGATAACATCTTCACCTGCAGTTTCATCAATAATTTCCTGTTCTGCAATCTGActatgt                                                                                                              <  2:130646/251‑3 (MQ=21)                                                                                                                                        aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGGATATCTGAGGCAGCGGGCTTGTAATGGGTTAAGCGATCACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGACTAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACACCTTCACCTGCAGTATCATCAATAATTACCTGGTCTGCAATCTGAATAttttttt                                                                                                          >  2:299602/1‑249 (MQ=37)                                                                                                                                               ccGCATAGGGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTGGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAAGAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTaa                                                                                                    <  2:172994/250‑1 (MQ=255)                                                                                                                                                     aGTGCGCTGGGTTATCGCTCGACACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTa                                                                                               >  1:214697/1‑249 (MQ=255)                                                                                                                                                             gggTTATGGCCCGCCACTGGAATTTCTTCGGCAGTGGGCTTGTAAAGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTTCCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTTCCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGa                                                                                      <  2:102298/250‑1 (MQ=38)                                                                                                                                                              ggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAAAATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGACAATCAATACCGGACCCTGAATAAAGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATAAATAATTACCTGGTCTGCAATATGAATATTTTTCTCTTTAACATTAAGCTTCAGAt                                                                                     >  1:29755/1‑250 (MQ=255)                                                                                                                                                                                        gCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGaa                                                             <  2:294908/248‑1 (MQ=255)                                                                                                                                                                                                    ttGTAATGGGTTAAGTGATAACAGATTTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTTCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTGTTTCAGTTCCTGTGACAGAAGGTTGGGaa                                                    <  1:356057/245‑1 (MQ=255)                                                                                                                                                                                                    ttGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGTGCAATTCCATTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGACACCATGACAAATCCGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGACTATTTTTCTCTTTAACATTAAGCTTCAGCTAGTTTTTCAGTACCTGTGACAGAAGGTTTGGaa                                                    >  2:356057/1‑245 (MQ=255)                                                                                                                                                                                                     tggaaggggTTAAGTGATAACAGATGTCTGGACATCTAGGGGCCAATCCAATTACCTATCAGGGACTTTCCGCCAATCACTACCGGACCCTCAATAACGGCACCATGTCAAATACGTGTCTTCGCGCCGATAACAACTTCCACTGCTGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTGCTCTTTAACATTAAGCTGCAGATAGTTGTGCAGTTCCTGTGACAGCAGGTTGGGAATAAGCg                                              <  2:105385/244‑1 (MQ=255)                                                                                                                                                                                                          tGGGTTAAGTGATAACAGATGTCTGGAAATCTAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGATTGGGAAAAAGCGGTTGAc                                        >  2:23979/1‑251 (MQ=255)                                                                                                                                                                                                                          aGATGTCTGGAAATATAGGGGCAAATACACTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCgg                        >  1:364605/1‑251 (MQ=255)                                                                                                                                                                                                                                        atatgggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTGCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtggg           <  2:431230/246‑1 (MQ=255)                                                                                                                                                                                                                                                aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTcgag  <  1:199937/251‑3 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     GCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCG  >  NC_000913/382345‑382834

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