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)	77 (1.350)	54/144	0.0	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)	69 (1.210)	48/144	0.1	100%	noncoding (1343/1343 nt)	IS186	repeat region
	?	REL606	4524522 =	0 (0.000)					coding (494/549 nt)	fimA	major type 1 subunit fimbrin (pilin)

GTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2772931‑2772854 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAAC  <  REL606/4524527‑4524453                                                                                                                                                           GTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATG                                                                               >  2:1009944/1‑76 GTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATG                                                                               <  2:286904/76‑1  TTATCGTGAGAGTAATTCATCGGCACGTTAAATCATCTCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGG                                                                              >  2:1246941/1‑76  TTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAGCCCTTAAGTTAGCGCTTATGG                                                                              <  1:1337587/76‑1  TTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGG                                                                              >  2:1706194/1‑76  TTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGG                                                                              >  2:1817495/1‑76  TTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGG                                                                              <  1:352007/76‑1   TATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG                                                                             >  2:832210/1‑76   TATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG                                                                             >  2:129224/1‑76   TATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG                                                                             <  1:706681/76‑1   TATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG                                                                             >  1:235953/1‑76   TATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG                                                                             >  1:1701304/1‑76   TATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG                                                                             >  1:1484044/1‑76    ATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGT                                                                            >  2:1578968/1‑76    ATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGT                                                                            >  2:1773928/1‑76    ATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGT                                                                            >  1:121154/1‑76    ATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGCAATACCACAACCCTTAAGTTAGCGCTTATGGGT                                                                            >  1:321175/1‑76      CGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTG                                                                          <  1:861912/76‑1       GTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGTGTAATACCACAACCCANAAGTTAGCGCTTATGGGTTGC                                                                         >  2:301569/1‑76       GTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGC                                                                         <  1:686962/76‑1       GTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGC                                                                         >  2:389291/1‑76       GTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGC                                                                         >  2:133840/1‑76        TGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCG                                                                        >  2:1477506/1‑76        TGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCG                                                                        >  1:574138/1‑76        TGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCG                                                                        <  2:1583787/76‑1          AGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGC                                                                      >  1:313689/1‑76           GAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCC                                                                     >  1:497945/1‑76               AATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGG                                                                 >  1:637240/1‑76               AATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGG                                                                 >  2:1358565/1‑76               AATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGG                                                                       <  1:555522/70‑1               AATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGG                                                                       >  2:555522/1‑70                ATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGT                                                                <  1:1385792/76‑1                 TTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTT                                                               <  2:1329915/76‑1                  TCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTG                                                              <  2:445828/76‑1                  TCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTG                                                              <  2:445818/76‑1                   CATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGC                                                             >  2:132534/1‑76                       GGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAA                                                         <  2:425626/76‑1                         CACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATA                                                       <  1:1459043/76‑1                         CACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATA                                                       >  1:810215/1‑76                         CACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATA                                                       <  2:173901/76‑1                           CGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAAC                                                     >  1:1571287/1‑76                              TAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCG                                                  <  2:780979/76‑1                               AAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGC                                                 >  2:1414459/1‑76                               AAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGC                                                 >  1:1244955/1‑76                                   CATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGG                                             <  1:816924/76‑1                                    ATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGA                                            >  2:795163/1‑76                                      ATCAGGCGTAATACCACACCCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAAC                                          <  1:888381/76‑1                                       TCAGGCGTAATCCCACACCCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACG                                         <  1:915561/76‑1                                       TCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAAag                                         >  1:93774/1‑74                                       TCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACG                                         <  2:435551/76‑1                                        CAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCTAAATAACGCGCCTGGAACGG                                        >  2:707401/1‑76                                        CAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGG                                        >  2:568594/1‑76                                        CAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGG                                        >  1:231332/1‑76                                         AGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGA                                       >  1:935897/1‑76                                         AGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGA                                       <  1:577724/76‑1                                            CGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATG                                    >  1:352025/1‑76                                             GTAATACCACAACCCTTAAGTTNGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGG                                   >  2:255592/1‑76                                             GTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGG                                   >  1:1284343/1‑76                                               AATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTA                                 >  1:437952/1‑76                                               AATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTA                                 >  2:1380228/1‑76                                               AATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAA                                      <  2:442967/71‑1                                               AATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAA                                      >  1:442967/1‑71                                                ATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGGATGGTAT                                >  1:585442/1‑76                                                ATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTAT                                <  1:1418273/76‑1                                                  ACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTG                              >  1:1076391/1‑76                                                  ACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTG                              <  1:1750694/76‑1                                                  ACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTG                              <  1:894893/76‑1                                                  ACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTG                              <  1:717007/76‑1                                                       AACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCC                         <  1:524117/76‑1                                                        ACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCG                        >  1:468734/1‑76                                                        ACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCG                        >  1:1179372/1‑76                                                         CCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGT                       >  1:535245/1‑76                                                          CCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTT                      <  2:1634107/76‑1                                                           CTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTA                     <  1:1389492/76‑1                                                            TTAAGTTAGCGCTGATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTAT                    <  1:648970/76‑1                                                             TAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATT                   <  2:1026019/76‑1                                                              AAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTC                  <  2:1514897/76‑1                                                              AAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTC                  <  2:1173224/76‑1                                                                GTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAG                >  2:679113/1‑76                                                                  TAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGG              >  2:168684/1‑76                                                                    GCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTT            <  1:937907/76‑1                                                                      GCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGT          >  2:823059/1‑76                                                                        TTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTT        <  2:909621/76‑1                                                                        TTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTT        <  2:1469589/76‑1                                                                        TTATGGGTTACGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTT        >  1:1383457/1‑76                                                                         TATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTC       <  2:649353/76‑1                                                                         TATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTC       >  2:281146/1‑76                                                                          ATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCT      <  1:1227636/76‑1                                                                           TGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTG     >  1:1787182/1‑76                                                                           TGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTG     >  2:1099280/1‑76                                                                            GGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGA    <  1:1105454/76‑1                                                                            GGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGA    >  1:1021965/1‑76                                                                             GGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAA   <  1:985021/76‑1                                                                              GTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAAC  <  2:798003/76‑1                                                                              GTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAAC  <  2:81700/76‑1                                                                                                                                                           GTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2772931‑2772854 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTCTGAAC  <  REL606/4524527‑4524453

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