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	annotation	gene	product
*	?	REL606	15386 =	NA (NA)	54 (1.080)	35/64	0.1	noncoding (1/1343 nt)	IS186	repeat region
	?	REL606	= 4524527	57 (1.040)				coding (499/549 nt)	fimA	major type 1 subunit fimbrin (pilin)
*	?	REL606	= 16728	NA (NA)	6 (0.120)	31/64	0.2	noncoding (1343/1343 nt)	IS186	repeat region
	?	REL606	4524522 =	57 (1.040)				coding (494/549 nt)	fimA	major type 1 subunit fimbrin (pilin)

GATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTA  >  REL606/15351‑15418                                    |                                 gATATGCGAGCACCTGAAAAGTGGCGGGGATCAATc                                  >  1:1875236/1‑36 (MQ=255)  atatGCGAGCACCTGAAAAGTGGCGGGGATCACTcc                                 <  1:5751475/36‑1 (MQ=255)  atatGCGAGCACCTGAAAAGTGGCGGGGATCACTcc                                 <  1:2349850/36‑1 (MQ=255)   tatGCGAGCACCTGAAAAGTGGCGGGGATCACTccc                                <  1:1888766/36‑1 (MQ=255)   tatGCGAGCACCTGAAAAGTGGCAGAGATCACTccc                                <  1:5531498/36‑1 (MQ=25)    atGCGAGCACCTGAAAAGTGGCGGGGGTCACACCCa                               >  1:806109/1‑36 (MQ=25)     tGCGAGCACCTGAAAAGTGGCGGGGGTCTCTCCCAt                              <  1:3931920/36‑1 (MQ=25)      gCGAGCACCTGAAAAGTGGCGGGGGTCACTTCCATa                             >  1:4942254/1‑36 (MQ=25)      gCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCCTa                             >  1:1065210/1‑36 (MQ=25)       cGAGCACCTGAAAAGTGGCGGGGGTCACTCCCCTaa                            >  1:4944976/1‑36 (MQ=25)         aGCACCTGAAAAGTGGCGGGGGTCACTCGCATTAgc                          >  1:5815570/1‑36 (MQ=21)           cACCTGAAAAGTGGCGGGGATCACTCCCATAAgcgc                        <  1:4817015/36‑1 (MQ=255)            agcTGAAAAGTGGCGGGGATCACTCCCATAAGCGCt                       <  1:2039304/34‑1 (MQ=255)             ccTGAAAAGTGGCGGGGGAGAATCCCATAAGCGGTa                      >  1:7244883/1‑36 (MQ=12)              cggAAAAGTGGCGGGGATCACTCCCATAAGCGCTaa                     <  1:4379684/34‑1 (MQ=255)              cTGAAAAGTGGCGGGGATCACCCCCCCAAGCGCttc                     >  1:3456208/1‑34 (MQ=18)                gAAAAGTGGCGGGGATCACTCCCATAAGCGCTaaat                   >  1:430605/1‑34 (MQ=255)                 aaaaGTGGCGGGGATCACTCCCATAAGCGCTAACtt                  <  1:3951085/36‑1 (MQ=255)                   aaGTGGCGGGGATCACTCCCATAAGCGCTGACTTaa                >  1:6621895/1‑36 (MQ=255)                   aaGTGGCGGGGAACACTCCCATAAGCGCTAACTTaa                >  1:5085748/1‑36 (MQ=255)                     gTGGCGGGGATCACTCCCATAAGCGCTAACTTAAgg              <  1:4433815/36‑1 (MQ=255)                     atGGCGGGGATCACTCCCATAAGCGCTAACTTAAgg              <  1:2915679/35‑1 (MQ=255)                       ggCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt            <  1:6087375/36‑1 (MQ=40)                       ggCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt            >  1:6437636/1‑36 (MQ=40)                        gCGGGGATCACTCCCATAAGCGCTAACTTATGGGtt           >  1:6430896/1‑36 (MQ=255)                        gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt           >  1:4304780/1‑36 (MQ=38)                        gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt           >  1:422646/1‑36 (MQ=38)                         cGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg          >  1:2162790/1‑36 (MQ=38)                           gggATCACTCCCATAAGCGCTAACTTAAGGGTtgtg        >  1:172928/1‑36 (MQ=255)                           gggATCACTCCCATAAGCGCTAACTTAAGGGTtgtg        <  1:3953775/36‑1 (MQ=255)                            ggATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg       <  1:5702961/36‑1 (MQ=255)                              aTCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTa     <  1:6498374/36‑1 (MQ=255)                              aTCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTa     <  1:6752435/36‑1 (MQ=255)                               tCACTCCCATAAGCGCTACCTTAAGGGTTGTGGTAt    >  1:5454557/1‑36 (MQ=255)                                cACTCCCATAAGCGCTAACTTACGGGTTGTGGTAtt   >  1:1630059/1‑36 (MQ=255)                                cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt   >  1:7125636/1‑36 (MQ=38)                                cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt   >  1:624417/1‑36 (MQ=38)                                cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt   >  1:6198710/1‑36 (MQ=38)                                 aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa  >  1:4893798/1‑36 (MQ=35)                                    |                                 GATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTA  >  REL606/15351‑15418 Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG