Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_017635	1,704,950	(GTTATTT)2→1	pseudogene (5322‑5328/5428 nt)	ECW_RS08400 ←	ECSE_1600 family autotransporter

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_017635	1704941	1704947	7	47 [0]	[0] 46	ECW_RS08400	ECSE_1600 family autotransporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_017635	1704948 =	0 (0.000)	40 (1.220)	32/80	0.1	100%	pseudogene (5330/5428 nt)	ECW_RS08400	ECSE_1600 family autotransporter
	?	NC_017635	= 1704940	0 (0.000)					pseudogene (5338/5428 nt)	ECW_RS08400	ECSE_1600 family autotransporter

TATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCG  >  NC_017635/1704891‑1704989                                                  |                                                  tataGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAgg                                                   >  1:344507/1‑50 (MQ=255)   taGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAgg                                                   <  1:82563/48‑1 (MQ=255)    aGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttg                                                <  1:1540611‑M2/50‑4 (MQ=255)     gCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgt                                               >  1:150988‑M2/1‑46 (MQ=255)          tATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatt                                           <  1:1970855‑M2/49‑9 (MQ=255)          tATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttat                                            >  1:1855836‑M2/1‑41 (MQ=255)           aTCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttattt                                          >  1:480708‑M2/1‑40 (MQ=255)            tCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatttac                                        >  2:1226276‑M2/1‑39 (MQ=255)            tCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttattta                                         >  2:115284‑M2/1‑39 (MQ=255)            tCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttattta                                         >  2:1100762‑M2/1‑39 (MQ=255)             cACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatttact                                       <  1:1724924‑M2/50‑13 (MQ=255)              aCCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatttactt                                      >  2:748771‑M2/1‑37 (MQ=255)              aCCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatttact                                       >  2:1295772‑M2/1‑37 (MQ=255)               ccTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatttacttt                                     >  1:737658‑M2/1‑36 (MQ=255)                cTAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatttactttc                                    <  1:1110806‑M2/50‑16 (MQ=255)                 tAGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatttactttcg                                   <  1:1768757‑M2/50‑17 (MQ=255)                  aGTTGCACACCGACATTCCCCGACAGGCTAAGGttgttatttactttcgc                                  <  1:1151567‑M2/50‑18 (MQ=255)                  aGTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatttactttcgc                                  <  1:365293‑M2/50‑18 (MQ=255)                   gTTGCACACCGACATTCCCCCACAGGCTAAGGttgttatttactttcgc                                  >  2:1357762‑M2/1‑32 (MQ=255)                      gCACACCGACATTCCCCCACAGGCTAAGGttgttatttactttcgcctc                               >  1:1903074‑M2/1‑29 (MQ=255)                       cacaCCGACATTCCCCCACAGGCTAAGGttgttatttactttcgcctcaa                             >  1:66936‑M2/1‑28 (MQ=255)                       cacaCCGACATTCCCCCACAGGCTAAGGttgttatttactttcgcctcaa                             <  2:734985‑M2/50‑23 (MQ=255)                        acacCGACATTCCCCCACAGGCTAAGGttgttatttactttcgcctcaac                            >  1:335122‑M2/1‑27 (MQ=255)                         cacCGACATTCCCCCACAGGCTAAGGttgttatttactttcgcctcaacc                           >  1:908245‑M2/1‑26 (MQ=255)                         cacCGACATTCCCCCACAGGCTAAGGttgttatttactttcgcctcaacc                           >  2:438151‑M2/1‑26 (MQ=255)                         cacCGACATTCCCCCACAGGCTAAGGttgttatttactttcgcctcaacc                           >  2:858359‑M2/1‑26 (MQ=255)                          acCGACATTCCCCCACAGGCTAAGGttgttatttactttcgcctcaaccc                          >  1:2075190‑M2/1‑25 (MQ=255)                           ccGACATTCCCCCACAGGCTAAGGttgttatttactttcgcctcaacccc                         >  1:1026719‑M2/1‑24 (MQ=255)                               cATTCCCCCACAGGCTAAGGttgttatttactttcgcctcaaccccggta                     <  2:22273‑M2/50‑31 (MQ=255)                                   cccccACAGGCTAAGGttgttatttactttcgcctcaaccccggtacgta                 <  1:704074‑M2/50‑35 (MQ=255)                                   cccccACAGGCTAAGGttgttatttactttcgcctcaaccccggtacgta                 <  2:640460‑M2/50‑35 (MQ=255)                                    ccccACAGGCTAAGGttgttatttactttcgcctcaaccccggtacgta                 >  2:67582‑M2/1‑15 (MQ=255)                                     cccACAGGCTAAGGttgttatttactttcgcctcaaccccggtacgtact               <  2:1893517‑M2/50‑37 (MQ=255)                                       cacaGGCTAAGGttgttatttactttcgcctcaaccccggtacgtacttc             >  1:1851922‑M2/1‑12 (MQ=255)                                        acaGGCTAAGGttgttatttactttcgcctcaaccccggtacgtacttca            <  2:1519762‑M2/50‑40 (MQ=255)                                         caGGCTAAGGttgttatttactttcgcctcaaccccggtacgtacttcac           <  2:1258776‑M2/50‑41 (MQ=255)                                         caGGCTAAGGttgttatttactttcgcctcaaccccggtacgtacttcac           <  1:974620‑M2/50‑41 (MQ=255)                                         caGGCTAAGGttgttatttactttcgcctcaaccccggtacgtacttcac           <  1:733384‑M2/50‑41 (MQ=255)                                         caGGCTAAGGttgttatttactttcgcctcaaccccggtacgtacttcac           >  2:21246‑M2/1‑10 (MQ=255)                                              tAAGGttgttatttactttcgcctcaaccccggtacgtacttcacc          <  2:1810670‑M2/46‑42 (MQ=255)                                               aaGGttgttatttactttcgcctcaaccccggtacgtacttcaccgagat     <  1:658050‑M2/50‑47 (MQ=255)                                                aGGttgttatttactttcgcctcaaccccggtacgtacttcaccgagatt    >  1:270952‑M2/1‑3 (MQ=255)                                                aGGttgttatttactttcgcctcaaccccggtacgtacttcaccgagat     >  2:993835‑M2/1‑3 (MQ=255)                                                aGGttgttatttactttcgcctcaaccccggtacgtacttcaccgagat     <  1:1981763‑M2/49‑47 (MQ=255)                                                aGGttgttatttactttcgcctcaaccccggtacgtacttcaccgagat     <  1:1089784‑M2/49‑47 (MQ=255)                                                 ggttgttatttactttcgcctcaaccccggtacgtacttcaccgagatta   <  1:993160‑M2/50‑49 (MQ=255)                                                  gttgttatttactttcgcctcaaccccggtacgtacttcaccgagattac  >  2:1566526‑M2/1‑1 (MQ=255)                                                  |                                                  TATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCG  >  NC_017635/1704891‑1704989

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑