Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,320,770	IS2 (+) +5 bp	coding (728‑732/1827 nt)	atoS →	sensor histidine kinase AtoS

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2068941 =	NA (NA)	37 (0.620)	31/480	0.8	97.3%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	2320770 =	1 (0.020)					coding (728/1827 nt)	atoS	sensor histidine kinase AtoS
*	?	NC_000913	= 2070271	NA (NA)	49 (0.810)	45/482	0.3	98.0%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 2320774	1 (0.020)					coding (732/1827 nt)	atoS	sensor histidine kinase AtoS

ACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACC  >  NC_000913/2068707‑2069189                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         aCTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCGTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTTGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTatat                                                                                                                                                                                                                                          <  1:605176/251‑1 (MQ=255)         cGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGAc                                                                                                                                                                                                                                  >  1:119656/1‑251 (MQ=255)                     aGAAACAAAAATCAGCAACGAAGAGACGATCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCAc                                                                                                                                                                                                                      >  1:588535/1‑251 (MQ=255)                      gAAACAAAAATCAGCAACGAAGAGACGCTAGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACACCTCCTTGGCAGGCCGTAAACAAATCAGCATACGTGGCGGGATTTGCCCCTATATTTCAAGACATCTGTTATCACt                                                                                                                                                                                                                     >  1:656424/1‑251 (MQ=255)                       aaaCAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAACGGTTAGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAAGGCATCCATGGCCGATAATTTGCAGGACATACCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACtt                                                                                                                                                                                                                    <  2:15937/251‑1 (MQ=255)                          cAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGCACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGACTGGATTGGCCCCTATATTTCCAGACCTCTATTATCACTTAAc                                                                                                                                                                                                                 >  2:306791/1‑251 (MQ=255)                           aaaaaTCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                                                                                                                                                <  1:87997/251‑1 (MQ=255)                             aaaTCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAAccc                                                                                                                                                                                                               >  1:344222/1‑250 (MQ=255)                                                                       cGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                                                                                                                                                                     <  1:341303/250‑1 (MQ=255)                                                                           acgcaaCCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAAACATCTGTTATCCCTTacc                                                                                                                                                                                                                 >  2:590593/3‑200 (MQ=255)                                                                           aagcaaCCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAAc                                                                                                                                                                                                                 <  1:590593/202‑1 (MQ=255)                                                                           aagcaaCCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCAtt                                                                                                                                                                                                            >  1:223753/1‑207 (MQ=255)                                                                           aagcaaCCGTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCAtt                                                                                                                                                                                                            <  2:223753/207‑1 (MQ=255)                                                                                     tGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCAATTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATg                                                                                                                                                      >  1:328107/1‑251 (MQ=255)                                                                                              ccGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                                                                                                                                              >  1:650100/1‑250 (MQ=255)                                                                                                     cTATACAGGTTATTGACCAACAAACATTGCAGCACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGtt                                                                                                                                       >  1:603521/1‑250 (MQ=255)                                                                                                      tataCAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACCAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTtat                                                                                                                                     >  1:389210/1‑251 (MQ=255)                                                                                                                         ccaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGAAGCCATTCGTTATAATGCTCGAACGCCTCTGc                                                                                                                  >  2:131940/3‑251 (MQ=255)                                                                                                                         caaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCCCTGc                                                                                                                  >  2:366067/1‑251 (MQ=255)                                                                                                                          aaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAACCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGc                                                                                                                  >  1:95582/1‑250 (MQ=255)                                                                                                                                   aGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCttt                                                                                                        <  1:366067/251‑1 (MQ=255)                                                                                                                                     aaCTCCACCACGGCGTCCATAGCCGCTAATTTGCAGGACATCCCCGGAGTAAAGCTAACAGACAACTCCTTGGCAGGCCGTCAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCCTTCCTAGCCCGCTGCCGCAGATATTCCCGTGGCCAGTGATCACCCAGCGCACTAGGCGGATGCCATTCGTTATAATGCTCGTACGCCACTGCAAGGTTCttacc                                                                                                      >  2:94352/1‑248 (MQ=18)                                                                                                                                     aaCTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTtgc                                                                                                      <  1:136118/251‑1 (MQ=255)                                                                                                                                           accaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACCCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATGCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGt                                                                                                <  1:602378/251‑1 (MQ=255)                                                                                                                                                  ccTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCCTTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAAACCg                                                                                         >  1:294719/1‑251 (MQ=255)                                                                                                                                                  ccTCCATAGCCGATAATTTGCAGGACAACCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACGTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCg                                                                                         <  2:233711/251‑1 (MQ=255)                                                                                                                                                          gCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGttt                                                                                 >  2:4476/1‑251 (MQ=255)                                                                                                                                                              aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTCTGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGc                                                                             >  2:320955/1‑251 (MQ=255)                                                                                                                                                                         ttcatccccgtcggagagATAACAGACACTTCCGGTGCAGGCCGTAAACAAATCCGCATGCGTGGCTGGATGGGCCCCTAGATTTCCAGACATCTGGTATCACTGACCCCGTTACAAGCCCGCGGCCGCAGAGATTCCGGGGGCGAGCGATAACCCAGAGCACTATGCTGATGCCATTCGTTATAATGCTCGAACGCCTCTTCAAGGTTCTTGGCTGCCGTTAACCCGTCTGGTTTGGGCATgtttctga                                                                   <  2:51690/236‑8 (MQ=9)                                                                                                                                                                               cccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCTAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGTTGTAGTc                                                            >  1:280794/1‑251 (MQ=255)                                                                                                                                                                               cccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTc                                                            <  1:654280/251‑1 (MQ=255)                                                                                                                                                                                ccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTc                                                            >  1:240953/1‑250 (MQ=255)                                                                                                                                                                                 cGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCAc                                                          >  1:280652/1‑251 (MQ=255)                                                                                                                                                                                       atagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATGCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTa                                                    <  1:396149/249‑1 (MQ=255)                                                                                                                                                                                            aaaaagaaaaccccTTGGCAGGCCGTAAACAAATCCTCAGTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCGCTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGtt                                               <  2:668356/239‑1 (MQ=255)                                                                                                                                                                                             aaCAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGtt                                               <  1:131940/250‑1 (MQ=255)                                                                                                                                                                                                           tGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCCTGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGctctct                                >  2:187034/1‑251 (MQ=255)                                                                                                                                                                                                           tGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGctctcc                                >  2:253195/1‑250 (MQ=255)                                                                                                                                                                                                            ggCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTg                               >  1:196385/1‑251 (MQ=255)                                                                                                                                                                                                                                     cGTGGCTGGCTTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATACCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTACACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCg      >  1:632589/1‑251 (MQ=14)                                                                                                                                                                                                                                     cGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCg      >  1:2937/1‑251 (MQ=31)                                                                                                                                                                                                                                       tGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGc     >  1:260362/1‑250 (MQ=31)                                                                                                                                                                                                                                         gCTGGTTTTGCCCCTATAGTTCCAGACAGCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGTGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCAcc  <  2:307258/251‑1 (MQ=2)                                                                                                                                                                                                                                         gCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGcccc  >  1:482386/1‑248 (MQ=11)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         ACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACC  >  NC_000913/2068707‑2069189

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