Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,867,359	Δ9 bp	coding (185‑193/993 nt)	rpoS ←	RNA polymerase, sigma S (sigma 38) factor

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2867359	2867367	9	51 [0]	[0] 50	rpoS	RNA polymerase, sigma S (sigma 38) factor

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2867358	0 (0.000)	50 (0.880)	47/288	0.4	100%	coding (194/993 nt)	rpoS	RNA polymerase, sigma S (sigma 38) factor
	?	NC_000913	2867368 =	0 (0.000)					coding (184/993 nt)	rpoS	RNA polymerase, sigma S (sigma 38) factor

CCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAT  >  NC_000913/2867223‑2867504                                                                                                                                                  |                                                                                                                                         ccacgattgccataacggcggggaatttttaccaacagacgcaaattactctcgatcatccggcggcgagaggcgacatctccacgcagggcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaacaaGGGaa                                                                                                                                      >  2:503096‑M2/146‑150 (MQ=255) ccacgattgccataacggcgggcaatttttaccaccagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTaaa                                                                                                                                     <  2:181286‑M2/6‑1 (MQ=255)        tgccataacggcgggcaatttttaccaccagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAg                                                                                                                               <  1:751988‑M2/12‑1 (MQ=255)            ataacggcgggcaatttttaccaccagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTcgcg                                                                                                                          <  1:86164‑M2/17‑1 (MQ=255)                      gcaatttttaccaccagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAacacac                                                                                                                >  1:648002‑M2/125‑151 (MQ=255)                       caatttttaccaccagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACg                                                                                                               <  2:573634‑M2/28‑1 (MQ=255)                       caatttttaccaccagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACg                                                                                                               >  1:283566‑M2/124‑151 (MQ=255)                         atttttaccaccagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCt                                                                                                             >  2:626796‑M2/122‑151 (MQ=255)                          tttttaccaccagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCtg                                                                                                            >  1:68235‑M2/121‑151 (MQ=255)                          tttttaccaacagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggggaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCt                                                                                                             >  1:253833‑M2/121‑150 (MQ=255)                                    cagacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTccc                                                                                                  >  1:400435‑M2/111‑151 (MQ=255)                                     agacgcaagttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCt                                                                                                 <  2:652023‑M2/42‑1 (MQ=255)                                            agttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa                                                                                          <  1:71576‑M2/49‑1 (MQ=255)                                             gttactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                                         >  2:873350‑M2/102‑151 (MQ=255)                                               tactctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACa                                                                                        <  2:685863‑M2/51‑1 (MQ=255)                                                actctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAg                                                                                       <  2:178088‑M2/52‑1 (MQ=255)                                                   ctcgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGtt                                                                                     <  2:684333‑M2/54‑1 (MQ=255)                                                     cgatcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCctc                                                                                 <  2:330078‑M2/58‑1 (MQ=255)                                                        tcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                                              >  2:280556‑M2/91‑151 (MQ=255)                                                         catccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                                              <  1:800292‑M2/61‑1 (MQ=255)                                                           tccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGc                                                                           >  1:543008‑M2/88‑151 (MQ=255)                                                              ggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCa                                                                         <  2:172875‑M2/66‑1 (MQ=255)                                                               gcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggcctttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCaa                                                                        >  1:503021‑M2/84‑150 (MQ=255)                                                                 ggcgagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATc                                                                     >  1:100404‑M2/82‑151 (MQ=255)                                                                    gagaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGtt                                                                  >  1:513092‑M2/79‑151 (MQ=255)                                                                     agaggcgacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGtt                                                                  >  2:211777‑M2/78‑150 (MQ=255)                                                                           gacatctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTg                                                           <  2:631621‑M2/80‑1 (MQ=255)                                                                               tctccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGtt                                                       >  1:374634‑M2/68‑151 (MQ=255)                                                                                  ccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGtt                                                       <  1:126836‑M2/84‑1 (MQ=255)                                                                                  ccacgcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGtt                                                       <  1:443095‑M2/84‑1 (MQ=255)                                                                                      gcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGtt                                                 >  1:852826‑M2/61‑150 (MQ=255)                                                                                      gcagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTc                                                >  1:605677‑M2/61‑151 (MQ=255)                                                                                       cagtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTct                                               <  1:56226‑M2/92‑1 (MQ=255)                                                                                         gtgcgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctac                                             <  2:584039‑M2/94‑1 (MQ=255)                                                                                            cgcgacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTaa                                          <  2:297034‑M2/97‑1 (MQ=255)                                                                                               gacgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGc                                       >  2:170905‑M2/52‑151 (MQ=255)                                                                                                acgcgcaaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGc                                       >  1:516660‑M2/51‑150 (MQ=255)                                                                                                  gcgcaaaataaacttctacttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCt                                     <  2:99242‑M2/102‑1 (MQ=255)                                                                                                     caaaataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTc                                 <  2:797866‑M2/106‑1 (MQ=255)                                                                                                         ataaacttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTc                              >  1:264747‑M2/42‑150 (MQ=255)                                                                                                             acttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaaa                          >  1:893093‑M2/38‑150 (MQ=255)                                                                                                               ttcttcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAAcc                       >  2:85684‑M2/36‑151 (MQ=255)                                                                                                                   tcttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCa                    >  1:753320‑M2/32‑150 (MQ=255)                                                                                                                    cttcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAAc                  >  1:914815‑M2/31‑151 (MQ=255)                                                                                                                      tcggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt                 <  1:512681‑M2/122‑1 (MQ=255)                                                                                                                       cggccgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTc                <  1:728915‑M2/123‑1 (MQ=255)                                                                                                                           cgttaacagtggtgaataaccaaGGTAAAGCTGAGTCGCTTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCGTCGGCCAAATCGTTATCACTGGGTTCCTGGTCTACTAATGCCTTTTCGTCAAAAACCTCAACTCCGtt            <  2:19343‑M2/127‑1 (MQ=255)                                                                                                                                  agtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatca     >  1:699959‑M2/17‑150 (MQ=255)                                                                                                                                  agtggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa    <  2:612740‑M2/135‑1 (MQ=255)                                                                                                                                    tggtgaataaccaaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt  >  2:402211‑M2/15‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                         CCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAT  >  NC_000913/2867223‑2867504

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