Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,093,481	IS2 (+) +5 bp	coding (606‑610/1359 nt)	dgcT →	putative diguanylate cyclase DgcT

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	381260 =	NA (NA)	47 (0.580)	43/486	0.8	93.9%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 1093485	3 (0.040)					coding (610/1359 nt)	dgcT	putative diguanylate cyclase DgcT
*	?	NC_000913	1093481 =	3 (0.040)	60 (0.740)	52/490	0.5	95.1%	coding (606/1359 nt)	dgcT	putative diguanylate cyclase DgcT
	?	NC_000913	= 1651548	NA (NA)					noncoding (706/706 nt)	IS2	repeat region

GCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAAAATCTTGCCTCATCG  >  NC_000913/1651303‑1651796                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         tcGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAAAGATAAAGCGTGACTACATCAGTATCATGCCCAAAACAGACGGGTTACCGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGACTAGCGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGGGATAACAGCTGTCTGGAAATATAGGGGGAAATCCACAggg                                                                                                                                                                                                                                                     >  2:531602/2‑250 (MQ=11)      gTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACt                                                                                                                                                                                                                                                >  1:213877/1‑251 (MQ=33)        ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTAc                                                                                                                                                                                                                                              <  1:136894/251‑1 (MQ=34)         cGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACg                                                                                                                                                                                                                                             <  2:196293/251‑1 (MQ=34)                 aaCGGAATAGGAGAGAACTTCGCGAAACTGATACAGCGTGAATACATCAGGATCATGCCCCACCCAGACGGGGTAATGGCAGGAAAGAGCCTTGCAGAGGCGTTCGAGCATTATAACGAATGTCTTCCGCATAGGGCGCTGGGTTATCGCTCGCCACGTGAATATCTGCGGCAGCGGGTTTGTAATTGGCTAAGTGATAACAGATGCCTGGAAAGATAGGGGCAAAGCCACAggggtatgactctccgct                                                                                                                                                                                                                                      <  2:104375/250‑17 (MQ=2)                  aCGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCtgt                                                                                                                                                                                                                                    >  1:630448/1‑251 (MQ=39)                   cGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTTCGCTCCTCtgtt                                                                                                                                                                                                                                   <  1:636809/251‑1 (MQ=38)                     gAATAGCAGAGAGCTACTGGAAAACGATAAAGCGTGACTACAGCAGTATCATGCCCAAACCAGGCGGGTGAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCAGTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCtgtttg                                                                                                                                                                                                                                 <  2:151461/251‑1 (MQ=14)                         aGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTAc                                                                                                                                                                                                                             <  1:17887/251‑1 (MQ=255)                                agCTTCGGGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATaa                                                                                                                                                                                                                      <  2:156371/251‑1 (MQ=255)                                        gAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTgg                                                                                                                                                                                                               <  1:853743/250‑1 (MQ=255)                                          aaaCGATAAAGCGTGACTACATCAGTAGCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGgtgt                                                                                                                                                                                                            <  2:377240/251‑1 (MQ=255)                                           aaCGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGtt                                                                                                                                                                                                           >  1:322501/1‑251 (MQ=255)                                                            aCATTAGTACCATGCCCAAACCCGAAGGGTTAACGGCAGCTAAGACCCTTGCAGAGGCGTACGATCATTATAACGAATGGCAGTCGCATAGTGCGCGGGGTTATCGATCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAAGGGGTTAAGTGATAACAGATGGCTGGAAATATAGGGGCAAATCCACAAGTGAACCACGCTCCGCTGGTTGTTACGCGATAAAAACTGGTGTTTTCGACTATGCAAAACa                                                                                                                                                                                          <  2:773590/251‑1 (MQ=17)                                                                    aTCATGCCCAATCCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGTCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCcct                                                                                                                                                                                  <  2:660068/251‑1 (MQ=255)                                                                           ccAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACg                                                                                                                                                                           >  1:442048/1‑251 (MQ=255)                                                                                   gACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTgg                                                                                                                                                                   <  1:191300/251‑1 (MQ=255)                                                                                        gTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGgt                                                                                                                                                                                                              <  2:417939/203‑1 (MQ=255)                                                                                        gTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGgt                                                                                                                                                                                                              >  1:417939/1‑203 (MQ=255)                                                                                            aCGGCAGCAAAGGACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGGGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGGCTGGAAATATAGGGGCAAATCCACAGGTTAACTGCGCTCCTCTGTTTGTTACGCAAGAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTg                                                                                                                                                          <  1:113701/251‑1 (MQ=255)                                                                                            aCGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTg                                                                                                                                                          <  2:164642/251‑1 (MQ=255)                                                                                                    aaaGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGGAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTTTTTTCGACTAGGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAAGTTGGTGAGGcc                                                                                                                                                  <  1:757406/251‑1 (MQ=255)                                                                                                              gCAGAGGCGGTCGAGCATTATAACGAATGGCATCCGCATAGTTCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGAtt                                                                                                                                        <  2:705908/251‑1 (MQ=255)                                                                                                                                tacaaCGAATGGCATCCGCATAGTGCGCTGGGTTATCGCGCGCCACGGGAAGATCTGCGGCAGCGGGCTGGTAATGGGTTAAGTGATAACAGATGTCTGGAAAGATAGGGGAAAAGCCTCAGGTGAACTACGCGCCTCTGGTTGTTACGCAATAAAAACGTGTGTGTTCGACTATGCAAAACAGTTGCACGCCTCACGCAAACTGGAAAAAGTTGGTGAGGACCGACAATATTTAACCATAAAAGAACTg                                                                                                                       <  2:71553/247‑1 (MQ=21)                                                                                                                                                      gTGCGCTGGGTGATCGCAGGCTACGGGAATATCTGCGGCAGCGGGGTTGTAATGGGTTAAGTGATAACAGATGGCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCACCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAAATTGAACCATAAAAGAACTTGCTGAAAAATTTCTGGCTCTgaa                                                                                                <  2:39951/251‑1 (MQ=255)                                                                                                                                                            ggggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAATGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAAATACGCTTCCCTGTTTGTTACGCAATAAAAACCGGTGTTTTCGACTATGCCAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAAGCCCGACCAGATTTAACCAAAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAAc                                                                                          <  1:769757/250‑1 (MQ=37)                                                                                                                                                             gggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACt                                                                                         <  2:637799/251‑1 (MQ=255)                                                                                                                                                              ggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTg                                                                                        >  1:756253/1‑251 (MQ=255)                                                                                                                                                                    cgctcgcCACGTGAATATCTGCGGCAGCGGGCTTGTAAGTGGTTAATTGAGAACAGATGTCTGTAATGATAGGGGCAAATCCACAGGGTAACTGCGCTCCTCGTTTTGTTTCGCAATAAAAACGGGTGTTTTCGTCTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAAGTTGGTGAGGCCCGACAAGAGTTAACCATAAAAGAACTGGCTGAAAAATTGCTGGCGCTGAAAGAAACTGAAGta                                                                                   <  2:529746/250‑2 (MQ=255)                                                                                                                                                                           cACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCCTTGGTTGTTACGCAATAAAAACTGGTGTTTTCGAATAAGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATATAACCATAAAAGAACTGGCTGAAAAATTTCTGTCACTGAAAGAAACTGAAGTCGCCAAAAc                                                                           <  2:213877/251‑1 (MQ=38)                                                                                                                                                                                    atCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTTAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCa                                                                   <  2:442048/250‑1 (MQ=255)                                                                                                                                                                                         cGGCAGCGGGTGTGTAATGGGTTAAGTTATAATAGTTGTCTGGATATATAGTGGCAAATCTATAGGTGAACTACGCACCTCTGTTGTTTACGCATTAAAAACTGGTGTCTTCCACTATGCAATACAGTTTACCTCCTCACGCAATCTTGAAAAATTTGGGGAGGCCCGCAAAGAATTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACGGAAAGAAACTGAAGTCGCTAAAACATCACTcttctcat                                                             <  2:613477/251‑8 (MQ=255)                                                                                                                                                                                              gCGGGCTTGTAATGGGTTAAGTGATCACAGATGTCTTGAAATATAGGGTCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGGGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGGCGCCAAAACATCACTCAACACAGAcc                                                          <  2:590910/249‑1 (MQ=255)                                                                                                                                                                                                                aaGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCGCAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATc                                      <  1:321144/251‑1 (MQ=255)                                                                                                                                                                                                                 aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCATTCTGGAAAAATTTGGTGAGGCCCGACAAGATTTa                                                                                                                                      <  2:266996/154‑1 (MQ=255)                                                                                                                                                                                                                 aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCATTCTGGAAAAATTTGGTGAGGCCCGACAAGATTTa                                                                                                                                      >  1:266996/1‑154 (MQ=255)                                                                                                                                                                                                                   tGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTg                                   >  1:207138/1‑251 (MQ=255)                                                                                                                                                                                                                   tGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTAAACTACGCACCTCTGGTTGTTACGCGATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCACGACAAGATTTAACCAGAAAAGAACAGGATGAAAAATTTCTGGCACTGAAAGAAACTTAAGTCACCAAAACATCACTCAACAAAAACCGTGCCGTCATCTAAAATCTccgg                                   >  2:39382/1‑249 (MQ=255)                                                                                                                                                                                                                          aGATGGCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATaa                            <  2:630448/251‑1 (MQ=255)                                                                                                                                                                                                                            aggTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATc                          <  2:473538/249‑1 (MQ=255)                                                                                                                                                                                                                                   gAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGATCATAATCGGTGaaa                   <  2:322501/251‑1 (MQ=255)                                                                                                                                                                                                                                        ataGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAAAAt              >  1:23892/1‑251 (MQ=255)                                                                                                                                                                                                                                                   tCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAAAATCTTGCCTCATc   >  1:433818/1‑251 (MQ=255)                                                                                                                                                                                                                                                    ccACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAAAATCTTGCCTCATCg  <  1:726607/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         GCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGAAAAAAATCTTGCCTCATCG  >  NC_000913/1651303‑1651796

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.