Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,468,298	IS3 (+) +3 bp	coding (85‑87/921 nt)	yfdH →	bactoprenol glucosyl transferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2170171 =	NA (NA)	20 (0.700)	19/484	0.4	95.2%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	= 2468300	1 (0.040)					coding (87/921 nt)	yfdH	bactoprenol glucosyl transferase
*	?	NC_000913	= 2171428	NA (NA)	9 (0.320)	9/482	1.1	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	2468298 =	0 (0.000)					coding (85/921 nt)	yfdH	bactoprenol glucosyl transferase

GACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGAAGTCATATGAAGTGGAAATCGTTTTCATAAATGACGGCAGCAAAGACGCTACGGAGTCAATCATTAATGCTCTGGCTGTTTCAGATCCTCTAGTTGTTCCGCTGTCATTTACACGCAACTTTGGTAAAGAACCAGCATTGTTTGCAGGGTTAGACCATGCAACCGGGGATGCGATAATCCCAATTGATGTTGACCTGCAAGACCCGATTGAGGTTATTCCTCATCTTATTGAAAAATGGCAAGCAG  >  NC_000913/2468079‑2468547                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          aatgcctgcgtggaaagcttctttcattcgctgaaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggagcaTTGAAGTCATATGAAGTGGAAATCGTTTTCAt                                                                                                                                                                                                                            <  2:57232‑M2/32‑1 (MQ=255)                              ctgATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggt                                                                                                                                                                                               <  1:11172‑M2/251‑59 (MQ=255)                               ctgaaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaaccagaacctcgcttaggcctgtgtccctattacgtgggtaggatcaTTGAAGTCATATGAAGTGGAAATCGtt                                                                                                                                                                                                                                 >  1:205737‑M2/190‑216 (MQ=255)                               ctgaaagtggaaggtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcgtcggcacagttggttgggcggcctcagtccggaacaatttgaaaaccagaacgtcgcttaggcctgtgtccatattacgttggtagtagcaGTGAAGTCATATGAACTGGAAATCGtt                                                                                                                                                                                                                                 <  2:205737‑M2/27‑1 (MQ=255)                                                gcccgATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTAAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgaggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaac                                                                                                                                                                             <  2:25577‑M2/250‑77 (MQ=255)                                                ccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaTTGAAGTCATATGAAGTGGAAATCGTTTTCATAAATGACGGCAGCAAAGACGCTACGGAGTCAATCATTAATGCTCTgg                                                                                                                                                                             <  1:42121‑M2/79‑1 (MQ=255)                                                     gagaacattttatcagcagggaaataatgcgggcaacggtgtttaattatatcgaaggtgattacaatcggtggcggtggcacagttggtgtggcggcttcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattccgtgggtaggatcaTTGAAGTCATATGAAGTGGAAATCGTTTTCATAAATGACGGCAGCAAAGACGCTACGGAGTCAATCATTAATGCTCTGGCTGtt                                                                                                                                                                        <  1:37384‑M2/84‑1 (MQ=255)                                                                            cGCCATCAGCCTGGTGTGCGGGTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGGCTTCAATGAAGAAGAAGCGATTCCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatctgacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaataaaccccgtaaa                                                                                                                                                 >  2:267297‑M2/1‑147 (MQ=255)                                                                               cATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaa                                                                                                                                                          <  1:167036‑M2/239‑96 (MQ=255)                                                                               cATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaa                                                                                                                                                          >  2:167036‑M2/1‑144 (MQ=255)                                                                               atgcgggcaacggtgtttaattatatcgtatgtgattacaatcggtggcggcggcacagttggtgtggaggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaTTGAAGTCATATGAAGTGGAAATCGTTTTCATAAATGACGGCAGCAAAGACGCTACGGAGTCAATCATTAATGCTCTGGCTGTTTCAAATCCTCTAGTTGTTCCGCTGTc                                                                                                                                              >  1:253377‑M2/142‑251 (MQ=255)                                                                                                  ttCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgc                                                                                                                            <  1:298937‑M2/250‑126 (MQ=255)                                                                                                     gTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgcagtg                                                                                                                        <  1:76122‑M2/251‑130 (MQ=255)                                                                                                               tgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccgggacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtggggaggatcaTTGAAGTCATATGAAGTGGAAATCGTTTTCATAAATGACGGCAGCAAAGACGCTACGGAGTCCATACTTAATGCTCTGGCTGTTTCAGATCCTCTAGTTGTTCCGCTGACATTTAAACGCAACTTTGGTAAAGAACCAGCAt                                                                                                              >  1:177975‑M2/110‑251 (MQ=255)                                                                                                               aGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagt                                                                                                                                                                >  1:182530‑M2/1‑112 (MQ=255)                                                                                                               aGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagt                                                                                                                                                                <  2:182530‑M2/201‑90 (MQ=255)                                                                                                                 attacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaTTTAAGTCATATGAAGTGGAAATCGTTTTCATAAATGACGGCAGCAAAGACGCTACGGAGTCAATCATTAATGCTCTGGCTGTTTCAGATCCTCTAGTTGTTCCGCTGTCATTTACACGCAACTTTGGTAAAGAAccagccttt                                                                                                            >  2:86437‑M2/108‑247 (MQ=255)                                                                                                                                           aaGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgcagtgaagccctgaagcttgctgaacgcatcggtgttactgcc                                                                                  <  1:30357‑M2/251‑168 (MQ=255)                                                                                                                                                     ttGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgcagtgaagccctgaagcttgctgaacgcatcggtgttactgccgcagcccgtg                                                                        <  2:82315‑M2/251‑178 (MQ=255)                                                                                                                                                       gTAGTTACTGGCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgcagtgaagccctgaagcttgctgaacgcatcggtgttactgccgcagcccgtgaa                                                                      >  2:247908‑M2/1‑72 (MQ=255)                                                                                                                                                              cTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgcagtgaagccctgaagcttgctgaacgcatcggtgttactgccgcagcccgtgaactcagcc                                                               <  1:158800‑M2/251‑187 (MQ=255)                                                                                                                                                                gTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgcagtgaagccctgaagcttgctgaacgcatcggtgttactgccgcagcccgtgaactcagcct                                                              >  2:22599‑M2/1‑63 (MQ=255)                                                                                                                                                                     cAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgccgtgaagccctgaagcttgctgaacgcagcgttgtttctgccgcagcccgtgaactcagcctgtatg                                                         >  1:190942‑M2/1‑58 (MQ=255)                                                                                                                                                                      aaTGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcacctgaatttcgcagtgaagccctgaagcttgatgaacgcatcggtgttactgccgcagcccgtgaactcagcctgtctga                                                        <  2:5437‑M2/250‑194 (MQ=255)                                                                                                                                                                                    gATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccggaaacagcattcgcctgaatttcgcagtgaagccctgaagcttgctgaacgcatcggtggtacggccgcagcccgtgaactcagcctgtatgaatcacaactctacaa                                         >  1:18029‑M2/1‑43 (MQ=255)                                                                                                                                                                                    gATACAAATTTTTTATAAAACGGTACGTGAATTCGAAGCATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgcagtgaagccctgaagcttgctgaacgcatcggtgttactgccgcagcccgtgcactcagcctgtatgaatcacaactctacac                                         <  2:291960‑M2/251‑209 (MQ=255)                                                                                                                                                                                                            tacgggggtaggatccTTGAAGTCATATGATGTGGAAATCGTTTTCATAAATGACGGCAGCAAAGACGCTACGGAGTCAATCATTAATGCTCTGGCTGTTTCAGATCCTCTAGTTGTTCCGCTGTCATGTACACGCAACTTTGGTAAAGAACCAGCATTGTTTGCAGGGTTAGACCATGCAACCGGGGATGCGATAATCCCAATTGATGTTGACCTGCAAGTCCCGATTGAGGTTATTCCTCATCTTAtt                  <  2:170779‑M2/234‑1 (MQ=255)                                                                                                                                                                                                             cGTGAATTCGAAGAATTGtgatcttacccagcaatagtggacacgcggctaagtgagtaaactctcagtcagaggtgactcacatgacaaaaacagtatcaaccagtaaaaaaccccgtaaacagcattcgcctgaatttcgcagtgaagccctgaagcttgctgaacgcatcggtgttactgccgcagcccgtgaactcagcctgtatgaatcacaactctacaactggcgcagtaaacagcaaaatcag                >  1:227790‑M2/1‑18 (MQ=255)                                                                                                                                                                                                                           aTTGAAGTCATATCACGTGGAAATCGTTTTCATAAATGACGGCAGCAAAGACGCTACGGAGTCAATCATTAATGCTCTGGCTGTTTCAGATCCTCTCGTTGTTCAGCTGTCATTTACCCGCGACTTTGGTAAAGAACCAGCATTGTTTGCAGGGTTCGACCATGCAACCGGGGATGCGATAATCCCAATTGATGTTGACCTGCACGACCCGCTTGAGGTTATTCCTCATCTTTTTGATAATTGGCaggctg  >  2:263147/1‑246 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          GACACTGAGTGCTACTGTTGGATGGGCTGCTGATAGATGCGCACTTCCCCCGATGATAACTCTTGTCACCTTCTCCGCCATCAGCCTGGTGTGCGGTTTCGTCTATTCAAAGTTCATTGTCTTTAGGGATGCGAAATGAAGATATCTCTTGTAGTTCCTGTCTTCAATGAAGAAGAAGCGATACCAATTTTTTATAAAACGGTACGTGAATTCGAAGAATTGAAGTCATATGAAGTGGAAATCGTTTTCATAAATGACGGCAGCAAAGACGCTACGGAGTCAATCATTAATGCTCTGGCTGTTTCAGATCCTCTAGTTGTTCCGCTGTCATTTACACGCAACTTTGGTAAAGAACCAGCATTGTTTGCAGGGTTAGACCATGCAACCGGGGATGCGATAATCCCAATTGATGTTGACCTGCAAGACCCGATTGAGGTTATTCCTCATCTTATTGAAAAATGGCAAGCAG  >  NC_000913/2468079‑2468547

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