Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp		[crl]	[crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908	258683	776	24 [0]	[0] 25	[crl]	[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	24 (1.060)	19/116	0.1	100%	intergenic (+8/‑769)	crl/crl	pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers/pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers

AGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTG  >  NC_000913/257841‑257965                                                                   |                                                           aGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAa                                                           >  1:760316‑M1/1‑67 (MQ=255)  gTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAag                                                          >  1:144346‑M1/1‑66 (MQ=255)   tGGACACCCGAAGAGCAGATTGATCATAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                         <  2:584107‑M1/68‑4 (MQ=255)   tGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                         <  2:685779‑M1/68‑4 (MQ=255)   tGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                         <  2:211298‑M1/68‑4 (MQ=255)    ggACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGTAGGTAagtg                                                        <  2:292071‑M1/68‑5 (MQ=255)          ccGAAGAGCTGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaaga                                                  <  1:445901‑M1/68‑11 (MQ=255)           cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagat                                                 >  2:612416‑M1/1‑57 (MQ=255)             aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                               <  2:851088‑M1/68‑14 (MQ=255)              agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcatagataat                                              >  1:194313‑M1/1‑54 (MQ=255)                  cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgat                                          <  2:331182‑M1/68‑19 (MQ=255)                   aGATTGATCAAAAAATTTACCGCACTAGGCACGTATATTCGTGAAGGTAagtgcaaagataatcgatt                                         >  2:732657‑M1/1‑49 (MQ=255)                      ttGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctt                                      >  2:352380‑M1/1‑46 (MQ=255)                      ttGATCAAAAAATTTACCGCACTAGGACCGTATATTCGTGAAGGTAagtgcaaagataatcgattctt                                      <  1:947576‑M1/68‑23 (MQ=255)                        gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagagcaaagataatcgattctttt                                    >  1:404595‑M1/1‑44 (MQ=255)                                aaTTTACCGCACTAGTCCAGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                            <  2:47792‑M1/68‑33 (MQ=255)                                aaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                            <  1:260199‑M1/68‑33 (MQ=255)                                        gCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgt                    <  1:482819‑M1/68‑41 (MQ=255)                                            tAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                <  2:760320‑M1/68‑45 (MQ=255)                                                cccGTATATTCGTGAAGGTAagtgcaaagaaaatcgattatttttcgattgtctggctgtatgagtca            >  2:791839‑M1/1‑20 (MQ=255)                                                   gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg         <  2:865787‑M1/68‑52 (MQ=255)                                                          cGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacc  >  2:100418‑M1/1‑10 (MQ=255)                                                          cGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacc  >  1:303873‑M1/1‑10 (MQ=255)                                                          cGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacc  >  2:708288‑M1/1‑10 (MQ=255)                                                                   |                                                           AGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTG  >  NC_000913/257841‑257965

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