Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	776,934	Δ90 bp	coding (593‑682/1266 nt)	tolA →	Tol‑Pal system protein TolA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	776934	777023	90	19 [0]	[0] 19	tolA	Tol‑Pal system protein TolA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 776933	0 (0.000)	12 (0.480)	12/254	0.9	100%	coding (592/1266 nt)	tolA	Tol‑Pal system protein TolA
	?	NC_000913	777024 =	0 (0.000)					coding (683/1266 nt)	tolA	Tol‑Pal system protein TolA

AGATGCTAAAGCTGCGGAAGAAGCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGAAGCGGCAGAAGCAGCTGCAGCTGAAGCAAGAAAGAAAGCGGCAACTGAAGCTGCTGAAAAAGCCAAAGCAGAAGCTGAGAAGAAAGCGGCTGCTGAAAAGGCTGCAGCTGATAAGAAAGCGGCAGCAGAGAAAGCTGC  >  NC_000913/776794‑777072                                                                                                                                            |                                                                                                                                            aGATGCTAAAGCTGCGGAAGAAGCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGG                                                                                                                                             <  1:246115‑M1/140‑1 (MQ=255)   aTGCTAAAGCTGCGGAAGAAGCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGct                                                                                                                                           >  1:344602‑M1/1‑138 (MQ=255)        aaaGCTGCGGAAGAAGCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctga                                                                                                                                      >  1:125832‑M1/1‑133 (MQ=255)         aaGCTGCGGAAGAAGCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaa                                                                                                                                     <  2:114248‑M1/140‑9 (MQ=255)              gCGGAAGAAGCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaaaaggc                                                                                                                                >  2:56353‑M1/1‑127 (MQ=255)               cGGAAGAAGCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaaaaggct                                                                                                                               >  1:154462‑M1/1‑126 (MQ=255)                ggAAGAAGCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaaaaggctg                                                                                                                              >  1:197991‑M1/1‑125 (MQ=255)                   agaagCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaaaaggctgcag                                                                                                                           <  2:190365‑M1/140‑19 (MQ=255)                                aaaGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaaaaggctgcagctgataagaaagc                                                                                                              >  2:117972‑M1/1‑109 (MQ=255)                                  aGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaaaaggctgcagctgataagaaagcgg                                                                                                            <  2:344602‑M1/140‑34 (MQ=255)                                                                          gccAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaaaaggctgcagctgataagaaagcggcagcagagaaagctgcagccgacaaaaaagcagcagaaaa                                                                    >  1:455776‑M1/1‑67 (MQ=255)                                                                               aGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaaaaggctgcagctgataagaaagcggcagcagagaaagctgcagccgacaaaaaagcagcagaaaaagcgg                                                               <  2:176480‑M1/140‑79 (MQ=255)                                                                                                            aGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGctgctgaaaaggctgcagctgataagaaagcggcagcagagaaagctgcagccgacaaaaaagcagcagaaaaagcggctgctgaaaaggcagcagctgataagaaa                                  >  2:217541‑M1/1‑33 (MQ=255)                                                                                                                                  gaagaaAGCGGctgctgaaaaggctgcagctgataagaaagcggcagcagagaaagctgcagccgacaaaaaagcagcagaaaaagcggctgctgaaaaggcagcagctgataagaaagcagcggcagaaaaagccgccg            <  1:56353‑M1/140‑130 (MQ=255)                                                                                                                                          cGGctgctgaaaaggctgcagctgataagaaagcggcagcagagaaagctgcagccgacaaaaaagcagcagaaaaagcggctgctgaaaaggcagcagctgataagaaagcagcggcagaaaaagccgccgcagacaaa    >  1:52150‑M1/1‑3 (MQ=255)                                                                                                                                          cGGctgctgaaaaggctgcagctgataagaaagcggcagcagagaaagctgcagccgacaaaaaagcagcagaaaaagcggctgctgaaaaggcagcagctgataagaaagcagcggcagaaaaagccgccgcagacaaa    >  2:250710‑M1/1‑3 (MQ=255)                                                                                                                                          cGGctgctgaaaaggctgcagctgataagaaagcggcagcagagaaagctgcagccgacaaaaaagcagcagaaaaagcggctgctgaaaaggcagcagctgataagaaagcagcggcagaaaaagccgccgcagacaaa    >  2:261418‑M1/1‑3 (MQ=255)                                                                                                                                          cGGctgctgaaaaggctgcagctgataagaaagcggcagcagagaaagctgcagccgacaaaaaagcagcagaaaaagcggctgctgaaaaggcagcagctgataagaaagcagcggcagaaaaagccgccgcagacaaa    >  2:301285‑M1/1‑3 (MQ=255)                                                                                                                                            gctgctgaaaaggctgcagctgataagaaagcggcagcagagaaagctgcagccgacaaaaaagcagcagaaaaagcggctgctgaaaaggcagcagctgataagaaagcagcggcagaaaaagccgccgcagacaaaaa  >  1:251670‑M1/1‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            AGATGCTAAAGCTGCGGAAGAAGCAGCGAAGAAAGCGGCTGCAGACGCAAAGAAAAAAGCAGAAGCAGAAGCCGCCAAAGCCGCAGCCGAAGCGCAGAAAAAAGCCGAGGCAGCCGCTGCGGCACTGAAGAAGAAAGCGGAAGCGGCAGAAGCAGCTGCAGCTGAAGCAAGAAAGAAAGCGGCAACTGAAGCTGCTGAAAAAGCCAAAGCAGAAGCTGAGAAGAAAGCGGCTGCTGAAAAGGCTGCAGCTGATAAGAAAGCGGCAGCAGAGAAAGCTGC  >  NC_000913/776794‑777072

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