Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	REL606	4,524,522	IS186 (+) +6 bp	coding (494‑499/549 nt)	fimA →	major type 1 subunit fimbrin (pilin)

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2772854 =	NA (NA)	65 (1.150)	43/144	0.2	100%	noncoding (1/1343 nt)	IS186	repeat region
	?	REL606	= 4524527	0 (0.000)					coding (499/549 nt)	fimA	major type 1 subunit fimbrin (pilin)
*	?	REL606	= 2774196	NA (NA)	50 (0.890)	41/144	0.3	100%	noncoding (1343/1343 nt)	IS186	repeat region
	?	REL606	4524522 =	0 (0.000)					coding (494/549 nt)	fimA	major type 1 subunit fimbrin (pilin)

ACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2774119‑2774196 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGA  >  REL606/4524522‑4524595                                                                                                                                                          ACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCGGATGGTTCAACCCTTAAGTTAGCGCTTATG                                                                              <  1:1703536/76‑1  CGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGG                                                                             <  1:114689/76‑1  CGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGG                                                                             >  1:1134843/1‑76  CGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGG                                                                             <  2:1423972/76‑1   GCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG                                                                            <  1:32186/76‑1   GCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG                                                                            >  2:1412496/1‑76   GCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG                                                                            <  2:466165/76‑1   GCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG                                                                            <  2:754786/76‑1   GCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG                                                                            <  2:77963/76‑1   GCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG                                                                            >  2:842615/1‑76    CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGC                                                                           >  2:412949/1‑76     GAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCC                                                                          >  1:501822/1‑76      AGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCG                                                                         <  2:1003619/76‑1      AGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCG                                                                         >  1:640984/1‑76       GGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGC                                                                        <  1:1478383/76‑1           TCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACC                                                                    <  1:1566612/76‑1            CGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCC                                                                   >  2:235654/1‑76             GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCC                                                                  >  1:770901/1‑76               CTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGG                                                                >  1:215817/1‑76               CTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGG                                                                <  1:673136/76‑1                TTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGG                                                               >  2:1862109/1‑76                 TAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGT                                                              >  1:460119/1‑76                  AACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTG                                                             >  2:1291713/1‑76                    CGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCT                                                           >  2:1045836/1‑76                    CGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCT                                                           >  1:998975/1‑76                    CGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCT                                                           >  1:894148/1‑76                         CACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAA                                                      <  1:344704/76‑1                          cCCGGGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAAT                                                     <  1:716436/75‑1                           CCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATG                                                    <  1:51435/76‑1                            CGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGC                                                   <  1:964574/76‑1                                  GTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGC                                             <  1:1679738/76‑1                                  GTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGC                                             >  1:1724078/1‑76                                  GTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGC                                             <  1:213932/76‑1                                    CGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGA                                           >  1:1517084/1‑76                                    CGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGA                                           >  2:204264/1‑76                                       ATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCT                                        >  2:92230/1‑76                                       ATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCT                                        <  2:599892/76‑1                                           TTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAA                                    <  1:953587/76‑1                                            TCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG                                   >  1:954475/1‑76                                               GATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTT                                <  1:1230911/76‑1                                                ATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTC                               <  2:1481512/76‑1                                                 TGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCA                              >  1:1762721/1‑76                                                   GTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGT                            <  1:449245/76‑1                                                   GTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGT                            <  1:1219893/76‑1                                                   GTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGT                            <  1:124436/76‑1                                                    TTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTA                           <  1:1579750/76‑1                                                     TCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAT                          >  2:858317/1‑76                                                     TCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAT                          <  2:1676888/76‑1                                                     TCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAT                          <  2:446017/76‑1                                                         CCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAT                      <  2:1356354/76‑1                                                             TAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACC                  <  1:579626/76‑1                                                              nAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCT                 >  1:544602/2‑76                                                              AAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCT                 >  2:1226858/1‑76                                                               AGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTA                <  1:528049/76‑1                                                                   AGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCA            <  2:1309319/76‑1                                                                       CTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTT        >  1:18873/1‑76                                                                         TATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCA      <  1:683092/76‑1                                                                         TATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCA      <  2:789052/76‑1                                                                          ATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAG     <  1:1305219/76‑1                                                                           TGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGG    >  2:1781061/1‑76                                                                            GGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGG   <  1:371092/76‑1                                                                            GGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGG   <  1:381821/76‑1                                                                             GGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGA  <  2:1702997/76‑1                                                                             GGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGA  <  1:1239405/76‑1                                                                             GGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGA  <  2:7597/76‑1                                                                             GGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGA  <  2:1374709/76‑1                                                                                                                                                          ACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2774119‑2774196 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGA  >  REL606/4524522‑4524595

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