Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,365,526	IS2 (–) +5 bp	coding (1244‑1248/1302 nt)	dcuA ←	C4‑dicarboxylate transporter DcuA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3186096 =	NA (NA)	27 (0.930)	25/488	0.2	86.9%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	4365526 =	4 (0.140)					coding (1248/1302 nt)	dcuA	C4‑dicarboxylate transporter DcuA
*	?	NC_000913	= 3187426	NA (NA)	24 (0.830)	21/486	0.4	85.5%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	= 4365530	4 (0.140)					coding (1244/1302 nt)	dcuA	C4‑dicarboxylate transporter DcuA

GTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATT  >  NC_000913/3187179‑3187668                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    gtgcgggttCCGGAGAGTAACGGAATAGCAGAGAGCGTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACTGGTTAACGGCAGCAAAGAACCGTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCGCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGGGATAACAGATGTCTGTAAATATAGTGGCAAAGCCAGtt                                                                                                                                                                                                                                                  <  2:51378/243‑1 (MQ=11)   gCGGAGTCCGGAGAGTAACGGAATAGCAGAGATCTTCGTGAAATCGATAAAGCTTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTat                                                                                                                                                                                                                                                <  1:372319/251‑1 (MQ=14)                                agagCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCg                                                                                                                                                                                                                  <  1:136368/251‑1 (MQ=255)                                                taTAAAGCGTGACTACATAAGTATCATGCCCAAAACAGACGGGTTAACGGAAGCAAAGAACCTTGCAGAGGCGTTCAAGCATTATAACGAATGGCATCCGCAAAGTGCGCTGGGTTATCGCTCGCCACGGGAAAATCATCGACAGCGGGATTGACATGGATTAAGTGATAACAGATGTCTGGAAATCTAGGGGCAAATCCAGTTATCATCATTGGATGT‑GCGTCCGCTTTGGCCGTTGAAtgacacaata                                                                                                                                                                                                   >  2:254881/2‑241 (MQ=14)                                                  tAAAGCGTTACTACATCAATATCAAGCACAAACCCGACGGGTTCACGGCAGCAAAACAACTTGCAGCGGCGATCGAGCATTATCACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGTGCATATCTGAGGCAGCTGACTTGTAATGGGTTAAGAGATCACAGATATCTGGCACTATAGGTGCAAATCCAGTTATCATCATAAGTTAT‑GCGACAGCTAATGCCGTTGAATTCCACAAAGAtt                                                                                                                                                                                                >  2:236465/1‑251 (MQ=9)                                                     aGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTaa                                                                                                                                                                                             <  1:202425/251‑1 (MQ=255)                                                       cGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATc                                                                                                                                                                                           >  1:123007/1‑251 (MQ=255)                                                                aTCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGa                                                                                                                                                                                   >  1:203432/1‑250 (MQ=255)                                                                            cccAAAACAGAAGTGTTAACGGAAGCGAAGAACCTTGCAGAGGCGTCCGAGCATTATAACGAATGGCTTCCGCATAGTTCTCTGGGTTATCGCTCGCCACGGGAATATCTGCGGTAGCGGGCTTGTAATGTGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATGTAACAAAGATTTAATCGAAGCCGACGATCGTGaaaa                                                                                                                                                                       <  2:256551/250‑1 (MQ=18)                                                                             ccAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACg                                                                                                                                                                     <  1:233461/251‑1 (MQ=255)                                                                               aaaCCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTCGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGtt                                                                                                                                                                   >  1:301898/1‑251 (MQ=255)                                                                                    ataCGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCt                                                                                                                                                              <  1:296822/249‑1 (MQ=255)                                                                                    aGACTGGTTCACGGCAGCAAAGAACCGTGCAGATGCGATCGAGCATTATAACGAATGGCATCCGCATAGTCCGCTTGGTTATCGCTCGCCACGTGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGAGAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTGATGCCTTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTaaa                                                                                                                                                                <  2:316498/249‑2 (MQ=18)                                                                                         ggTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTccc                                                                                                                                                         <  1:309906/251‑1 (MQ=255)                                                                                           ttAACGGCAGCAAAGAACCTTGCAGCGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCa                                                                                                                                                        >  1:258696/1‑250 (MQ=255)                                                                                                ggCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGa                                                                                                                                                  <  1:181419/251‑1 (MQ=255)                                                                                                  cagcaAAGAACCTTGCAGAGGCGTTCGTGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACTGGAATATCTGCGGCAGCGTGCTTGTAAGGGGTTAAGTGATAACAGATGTCTGGAAATAAAGGGGCAAATCTAGTTATCATCATTTGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGATGCCGAAGATCGTGAAATCGTTAACCTTTTCCAATTTGaaa                                                                                                                                                <  2:193617/251‑1 (MQ=21)                                                                                                      aaaGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGAGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCACCATTGGTTGT‑GCGTCAGCTTATGCAGTTGAATTCAACAAAGATTTAATCGAAGCCGACGATCGTGAAAACGTTAACCTTTCCCAATTTGACACTg                                                                                                                                             >  2:289555/1‑250 (MQ=255)                                                                                                       aaGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGAt                                                                                                                                           <  1:201623/251‑1 (MQ=255)                                                                                                                 cAGAGGCGTTAGAGCATTATAAATAGTGGCATCAGCATAGTTAGCTGGGTTATCGCTCGCAACGGGAATATCTAAGGCAGCGGGCTTGTAATGGGTTAAGAGATAACAGATGTATGGAAATATAGGGGCAAATCCAGTTATCATCTTTGGG‑GTAGCGTCAGCTTATGCCGTAGAATTCAACAACGATTTACTCGAAGCCGTAGATCGGGAAAACGTTACACTTTCACAATTTAATACTTACGGCCAAtaa                                                                                                                                  >  2:286857/1‑248 (MQ=255)                                                                                                                          tCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGc                                                                                                                         <  1:41973/250‑1 (MQ=255)                                                                                                                                     aaCGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACt                                                                                                              <  1:55506/250‑1 (MQ=255)                                                                                                                                       cGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAg                                                                                                           <  1:106982/251‑1 (MQ=255)                                                                                                                                            ggCATCCGCATAGTGCGCGTGTTTATCGCGCGCCACGGGAATATCGGAGGAAGAGTGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATAGGGAAAACGTTAACCGTTCCAAAGGGGAAAGGGATGGCAAATTACCCGTCGGCAAATATTCACTAAGCAct                                                                                                       <  1:318699/250‑1 (MQ=25)                                                                                                                                                  cGCATAGTGCGCTGGGATATCGCTCGCCACGGGAATCTCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGAtt                                                                                                 >  2:66872/1‑250 (MQ=255)                                                                                                                                                        gTGCGCTGGGTTATCGCTCGCCACTGGAATATCTGCTGCTGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTTAATTCAACAAAGATTTAATCGAAGCTGAAGATCGTGAAATCTTTAACCTTTCCCAATTTGAAACTTATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATtaataata                                                                                          <  1:347597/251‑1 (MQ=37)                                                                                                                                                            gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAagag                                                                                      <  1:59162/251‑1 (MQ=255)                                                                                                                                                                  ttATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGaaa                                                                                                                                                <  2:96652/187‑1 (MQ=255)                                                                                                                                                                  ttATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGaaa                                                                                                                                                >  1:96652/1‑187 (MQ=255)                                                                                                                                                                                                  gggattgtaggggattAAGTGAAAACAGATGTATGTAGATATAGGGGCAAATCAAGTTATCATCATTGGAAGT‑GCGACAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGATGCCGAAGAGCGTGAAAACGTTAACCACTCCCAATTTAAAAATGATGGCCAATTACCCGAAGACCAATATTCACTAAACAATCTGATTAATAATAAGAGACCGCAAATCTACCTTGACCTCCATTGGATATAAAttt                                                >  2:208131/15‑250 (MQ=255)                                                                                                                                                                                                        gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAAcc                                          <  1:105161/251‑1 (MQ=255)                                                                                                                                                                                                                   aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGATGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                               <  1:254881/251‑1 (MQ=255)                                                                                                                                                                                                                      gATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAAATCGTGAAAACGTTAACCTTTCCCACTTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCACTGGGTATTAATTGATAACCAAACTGCAGTTTGc                            >  1:381713/1‑251 (MQ=255)                                                                                                                                                                                                                                                gCAAATCCAGGTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattatt  <  1:68308/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGT‑GCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATT  >  NC_000913/3187179‑3187668

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.