Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	CP009273	3,178,254	IS5 (+) +4 bp	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)	39 (0.660)	31/290	1.3	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	CP009273	3178254 =	0 (0.000)					coding (56/291 nt)	yqiC	BMFP family putative fusogenic protein
*	?	CP009273	= 2283592	NA (NA)	52 (0.870)	43/292	0.4	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	CP009273	= 3178257	0 (0.000)					coding (59/291 nt)	yqiC	BMFP family putative fusogenic protein

CGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAG  >  CP009273/2283445‑2283735                                                                                                                                                    |                                                                                                                                                cGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCa                                                                                                                                                <  2:775594/149‑1 (MQ=11)         gtaaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATaa                                                                                                                                             <  1:755831/142‑1 (MQ=17)                tctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTa                                                                                                                                      <  1:354463/144‑1 (MQ=38)                tctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccac                                                                                                                                 <  2:123026/149‑1 (MQ=255)                        ggTCCGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGTCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACTCTTTAGCCACCAg                                                                                                                              <  1:343488/144‑1 (MQ=255)                             gACGGCGCTTACTGCTGCATTCACTGCCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATTATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAg                                                                                                                         <  1:423662/144‑1 (MQ=255)                               cGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTACGGACTTACTTTAACGCTGTAGCCACCAGAACAGATATTGc                                                                                                                  <  2:761177/149‑1 (MQ=17)                                ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCg                                                                                                                 >  2:409325/1‑149 (MQ=255)                                        ctggtgAAGTCACTGTCGGCGTAGGTAAGTTGATGACTCATGATGAACCCTGTACTAAGCCCTCAGGTGACACACATGTTCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGaa                                                                                                              <  1:808130/144‑1 (MQ=18)                                               aTTCCCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGaacaa                                                                                                       <  1:12591/144‑5 (MQ=255)                                                ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaaa                                                                                                 <  2:113358/149‑1 (MQ=255)                                                   aCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAgaga                                                                                                   >  1:520026/1‑144 (MQ=255)                                                     tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCCGAACGACAAAGAGaaa                                                                                                 >  1:346031/1‑144 (MQ=255)                                                          gCGAAGGTAAGTTGATGACTCATGATGAACCGTGTACTATGGCTCCAGATGACAAACATGATCTCATATCAGGTACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAAGCAAACAGaa                                                                                            >  1:627632/1‑144 (MQ=18)                                                            gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGACACAGAACCAGAtt                                                                                     >  2:5584/1‑149 (MQ=255)                                                               ggTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGa                                                                                       <  1:816266/144‑1 (MQ=255)                                                                gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCGCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAATGACAAAGAGAAACAGAACCAGATTGATg                                                                                 <  2:674926/149‑1 (MQ=255)                                                                         tgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCAtt                                                                             <  1:358413/144‑1 (MQ=255)                                                                               atgatgAACCCTGTTCTAAGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATc                                                                  >  2:759989/1‑149 (MQ=255)                                                                                    gaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATc                                                                  >  1:181156/1‑144 (MQ=255)                                                                                        ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaat                                                         >  2:802801/1‑149 (MQ=255)                                                                                         cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCCGATTGATGCATTGAGCTTTCATCCtcag                                                             >  1:504017/1‑141 (MQ=255)                                                                                             tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTTAGCTTTCATCCTATGAaattaatt                                                    <  2:770294/149‑1 (MQ=255)                                                                                                aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCt                                                 <  2:313475/149‑1 (MQ=255)                                                                                                aTGGCTCCAGAGCACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCACACCACTCCATTGAGCTTTCATCCTATGAAATTAATTGCt                                                 <  2:1020698/149‑1 (MQ=37)                                                                                                 tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAATCAGATTGATGCATTGAGCTTTCATCCTATGAaattaat                                                     <  1:518525/144‑1 (MQ=255)                                                                                                  ggCTCCAGATGACAAACATGATCTCATATCAGGGACTAGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCGTTCATCCTCTGAAATTAATTGCTGt                                               >  2:658369/1‑149 (MQ=255)                                                                                                   gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGAACCTTCCATCACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTCATGCATTGAGCTTTCATCCTATGAAATTAATTg                                                   <  1:245182/144‑1 (MQ=255)                                                                                                     tctagttggcAAACATGATCTCATATCAGGGACTTGTTCTCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCt                                                 <  1:897936/135‑1 (MQ=255)                                                                                                     tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGGATTGAGCTTTTATCCTGTGACACTAATTTCTGTTaa                                            >  2:799210/1‑149 (MQ=38)                                                                                                       cAGATGACAAACAGGATCACATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTaaaa                                          <  2:78961/149‑1 (MQ=255)                                                                                                          aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCCGTTaa                                            >  1:371044/1‑144 (MQ=255)                                                                                                            gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTaaaa                                          <  1:821149/144‑1 (MQ=255)                                                                                                            gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAtt                                     >  2:7107/1‑149 (MQ=255)                                                                                                               aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCa                                       >  1:111825/1‑144 (MQ=255)                                                                                                               aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCGGTAGCCTCCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATTCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCa                                       <  1:621476/144‑1 (MQ=255)                                                                                                               aaaCATGATCTCATATCAGGGACTTGTTCGCACATTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATACATGGAGCTTTCATCCTAGAAAATAAATTGCTGTTAAAAGCa                                       >  1:362872/1‑144 (MQ=39)                                                                                                                aaCATGATCTCATATCAGGGACTTGTTCGCACTTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGt                                 <  2:703396/149‑1 (MQ=255)                                                                                                                   aTGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTgg                                   <  1:941221/144‑1 (MQ=255)                                                                                                                    tGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTggg                                  <  1:666324/144‑1 (MQ=255)                                                                                                                     gATCTCATATCAAGGACTTGTTCGAACCCTCCATAACGCTGAAGCCACCAGCACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAACTGCTGTTAAAAGCATTGGGt                                 <  1:54789/144‑1 (MQ=255)                                                                                                                         tcATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAAt                        >  2:947755/1‑149 (MQ=255)                                                                                                                          cATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGa                            >  1:883292/1‑144 (MQ=255)                                                                                                                          cATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGa                            <  1:555808/144‑1 (MQ=255)                                                                                                                          cATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGa                            <  1:555688/144‑1 (MQ=255)                                                                                                                            tatCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGaaa                          >  1:180174/1‑144 (MQ=255)                                                                                                                                 gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAGTTGCTGTTAAAAGCCTTGGGTACAGAAAATACCCATAg                <  2:556267/149‑1 (MQ=255)                                                                                                                                    aCTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCAt                  <  1:781767/144‑1 (MQ=255)                                                                                                                                      ttGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATAGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGCGCTTTCACCCAAAGAAATTAATTGCTGTTAAAAGAATTGGGTACAGAAAATACCCaacg                >  1:25945/1‑141 (MQ=255)                                                                                                                                          tCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATAcc       >  2:451853/1‑149 (MQ=255)                                                                                                                                            gCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCg     <  2:566277/149‑1 (MQ=255)                                                                                                                                                    cATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGAGTACAGAAAATACCCCTAGCTCAATACCCGGAg  >  1:944937/1‑144 (MQ=255)                                                                                                                                                    |                                                                                                                                                CGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAG  >  CP009273/2283445‑2283735

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