Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,893,961	Δ3 bp	coding (43‑45/261 nt)	ygcO →	putative 4Fe‑4S cluster‑containing protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2893961	2893963	3	47 [0]	[0] 47	ygcO	putative 4Fe‑4S cluster‑containing protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2893960	0 (0.000)	47 (1.190)	34/238	0.2	100%	coding (42/261 nt)	ygcO	putative 4Fe‑4S cluster‑containing protein
	?	NC_000913	2893964 =	0 (0.000)					coding (46/261 nt)	ygcO	putative 4Fe‑4S cluster‑containing protein

TAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAGCCCGTAATCTCTGGCGCGTTGCTGATGCGCCGCACATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCAGTTGCCTGGAGTGTGGCACCTGCC  >  NC_000913/2893826‑2894099                                                                                                                                           |                                                                                                                                        tgataaacctgttccaccgctgcgccagttattctggcatcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccga                                                                                                                                         <  1:455490‑M2/1‑1 (MQ=255)   ataaacctgttccaccgctgcgccagttattctggcatcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgAtt                                                                                                                                       >  2:774914‑M2/137‑139 (MQ=255)    taaacctgttccaccgctgcgccagttattctggcatcatttgcgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTg                                                                                                                                      >  2:523005‑M2/136‑139 (MQ=255)    taaacctgttccaccgctgcgccagttattctggcatcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTg                                                                                                                                      >  2:490187‑M2/136‑139 (MQ=255)           gttccaccgctgcgccagttattctggcatcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGc                                                                                                                               >  2:200637‑M2/129‑139 (MQ=255)           gttccaccgctgcgccagttattctggcatcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGc                                                                                                                               >  1:631573‑M2/129‑139 (MQ=255)                      gcgccagttattctggcatcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTg                                                                                                                    >  1:154056‑M2/118‑139 (MQ=255)                        gccagttattctggcatcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAg                                                                                                                  <  1:139365‑M2/24‑1 (MQ=255)                            gttattctggcatcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGcc                                                                                                              >  1:799522‑M2/112‑139 (MQ=255)                                 tctggcatcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACg                                                                                                         >  2:182627‑M2/107‑139 (MQ=255)                                      catcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGa                                                                                                    >  1:504723‑M2/102‑139 (MQ=255)                                       atcatttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGaa                                                                                                   <  1:503746‑M2/39‑1 (MQ=255)                                          atttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACgg                                                                                                >  2:136588‑M2/98‑139 (MQ=255)                                           tttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGt                                                                                               >  2:52241‑M2/97‑139 (MQ=255)                                           tttacgtcgtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGt                                                                                               >  2:165579‑M2/97‑139 (MQ=255)                                                   gtcacggcctgtggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgcggttgcgacgATTGTTCCGGCTGATTACCTTGAGCGCCAGACGGCAGCACGGTTGATTAAc                                                                                       <  1:97457‑M2/51‑1 (MQ=255)                                                              tggcatctggcgggcgatgttatcaggagtctgcgatgtctgtaacccggaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGc                                                                            <  1:526456‑M2/62‑1 (MQ=255)                                                                    ctggcgggcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCt                                                                      >  2:485636‑M2/72‑139 (MQ=255)                                                                           gcgatgttatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCg                                                               <  1:485636‑M2/75‑1 (MQ=255)                                                                                  tatcaggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCacac                                                        <  1:200637‑M2/82‑1 (MQ=255)                                                                                     caggagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCgg                                                     <  2:504723‑M2/85‑1 (MQ=255)                                                                                        gagtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAg                                                  <  1:69219‑M2/88‑1 (MQ=255)                                                                                         agtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACgg                                                                                                <  1:134313‑M2/42‑1 (MQ=255)                                                                                         agtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACgg                                                                                                >  2:134313‑M2/51‑92 (MQ=255)                                                                                         agtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCgg                                                     <  1:746488‑M2/85‑1 (MQ=255)                                                                                         agtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCgg                                                     >  2:746488‑M2/51‑135 (MQ=255)                                                                                         agtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAgg                                                 >  1:540546‑M2/51‑139 (MQ=255)                                                                                          gtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGt                                                <  1:165579‑M2/90‑1 (MQ=255)                                                                                          gtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGt                                                <  2:398098‑M2/90‑1 (MQ=255)                                                                                          gtctgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGt                                                <  1:490187‑M2/90‑1 (MQ=255)                                                                                             tgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAAc                                             >  1:252300‑M2/47‑139 (MQ=255)                                                                                             tgcgatgtctgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAAc                                             >  2:644869‑M2/47‑139 (MQ=255)                                                                                                      tgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGt                                                <  1:702277‑M2/90‑1 (MQ=255)                                                                                                      tgtaacccgtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGt                                                >  2:702277‑M2/38‑127 (MQ=255)                                                                                                              gtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCg                            <  1:590174‑M2/110‑1 (MQ=255)                                                                                                              gtaatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCg                            >  2:302992‑M2/30‑139 (MQ=255)                                                                                                                aatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCa                          <  1:440567‑M2/112‑1 (MQ=255)                                                                                                                aatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCa                          <  2:729319‑M2/112‑1 (MQ=255)                                                                                                                aatctctggcgcgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCa                          <  2:219284‑M2/112‑1 (MQ=255)                                                                                                                        gcgcgttgctgatgcgccgATTGTTCCGGCTGACTGCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCAGTTGCCTg                  <  1:720572‑M2/120‑1 (MQ=255)                                                                                                                           cgttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCAGTTGCCTGGAg               >  2:647128‑M2/17‑139 (MQ=255)                                                                                                                             ttgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCAGTTGCCTGGAgtg             >  1:312247‑M2/15‑139 (MQ=255)                                                                                                                              tgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGt                                                <  1:389087‑M2/90‑1 (MQ=255)                                                                                                                              tgctgatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGt                                                >  2:389087‑M2/14‑103 (MQ=255)                                                                                                                                  gatgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCAGTTGCCTGGAGTGTGGCa        >  2:341254‑M2/10‑139 (MQ=255)                                                                                                                                   atgcgccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTCTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCAGTTGCCTGGAGTGTGGCAc       <  1:52241‑M2/131‑1 (MQ=255)                                                                                                                                        ccgATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCAGTTGCCTGGAGTGTGGCACCTGcc  >  1:686788‑M2/4‑139 (MQ=255)                                                                                                                                           |                                                                                                                                        TAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAGCCCGTAATCTCTGGCGCGTTGCTGATGCGCCGCACATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAACTTACGTATTGACTATCGCAGTTGCCTGGAGTGTGGCACCTGCC  >  NC_000913/2893826‑2894099

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