Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,213,568	Δ48 bp	coding (524‑571/1842 nt)	rpoD →	RNA polymerase, sigma 70 (sigma D) factor

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3213568	3213615	48	51 [0]	[1] 52	rpoD	RNA polymerase, sigma 70 (sigma D) factor

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3213567	0 (0.000)	47 (0.730)	41/264	0.5	100%	coding (523/1842 nt)	rpoD	RNA polymerase, sigma 70 (sigma D) factor
	?	NC_000913	3213616 =	0 (0.000)					coding (572/1842 nt)	rpoD	RNA polymerase, sigma 70 (sigma D) factor

CTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTCAGTACGTTGTAACGCGTGACACCATCAAAGCGA  >  NC_000913/3213474‑3213756                                                                                                                                               |                                                                                                                                             aaccaggttcaatgctccgttgctgaatatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagccgaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaaggcagaagaagagctgga                                                                                                                                              >  2:556915‑M2/143‑143 (MQ=255) aaccaggttcaatgctccgttgctgaatatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctgga                                                                                                                                              >  1:896674‑M2/143‑143 (MQ=255)  accaggttcaatgctccgttgctgaatatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggAc                                                                                                                                             >  2:459399‑M2/142‑143 (MQ=255)     aggttcagtgcgccgttgatggatgaccggaagcgatcacctagctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGAt                                                                                                                                          <  1:493687‑M2/5‑1 (MQ=255)      ggttcaatgctccgttgctgaatatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgactcgaacgcagaaggagttctggTCGATg                                                                                                                                         >  1:640548‑M2/138‑143 (MQ=255)       gttcaatgctccgttgctgaatatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagcagcgcgtcggtccgatctgatcaccggcttggtcgacccgaacgcagcagaagaacgggACCATGa                                                                                                                                        >  2:995723‑M2/137‑143 (MQ=255)       gttcaatgctccggtgctgaatatccggaagcgagcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGa                                                                                                                                        <  1:1087137‑M2/7‑1 (MQ=255)             tgctccgttgctgaatatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGa                                                                                                                                  >  2:314446‑M2/131‑143 (MQ=255)               ctccgttgctgaatatccggacgcgattacctatctgccggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATg                                                                                                                                <  1:320789‑M2/15‑1 (MQ=255)                  cgttgctgaatatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtcggcaccatcttatctccgtttttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAg                                                                                                                             <  1:555112‑M2/18‑1 (MQ=255)                   gctgctgaatatccgggagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGa                                                                                                                            <  2:58575‑M2/19‑1 (MQ=255)                       ctgaatatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagacgaagatctggACGATGACGAAGATGAAGACgaa                                                                                                                        >  2:715382‑M2/121‑143 (MQ=255)                          aatatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACgaagaa                                                                                                                     >  1:801894‑M2/118‑143 (MQ=255)                           atatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgcctgtccgatctgatcaccggctacgttgacccgaacgcagaagaagatctggACGATGACGAAAATGAAAACGaagaag                                                                                                                    >  2:911483‑M2/117‑143 (MQ=255)                            tatccggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggcattgttggcccgaacgccgaagaagatctcgACGATCACGAAGAAGAAGACCAcgaaga                                                                                                                   >  2:815376‑M2/116‑143 (MQ=255)                               ccggaagcgatcacctatctgcgggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgagctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATgg                                                                                                                <  2:942787‑M2/31‑1 (MQ=255)                                 ggaagcgatcacctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCg                                                                                                              >  2:102486‑M2/111‑143 (MQ=255)                                      cgatcacctatctgctggaacaggacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATgac                                                                                                         <  1:631142‑M2/38‑1 (MQ=255)                                             ctatctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAgcg                                                                                                  >  1:849265‑M2/99‑143 (MQ=255)                                                tctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGGTGGCGAGGACGACAGCGCCg                                                                                               >  2:5786‑M2/96‑143 (MQ=255)                                                tctgctggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggAACGTGTCTTAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCg                                                                                               >  1:868433‑M2/96‑143 (MQ=255)                                                 gtgctggagcagtacgatcgtgttgatgcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGTCGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCga                                                                                              <  2:988201‑M2/49‑1 (MQ=255)                                                     tggacccgtagggtcgtgttgcagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCgatgat                                                                                          <  2:318385‑M2/53‑1 (MQ=255)                                                     tggaacagtacgatcgggttgaagcagaagaagggcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCgatgat                                                                                          <  2:102107‑M2/53‑1 (MQ=255)                                                      ggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGatgatg                                                                                         >  1:1111607‑M2/90‑143 (MQ=255)                                                      ggaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcacaggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGatgctg                                                                                         >  1:626054‑M2/90‑143 (MQ=255)                                                       gaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGAtgatga                                                                                        <  1:941685‑M2/55‑1 (MQ=255)                                                       gaacagtacgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgtggacccgaacgcagaggaagatctggACGCTGACGAGCATGGAGACGAAGACAATGGCGAGGACGACAgcgcggatagtga                                                                                        >  2:92092‑M2/89‑138 (MQ=255)                                                               cgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGCAGAAGATGGCGCTGGCGACAGCGCCGATGATGACAACAGCa                                                                                >  1:1151408‑M2/81‑143 (MQ=255)                                                               cgatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcacaggctttgttgacccgaacgcagaagacgacctggACGAAGACGACGATGACGACGAACCAGATGGCGACGACGACCGCGCCGACGACGACAACAGCa                                                                                >  2:455329‑M2/81‑143 (MQ=255)                                                                gatcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCAt                                                                               <  1:1174327‑M2/64‑1 (MQ=255)                                                                  tcgtgttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCg                                                                             >  2:191596‑M2/78‑143 (MQ=255)                                                                   cgggttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGa                                                                            <  1:370692‑M2/67‑1 (MQ=255)                                                                      gttgaagcagaagaagcgcgtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGAccc                                                                         >  1:158037‑M2/74‑143 (MQ=255)                                                                                 agaagcgcgtctgtccgatgtgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTc                                                              <  2:874619‑M2/81‑1 (MQ=255)                                                                                       gcgtgtgtccgatcggatcacgggctttggggacccgaacgcggaagaagatctggACGATGACGAAGATGAAGACGAAGAGGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGTCCCGGAACTGGCTCGCGaaa                                                        <  2:148633‑M2/87‑1 (MQ=255)                                                                                         gtctgtccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACCGCGCCGATGATGCCAACAGCCTCGACCCGGAACGGGCTCGCGaacaa                                                      >  2:238540‑M2/55‑143 (MQ=255)                                                                                               ccgatctgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCg                                                <  1:459399‑M2/95‑1 (MQ=255)                                                                                                    ctgatcacgggctttgtggatccgaacgcagaagaagagctggAGGATGACGAAGATGAAGACGAAGAAGAGGGCAATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACt                                           <  1:500614‑M2/100‑1 (MQ=255)                                                                                                     tgatcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTa                                          >  2:171353‑M2/43‑143 (MQ=255)                                                                                                       atcaccggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGGAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAAGTGGCTCGCGAAAAATTTGCGGAACTAcg                                        >  2:672976‑M2/41‑143 (MQ=255)                                                                                                            gggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTc                                   <  1:1033498‑M2/108‑1 (MQ=255)                                                                                                             ggctttgttgacccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGAGCCGGAACTGGCTCGCGAACAATTTGCGGAACTACGCGCtcc                                  >  1:408869‑M2/35‑142 (MQ=255)                                                                                                                       acccgaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTCAGTACGTTGTa                        >  2:497803‑M2/25‑143 (MQ=255)                                                                                                                           taacgcagaagaatatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTCAGTACGTTGTAAcgc                    <  1:1070243‑M2/123‑1 (MQ=255)                                                                                                                           gaacgcagaagaagatctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCCCGCGAAAACTTTGCGGAACTACGCGCTCAGCACGGTGTAAcgc                    >  2:52352‑M2/21‑143 (MQ=255)                                                                                                                             acgcagaagaagacctggACGAGGACGAAGACGAAGACGAAGAGGATGGCAATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTCAGTACGTTGTAACGCGt                  <  2:510549‑M2/125‑1 (MQ=255)                                                                                                                                      agtatcttgGCGATGACGAAGATGAAGACGAAGAAGAGGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTCAGTACGTTGTAACGCGTGACACCATc         <  2:338101‑M2/134‑1 (MQ=255)                                                                                                                                         atctggACGATGACGAAGACGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTCAGTACGTTGTAACGCGTGACACCATCaaa      <  1:53114‑M2/137‑1 (MQ=255)                                                                                                                                         acctggACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTCAGTACGTTGTAACGCGTGACACCATCaaa      <  1:123797‑M2/137‑1 (MQ=255)                                                                                                                                            tGGACGATGACGAAGATGAAGACGCAGACGCTGGCGATGACGACAGCGCAGATGATGACAACAGCGTCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAGCTAGCCCGTCAGTACCTTCGAACGGGTCGCACCagcccaccc   >  2:361312/1‑135 (MQ=255)                                                                                                                                             ggACGATGACGAAGATGAAGCCGAAAAAGATGGCGTTGACGACAGCGGCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTCAGTACGTTGTAACGCGTGACACCATCAAAGCGa  <  2:1010059‑M2/141‑1 (MQ=255)                                                                                                                                               |                                                                                                                                             CTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTGCGGAACTACGCGCTCAGTACGTTGTAACGCGTGACACCATCAAAGCGA  >  NC_000913/3213474‑3213756

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