Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,560,930	IS5 (+) +4 bp	coding (520‑523/831 nt)	glpG ←	rhomboid intramembrane serine protease

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3560930 =	0 (0.000)	47 (1.270)	39/254	0.1	100%	coding (523/831 nt)	glpG	rhomboid intramembrane serine protease
	?	NC_000913	3652036 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region
*	?	NC_000913	= 3560933	0 (0.000)	38 (1.020)	35/254	0.2	100%	coding (520/831 nt)	glpG	rhomboid intramembrane serine protease
	?	NC_000913	= 3653230	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

CAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCC  >  NC_000913/3651899‑3652174                                                                                                                                          |                                                                                                                                           cAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGaa                                                                                                                                 <  8:152514/149‑1 (MQ=255)    gATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAAt                                                                                                                                >  1:180757/1‑147 (MQ=255)       gAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCg                                                                                                                           >  5:30811/1‑149 (MQ=255)          aGTCTTTGCCGTAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTTCTTAATATAAATTTCACTGACGCAAATTTAGGGAATGTGCGAATAAGCgggg                                                                                                                        <  6:57455/149‑1 (MQ=255)                   cGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                               <  6:30811/149‑1 (MQ=255)                    ggAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                               >  5:135895/1‑148 (MQ=255)                    ggAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATtcttct                                                                                                              <  2:180757/149‑1 (MQ=255)                               tATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCg                                                                                                                           >  6:104567/1‑125 (MQ=255)                               tATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCg                                                                                                                           <  5:104567/125‑1 (MQ=255)                                             tGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttct                                                                                     >  5:140519/1‑149 (MQ=255)                                               acaggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcttc                                                                                   <  6:140519/149‑1 (MQ=255)                                                        gttcTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATaa                                                                                                                              <  3:39016/94‑1 (MQ=255)                                                        gttcTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATaa                                                                                                                              >  4:39016/4‑97 (MQ=255)                                                        aaCCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                                                                          <  5:83965/149‑1 (MQ=255)                                                                              ccTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                             <  5:135291/140‑1 (MQ=255)                                                                              ccTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                             >  6:135291/1‑140 (MQ=255)                                                                               cTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                   >  8:120535/1‑149 (MQ=255)                                                                                   tCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTGGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt                                                <  6:135895/148‑1 (MQ=255)                                                                                      aGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                             >  4:28807/1‑132 (MQ=255)                                                                                      aGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                             <  3:28807/132‑1 (MQ=255)                                                                                        gTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtc                                                               <  5:79324/128‑1 (MQ=255)                                                                                        gTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtc                                                               >  6:79324/1‑128 (MQ=255)                                                                                        gTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGACTTTTtc                                                               <  8:130344/128‑1 (MQ=255)                                                                                        gTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGACTTTTtc                                                               >  7:130344/1‑128 (MQ=255)                                                                                              tttCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAt                                    >  4:122969/1‑149 (MQ=255)                                                                                                  gAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                         <  2:9839/108‑1 (MQ=255)                                                                                                  gAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                         >  1:9839/1‑108 (MQ=255)                                                                                                      aTCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCg                            >  2:15543/1‑149 (MQ=255)                                                                                                          aCTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttc                        >  2:55196/1‑149 (MQ=255)                                                                                                           cTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttct                       >  2:6168/1‑149 (MQ=255)                                                                                                              aaTATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                     <  3:122969/148‑1 (MQ=255)                                                                                                              aaTATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                     <  2:19220/148‑1 (MQ=255)                                                                                                                tataAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGc                   >  2:87040/1‑148 (MQ=255)                                                                                                                    aaTTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt              <  7:120535/149‑1 (MQ=255)                                                                                                                      tttCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTaa            <  1:55196/149‑1 (MQ=40)                                                                                                                          aCTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcat                                                                                              <  7:143648/63‑1 (MQ=255)                                                                                                                          aCTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcat                                                                                              >  8:143648/1‑63 (MQ=255)                                                                                                                          aCTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg        >  7:105159/1‑149 (MQ=38)                                                                                                                                gCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGa                                                                           <  4:140226/76‑1 (MQ=35)                                                                                                                                gCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGa                                                                           >  3:140226/1‑76 (MQ=35)                                                                                                                                gCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcc  >  3:11778/1‑149 (MQ=32)                                                                                                                                          |                                                                                                                                           CAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCC  >  NC_000913/3651899‑3652174

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