Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	REL606	3,015,771	IS150 (–) +4 bp	coding (1192‑1195/1677 nt)	ECB_02816	KpsD protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	annotation	gene	product
*	?	REL606	588495 =	NA (NA)	49 (0.950)	33/66	0.2	noncoding (1/1977 nt)	IS150	repeat region
	?	REL606	3015771 =	57 (1.040)				coding (1192/1677 nt)	ECB_02816	KpsD protein
*	?	REL606	= 3015774	58 (1.120)	47 (0.910)	32/66	0.2	coding (1195/1677 nt)	ECB_02816	KpsD protein
	?	REL606	= 3652533	NA (NA)				noncoding (1443/1443 nt)	IS150	repeat region

CATTTCATCGTCCAACAAAATGGGTGCAGTACA...................................  <  REL606/588527‑588495 ...............................CACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  REL606/3015771‑3015805                                                                      CATTTCATCGTCCAACAAAATGGGCGCAGTACACTT                                  >  1:3707352/1‑36 (MQ=255)  ATTTCATCGTCCAACAAAATGGGTGCAGTACACcgg                                 <  1:6170468/36‑4 (MQ=1)  ATTTCATCGTCCAACAAAATGGGTGCAGTACACcgg                                 >  1:6964949/1‑33 (MQ=1)  ATTTCATCGTCCAACAAAATGGGTGCAGTACACaac                                 >  1:5425897/1‑33 (MQ=0)    TTCATCGTCCAGCAAAATGGGTGCAGTACACTTTCT                               <  1:3007280/36‑1 (MQ=255)     TCATCGTCCGACAAAATGGGTGCAGTACACTTTCTG                              >  1:1294855/1‑36 (MQ=255)     TCATCGTCCAACAAAATGGGTGCAGTACACTTTCTG                              <  1:5269797/36‑1 (MQ=255)      CATCGTCCAACAAAATGGGTGCAGTACACTTTCTGC                             >  1:791289/1‑36 (MQ=255)      CATCGTCCAACAAAATGGGTGCAGTACACTTTCTGC                             <  1:3212360/36‑1 (MQ=255)       ATCGTCCAACAAAATGGGTGCAGTACACTTTCTGCC                            >  1:7402286/1‑36 (MQ=255)        TCGTCCAACAAAATGGGTGCAGTACACTTTCTGCCC                           >  1:2109187/1‑36 (MQ=255)          GTCCAACAAAATGGGTGCCGTACACTTTCTGCCCAG                         >  1:5746435/1‑36 (MQ=255)          GTCCAACAAAATGGGTGCAGTACACTTTCTGCCCAG                         <  1:5365234/36‑1 (MQ=255)           TCCAACAAAATGGGTGCAGTACACTTTCTGCTCAGT                        <  1:6396544/36‑1 (MQ=255)           TCCAACAAAATGGGTGCAGTACACTTTCTGCCCAGT                        >  1:3464077/1‑36 (MQ=255)           TCCAACAAAATGGGTGCAGTACACTTTCTGCCCAGT                        <  1:312984/36‑1 (MQ=255)            CCAACAAAATGGGTGCAGTACACTTTCTGCCCAGTC                       <  1:4849083/36‑1 (MQ=255)             CAACAAAATGGGTGCAGTACACTTTCTGCCCAGTCA                      >  1:6048984/1‑36 (MQ=255)             CAACAAAATGGGTGCAGTACACTTTCTGCCCAGTCA                      >  1:6987641/1‑36 (MQ=255)              AACAAAATGGGTGCAGTACACTTTCTGCCCAGTCAT                     <  1:1538542/36‑1 (MQ=255)               ACAAAATTGGTGCAGTACACTTTCTGCCCAGTCATC                    <  1:3912889/36‑1 (MQ=255)                 AAAATGGGTGCAGTACACTTTCTGCCCAGTTATCCA                  >  1:3118564/1‑36 (MQ=255)                 AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  >  1:26606/1‑36 (MQ=255)                 AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  >  1:2688366/1‑36 (MQ=255)                 AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  <  1:2841118/36‑1 (MQ=255)                 AAAATGGGTGCAGTACACTTTCTGCCCAGTCATCCA                  >  1:3012721/1‑36 (MQ=255)                  AAATGGGTGCAGTACACTTTCTGCCCAGTCATCCAC                 >  1:2235403/1‑36 (MQ=255)                  AAATGGGTGCAGTACACTTTCTGCCCAGTCATCCAC                 <  1:1787012/36‑1 (MQ=255)                  AAAGGGGTGCAGTACACTTTCTGCCCAGTCATCCAC                 <  1:5811601/36‑1 (MQ=255)                   AATGGGCGCAGTACACTTTCTGCCCAGTCATCCACC                >  1:6367861/1‑36 (MQ=255)                    ATGGGTGCAGTACACTTTCTGCCCAGTCATCCACta               >  1:2764810/1‑34 (MQ=255)                     TGGGTGCAGTACCCTTTCTGCCCAGTCATCCACCAA              >  1:5523567/1‑36 (MQ=255)                     TGGGTGCAGTACACTTTCTGCCCAGTCATCCACCAA              <  1:4169868/36‑1 (MQ=255)                      GGGTGCAGTACACTTTCTGCCCAGTCATCCACCAAA             <  1:892269/36‑1 (MQ=255)                       GGTGCAGTACACTTTCTGCCCAGTCATCCACCAAcg            >  1:455081/1‑34 (MQ=255)                        GTGCAGTACACTTTCTGCCCAGTCATCCACCAAAGA           >  1:1122133/1‑36 (MQ=255)                        GTGCAGTACACTTTCTGCCCAGTCATCCACCAAAGA           >  1:764288/1‑36 (MQ=255)                        GTGCAGTACACTTTCTGCCCAGTCATCCACCAAAGA           >  1:289823/1‑36 (MQ=255)                        GTGCAGTACACTTTCTGCCCAGTCATCCACCAAAGA           >  1:7416139/1‑36 (MQ=255)                          GCAGTACACTTTCTGCCCAGTCATCCACCAAAGAAG         <  1:6414640/36‑1 (MQ=255)                           CAGTACACTTTCTGCCCAGTCATCCACCAAAGAAGA        <  1:91716/36‑1 (MQ=255)                           CAGTACACTTTCTGCCCAGTCATCCACCAAAGAAGA        >  1:4150525/1‑36 (MQ=255)                            AGTACACTTTCTGCCCAGTCATCCACCAAAGAAGAA       <  1:3279116/36‑1 (MQ=255)                              TACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGC     >  1:3111387/1‑36 (MQ=255)                              TACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGC     >  1:1779211/1‑36 (MQ=255)                              TACACTTGCTGCCCAGTCATCCACCAAAGAAGAAGC     <  1:7300906/36‑1 (MQ=255)                               ACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCC    >  1:3060729/1‑36 (MQ=255)                               ACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCC    <  1:2628850/36‑1 (MQ=255)                               ACACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCC    <  1:1645864/36‑1 (MQ=255)                                CACTTTCTGCCCAGTCATCCACCAAAGACGAAGCCA   >  1:4535629/1‑36 (MQ=255)                                 ACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  1:6473681/1‑36 (MQ=255)                                 ACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  <  1:77937/36‑1 (MQ=255)                                 ACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  1:1928816/1‑36 (MQ=255)                                                                      CATTTCATCGTCCAACAAAATGGGTGCAGTACA...................................  <  REL606/588527‑588495 ...............................CACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAG  >  REL606/3015771‑3015805 Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG