Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	CP009273	3,178,254	IS5 (+) +4 bp	28.4%	coding (56‑59/291 nt)	yqiC →	BMFP family putative fusogenic protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	2282398 =	NA (NA)	14 (0.200)	11/274	NT	25.8%	noncoding (1195/1195 nt)	IS5	repeat region
	?	CP009273	3178254 =	39 (0.580)					coding (56/291 nt)	yqiC	BMFP family putative fusogenic protein
*	?	CP009273	= 2283592	NA (NA)	17 (0.240)	14/276	NT	29.7%	noncoding (1/1195 nt)	IS5	repeat region
	?	CP009273	= 3178257	39 (0.570)					coding (59/291 nt)	yqiC	BMFP family putative fusogenic protein

ATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCA  >  CP009273/2282253‑2282539                                                                                                                                                  |                                                                                                                                              gttCTTTTATATCCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATaaggaa                                                                                                                                            <  1:392389/148‑1 (MQ=255)     ttttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGAGATCTTTTTCCAGGACAACACCATCAACTGTCTCCCCCTGATCAAAAGCTGGGGGCACGGCTATGGGCGACTATAAGGAAGGTg                                                                                                                                        >  1:1152957/1‑149 (MQ=255)      tttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTCTGGGTGAATATAAGGAAGGtct                                                                                                                                       >  1:522072/1‑147 (MQ=255)      tttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGc                                                                                                                                       <  1:1238238/149‑1 (MQ=255)      tttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGc                                                                                                                                       >  1:571480/1‑149 (MQ=255)      tttATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGc                                                                                                                                       <  1:798660/149‑1 (MQ=255)          tatGCTGTTCTACGCCAGGATTGATTATCGGATTATTGTAGGTCCCATAAACCAACTGGATATCGTTTTCTAGGAAAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGATATGGGTGAATATAAGTAAGGGGCGAAt                                                                                                                                   >  2:376914/1‑149 (MQ=255)                   ctcCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAAAGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGa                                                                                                                          >  2:20819/1‑149 (MQ=255)                       gCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAtt                                                                                                                      <  2:1101065/149‑1 (MQ=255)                        ccAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttc                                                                                                                     <  1:556569/149‑1 (MQ=255)                          aGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                                   <  1:988380/149‑1 (MQ=255)                            gattgattTTCGTATTATTGGCGGTCACATAGATCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATtcttct                                                                                                                 >  2:578132/1‑149 (MQ=255)                                        tattatGGGCGGTCCCATCAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGCGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAg                                                                                                     <  1:279468/149‑1 (MQ=255)                                                        aTAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTCTGGGTGAATATAAGGGAGGTGCGAATAAGCGGGGGCATTCTTCTCGGCTGACTCCGTCCTTTCATttcttct                                                                                     >  2:651175/1‑148 (MQ=255)                                                               aGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTCCATGTTTGa                                                                              <  1:328723/149‑1 (MQ=255)                                                                     atatCTTTTTCCAGGAAAACACCATCAACTGTTCCACCATACTCAACAGTTGATGCCAATGCTTTGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAt                                                                        <  2:700747/149‑1 (MQ=255)                                                                              tCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTcc                                                               <  1:313687/149‑1 (MQ=255)                                                                                 aGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGGGCCGATTTTTTTTCCCGt                                                            >  2:24723/1‑149 (MQ=255)                                                                                 aGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATAGGAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGGGCCGATTGTTTCTCCCGt                                                            >  2:316564/1‑149 (MQ=255)                                                                                       acacCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                      >  2:956495/1‑149 (MQ=255)                                                                                          ccATCCGCTGTCTCACCATACGCAACAGTTGACGCCAATGCTATTGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt                                                   <  2:561523/149‑1 (MQ=18)                                                                                            aTCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGa                                                 <  1:1085580/149‑1 (MQ=255)                                                                                             tCAACTGTCTCACCACACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                                <  1:561598/149‑1 (MQ=255)                                                                                              taaCTGTGTCACCATACGTCACTGTTGTTGCGAATGCTACGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTTTTCTGGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAt                                               <  2:2443/148‑1 (MQ=17)                                                                                               aaCTGTCTCACCATACTCCCCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATCTTTTCTCCCGTAAATGCCTTGAAta                                              >  1:910668/1‑148 (MQ=255)                                                                                                aCTGTCTCACCATACTCAACAGTTGCTGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCCATTTTTTCTCCCGTAAATGCCTTGAATCa                                             <  1:759572/149‑1 (MQ=255)                                                                                                 gtGTCCCACCAGACTCAACGGGTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                            <  2:722326/148‑1 (MQ=255)                                                                                                 ctggcccaacATACGCAACAGTTGATGCCAATGCTATGGGTGAATACAAGGAAGGTGCGAATAAGCGGGCAAATTCTTCTCGGCTGACTCAGTCATTTGATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                            <  2:317859/140‑1 (MQ=17)                                                                                                 cTGTCTCACCCTACTCAACCGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                            <  2:415079/149‑1 (MQ=255)                                                                                                  tgcttcACCATACTCAACTGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGCGCTAATAAGCTGTGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                           <  2:369262/145‑1 (MQ=18)                                                                                                  tGTCTCACCATCCTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                           <  1:578132/149‑1 (MQ=255)                                                                                                        aCCATACTCAACAGTTGATGCCAATGCTATCGTTGTATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTGTTTCTCCCGTAAATGCCTTGAATCAGCCTAttt                                     >  2:575694/1‑149 (MQ=255)                                                                                                               tCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTGATGTTTGAGCCGATTTTTTCTCGTGTAAAGGCCTTGAATCAGCCTATGTCGCCCGt                              >  1:708783/1‑149 (MQ=17)                                                                                                                   cAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttct                          >  2:285668/1‑149 (MQ=255)                                                                                                                   cAGTTGATGCCAATGCTATGGGTGAAGATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCTGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttct                          >  2:293424/1‑149 (MQ=255)                                                                                                                     gTTGCTGCCAATGCTATGGGTGAATATAAGGCAGGTGCGAATCAGCTGGTAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGTCTTGAATCCGCCTATTTAGACCGTttcttc                        <  1:648351/149‑1 (MQ=17)                                                                                                                      ttGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg                       >  1:849311/1‑149 (MQ=255)                                                                                                                         aTGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTCTTTAGAGCGTTTCTTCGCCa                    >  1:430626/1‑149 (MQ=255)                                                                                                                                tGCTATGGGTGACTATAAGGAAGGTGCGAATAAGCGGGGAAGTTCTTCTCGGCTGACTCAGTCGTTTACTTTCGTCATGTTTGAGCCGATTTTTTCTCCCGTACATGCCTTGAATCAGCCGCTGTTGACCGTTGCTTCGCCAGTTagac             >  2:1152903/1‑146 (MQ=14)                                                                                                                                     tGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGTAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTa        <  1:1084130/149‑1 (MQ=38)                                                                                                                                        gTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGATGACTCAGGCGTGTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcc     >  1:862349/1‑149 (MQ=255)                                                                                                                                         tGAATATTAGTAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGGCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATccc    <  1:376914/149‑1 (MQ=16)                                                                                                                                           aaTATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCTGTTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGGTATCCCCa  <  1:48085/149‑1 (MQ=12)                                                                                                                                           aaTATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCa  >  1:509111/1‑149 (MQ=32)                                                                                                                                                  |                                                                                                                                              ATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCA  >  CP009273/2282253‑2282539

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