Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,464,157	IS5 (–) +4 bp	intergenic (+87/+92)	yfdF → / ← mlaA	protein YfdF/intermembrane phospholipid transport system ‑ outer membrane lipoprotein MlaA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	24 (0.660)	24/486	0.5	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 2464160	0 (0.000)					intergenic (+90/+92)	yfdF/mlaA	protein YfdF/intermembrane phospholipid transport system ‑ outer membrane lipoprotein MlaA
*	?	NC_000913	= 2290113	NA (NA)	35 (0.960)	32/488	0.2	97.2%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	2464157 =	1 (0.030)					intergenic (+87/+95)	yfdF/mlaA	protein YfdF/intermembrane phospholipid transport system ‑ outer membrane lipoprotein MlaA

AACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCA  >  NC_000913/2288684‑2289159                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 aaCCGTCATTTAATACTGAGTATTCTGCGGCGCCACTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCGATGGGTGAATATAAGGAAGGTGCGAATAAg                                                                                                                                                                                                                                   <  1:350357/251‑1 (MQ=255)               aCTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCAAAACCAGCTGGATATCTTTTT‑TCAGGACAACACC‑CT‑CAAC‑TCTCTCACCA‑TACTCAAACGTTGATGACAATGCTATGGGGGGATATAAGGAAGGTGCGTATAAGCGGGGTAAttcttt                                                                                                                                                                                                                     >  1:44042/1‑250 (MQ=255)                  gAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCATGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTtctc                                                                                                                                                                                                                   <  1:18614/250‑1 (MQ=255)                   attATTCTGCGGCGCCATTCATTTCAATGTTCTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                                                                                                                                                                                                                 <  1:171963/249‑1 (MQ=255)                             ggcgCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑TCAGGACAACACC‑AT‑CAAAATG‑CTCACCAATA‑TCAAAAGTTGATGCCCATGCTCTGGGGGAATATAAGGAAGGTGCGAATAAGCGGGGGAATTCTTCTTGGCGGAccca                                                                                                                                                                                                         >  1:213973/1‑246 (MQ=255)                                 ccATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCCACTGGATATCTTTTT‑TCCGG‑TAACACCCAT‑CCAC‑TGTATCACCA‑TACTCCAAAGTTTATTCCAATTCCATGGGTGGATATAAGGGAGGTGGGAATAAGCGGGGAAATTCTTCTCGGCTGACTAAGacatt                                                                                                                                                                                                   >  2:211840/1‑251 (MQ=255)                                                     aTACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGGGGATATAAGGAAGGTGCGAAAAAGCGGGGGAATTCTTTTCGGCTGACTCCGTCATTTCATTTCTTCATGTTTTAGc                                                                                                                                                                               >  1:167722/1‑251 (MQ=255)                                                              gTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCAGTACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                                                                                                                      <  1:58834/251‑1 (MQ=255)                                                               ttAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CTAGGACAAAACC‑CT‑CAAA‑TGTCTCACCC‑AACTCAACAGGTGATGCCAATGATATGGGTGGATATAAGGAAGGGGCGGAAAAACGGGGGAATTTTTCTCGGCTGTCTCAGTCATTTCATTTCTTTATGTTTTTGCTGAttttttt                                                                                                                                                                     >  1:106922/1‑250 (MQ=255)                                                                     tCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCAAGAAAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGGGAATATaagggaggg                                                                                                                                                                                                                                            >  1:144908/1‑173 (MQ=255)                                                                        gTATTATTACTTACCCCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCATCATGTTTGAGCCGATTTTTTCTCCCGTaaa                                                                                                                                                            <  2:448997/251‑1 (MQ=255)                                                                                acttacACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGa                                                                                                                                                    <  1:111940/251‑1 (MQ=255)                                                                                     cacCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGCTTTTCGTATTATTGGCGGTCCCAGCAACCCTCCGGATATTTTTTTTCCAGAAAAACAC‑‑ATACAACGTTCCTCAC‑‑‑TAATCAACAGTTGATTCCAATGGTCTGGGGGAATATAAGGAAGGTGCGCCAAATCGCGGTAATTCTTCTCCGGTGGATCAATCCTTTCAATTGTTCATGATTGTACCGATTTTTTTTCCCGTacattactttagttcc                                                                                                                                               >  2:94774/1‑236 (MQ=255)                                                                                               ttttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGGTATTTTTTT‑CCAGGACAACACC‑CT‑CAAC‑TGTCTCATCA‑TACTCAACAGTTGATGCCAATGCTATGGGGGAATATAAGGAAGGTGCGAATAAGCGGGGGAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTTAGCCGATTTTTTTTCCCGTAAATGCCTTGAATCAGCCTtttttg                                                                                                                                      >  1:119582/1‑248 (MQ=37)                                                                                                tttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTa                                                                                                                                       <  1:265284/248‑1 (MQ=255)                                                                                                                             tttCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑ACCGGACAACACC‑AT‑CCAC‑TGTCTCACCA‑TACTCAACAGGTGGTGCCAATGCTATGGGTGAATATAAGGAAGGAGCGAATAAGCGGGGACAATCTTCTCGGCTGGCTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTGAATGCCTTGAATCCACCTATTTTGACCGTTTCTTCGCCATTTAAGGCGTTAtagc                                                                                                       >  1:179316/1‑248 (MQ=21)                                                                                                                                                ccATAAACCAGCTGGATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGTATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGAtctc                                                                                    <  1:277810/251‑1 (MQ=255)                                                                                                                                                             ggATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACGCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTGTTTAGTGAGATCTCTCCCACTGACg                                                                         <  1:94774/249‑1 (MQ=255)                                                                                                                                                                          ccAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGCCGTCTCATTTGGTCCg                                                           >  2:74485/1‑251 (MQ=255)                                                                                                                                                                                        aT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGGGAATATAAGGAAGGTGCGAATAAGCGGGGACATTCTTCTCGACTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGtt                                              >  1:27321/1‑251 (MQ=255)                                                                                                                                                                                                     cACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGAtctct                                                                                   >  1:348680/1‑203 (MQ=255)                                                                                                                                                                                                     cACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGAtctct                                                                                   <  2:348680/203‑1 (MQ=255)                                                                                                                                                                                                       ccA‑TACTCAACAGTTGTTGCCACTGCTATGTGTGAATATAAGGAAGGTGCGAATAAGCGGGGAACTTCTTCTCGGCTGACTCAGGCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTGTTTCGACCTTTTCTTCGCCATTTAAGTCGTTATCCCCCGTTTTTGGGGAGATATATCCCACTGAAGTATCATTTGGTCCGCCCAAAACAGGTTGGCCAGCgctaa                                  >  2:221329/1‑246 (MQ=14)                                                                                                                                                                                                                accAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGc                         <  2:96741/249‑1 (MQ=255)                                                                                                                                                                                                                              aTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGt           <  1:170557/251‑1 (MQ=37)                                                                                                                                                                                                                                        tGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCGCGGCTGACTCAGACATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGACTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCAGCAACAGGTTGGCCAGCGTGCATAACATCGCCAGTTGGTTATTGtttttacgaa  >  2:251181/1‑245 (MQ=11)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 AACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTT‑CCAGGACAACACC‑AT‑CAAC‑TGTCTCACCA‑TACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCA  >  NC_000913/2288684‑2289159

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