Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,213,570	Δ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	3213570	3213617	48	51 [0]	[0] 51	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	= 3213569	0 (0.000)	47 (0.760)	41/264	0.5	100%	coding (523/1842 nt)	rpoD	RNA polymerase, sigma 70 (sigma D) factor
	?	NC_000913	3213618 =	0 (0.000)					coding (572/1842 nt)	rpoD	RNA polymerase, sigma 70 (sigma D) factor

AACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTG  >  NC_000913/3213428‑3213710                                                                                                                                              |                                                                                                                                              aaCCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCCGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAAGGCAGAAGAAGAGCTGGa                                                                                                                                              >  2:556915‑M1/1‑142 (MQ=255) aaCCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGa                                                                                                                                              >  1:896674‑M1/1‑142 (MQ=255)  aCCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGac                                                                                                                                             >  2:459399‑M1/1‑141 (MQ=255)     aGGTTCAGTGCGCCGTTGATGGATGACCGGAAGCGATCACCTAGCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgat                                                                                                                                          <  1:493687‑M1/143‑6 (MQ=255)      ggTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACTCGAACGCAGAAGGAGTTCTGGtcgatg                                                                                                                                         >  1:640548‑M1/1‑137 (MQ=255)       gTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGCAGCGCGTCGGTCCGATCTGATCACCGGCTTGGTCGACCCGAACGCAGCAGAAGAACGGGaccatga                                                                                                                                        >  2:995723‑M1/1‑136 (MQ=255)       gTTCAATGCTCCGGTGCTGAATATCCGGAAGCGAGCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatga                                                                                                                                        <  1:1087137‑M1/143‑8 (MQ=255)             tGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaaga                                                                                                                                  >  2:314446‑M1/1‑130 (MQ=255)               cTCCGTTGCTGAATATCCGGACGCGATTACCTATCTGCCGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatg                                                                                                                                <  1:320789‑M1/143‑16 (MQ=255)                  cGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCGGCACCATCTTATCTCCGTTTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaag                                                                                                                             <  1:555112‑M1/143‑19 (MQ=255)                   gctgcTGAATATCCGGGAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaaga                                                                                                                            <  2:58575‑M1/141‑20 (MQ=255)                       cTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGACGAAGATCTGGacgatgacgaagatgaagacgaa                                                                                                                        >  2:715382‑M1/1‑120 (MQ=255)                          aaTATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaa                                                                                                                     >  1:801894‑M1/1‑117 (MQ=255)                           atatCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGCCTGTCCGATCTGATCACCGGCTACGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaaaatgaaaacgaagaag                                                                                                                    >  2:911483‑M1/1‑116 (MQ=255)                            tatCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCATTGTTGGCCCGAACGCCGAAGAAGATCTCGacgatcacgaagaagaagaccacgaaga                                                                                                                   >  2:815376‑M1/1‑115 (MQ=255)                               ccGGAAGCGATCACCTATCTGCGGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGAGCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatgg                                                                                                                <  2:942787‑M1/143‑32 (MQ=255)                                 ggAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcg                                                                                                              >  2:102486‑M1/1‑110 (MQ=255)                                      cGATCACCTATCTGCTGGAACAGGACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgac                                                                                                         <  1:631142‑M1/143‑39 (MQ=255)                                             cTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcg                                                                                                  >  1:849265‑M1/1‑98 (MQ=255)                                                tCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaaggtggcgaggacgacagcgccg                                                                                               >  2:5786‑M1/1‑95 (MQ=255)                                                tCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGaacgtgtcttagatgaagacgaagaagatggcgatgacgacagcgccg                                                                                               >  1:868433‑M1/1‑95 (MQ=255)                                                 gtgctggAGCAGTACGATCGTGTTGATGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgtcgaagatgaagacgaagaagatggcgatgacgacagcgccga                                                                                              <  2:988201‑M1/142‑50 (MQ=255)                                                     tgGACCCGTAGGGTCGTGTTGCAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgat                                                                                          <  2:318385‑M1/143‑54 (MQ=255)                                                     tgGAACAGTACGATCGGGTTGAAGCAGAAGAAGGGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgat                                                                                          <  2:102107‑M1/143‑54 (MQ=255)                                                      ggAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatg                                                                                         >  1:1111607‑M1/1‑89 (MQ=255)                                                      ggAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACAGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgctg                                                                                         >  1:626054‑M1/1‑89 (MQ=255)                                                       gAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatga                                                                                        <  1:941685‑M1/143‑56 (MQ=255)                                                       gAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTGGACCCGAACGCAGAGGAAGATCTGGacgctgacgagcatggagacgaagacaatggcgaggacgacagcgcggatagtga                                                                                        >  2:92092‑M1/1‑88 (MQ=255)                                                               cGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgcagaagatggcgctggcgacagcgccgatgatgacaacagca                                                                                >  1:1151408‑M1/1‑80 (MQ=255)                                                               cGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACAGGCTTTGTTGACCCGAACGCAGAAGACGACCTGGacgaagacgacgatgacgacgaaccagatggcgacgacgaccgcgccgacgacgacaacagca                                                                                >  2:455329‑M1/1‑80 (MQ=255)                                                                gATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcat                                                                               <  1:1174327‑M1/143‑65 (MQ=255)                                                                  tCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcg                                                                             >  2:191596‑M1/1‑77 (MQ=255)                                                                   cgggttgaAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcga                                                                            <  1:370692‑M1/140‑68 (MQ=255)                                                                      gtTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgaccc                                                                         >  1:158037‑M1/1‑73 (MQ=255)                                                                                 agaagCGCGTCTGTCCGATGTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggctc                                                              <  2:874619‑M1/143‑82 (MQ=255)                                                                                       gcgTGTGTCCGATCGGATCACGGGCTTTGGGGACCCGAACGCGGAAGAAGATCTGGacgatgacgaagatgaagacgaagaggatggcgatgacgacagcgccgatgatgacaacagcatcgtcccggaactggctcgcgaaa                                                        <  2:148633‑M1/143‑88 (MQ=255)                                                                                         gTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgaccgcgccgatgatgccaacagcctcgacccggaacgggctcgcgaacaa                                                      >  2:238540‑M1/1‑54 (MQ=255)                                                                                               ccGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcg                                                <  1:459399‑M1/143‑96 (MQ=255)                                                                                                    cTGATCACGGGCTTTGTGGATCCGAACGCAGAAGAAGAGCTGGaggatgacgaagatgaagacgaagaagagggcaatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaact                                           <  1:500614‑M1/143‑101 (MQ=255)                                                                                                     tGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaacta                                          >  2:171353‑M1/1‑42 (MQ=255)                                                                                                       aTCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaaggagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaagtggctcgcgaaaaatttgcggaactacg                                        >  2:672976‑M1/1‑40 (MQ=255)                                                                                                            gggcTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaactacgcgctc                                   <  1:1033498‑M1/142‑109 (MQ=255)                                                                                                             ggCTTTGTTGACCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgagccggaactggctcgcgaacaatttgcggaactacgcgctcc                                  >  1:408869‑M1/1‑34 (MQ=255)                                                                                                                       aCCCGAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaactacgcgctcagtacgttgta                        >  2:497803‑M1/1‑24 (MQ=255)                                                                                                                           taacGCAGAAGAATATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaactacgcgctcagtacgttgtaacgc                    <  1:1070243‑M1/142‑124 (MQ=255)                                                                                                                           gAACGCAGAAGAAGATCTGGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggcccgcgaaaactttgcggaactacgcgctcagcacggtgtaacgc                    >  2:52352‑M1/1‑20 (MQ=255)                                                                                                                             aCGCAGAAGAAGACCTGGacgaggacgaagacgaagacgaagaggatggcaatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaactacgcgctcagtacgttgtaacgcgt                  <  2:510549‑M1/143‑126 (MQ=255)                                                                                                                                      agtatctTGgcgatgacgaagatgaagacgaagaagagggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaactacgcgctcagtacgttgtaacgcgtgacaccatc         <  2:338101‑M1/140‑135 (MQ=255)                                                                                                                                         acctgGacgatgacgaagatgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaactacgcgctcagtacgttgtaacgcgtgacaccatcaaa      <  1:123797‑M1/141‑138 (MQ=255)                                                                                                                                         aTCTGGacgatgacgaagacgaagacgaagaagatggcgatgacgacagcgccgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaactacgcgctcagtacgttgtaacgcgtgacaccatcaaa      <  1:53114‑M1/143‑138 (MQ=255)                                                                                                                                             ggacgatgacgaagatgaagccgaaaaagatggcgttgacgacagcggcgatgatgacaacagcatcgacccggaactggctcgcgaaaaatttgcggaactacgcgctcagtacgttgtaacgcgtgacaccatcaaagcga  <  2:1010059‑M1/143‑142 (MQ=255)                                                                                                                                              |                                                                                                                                              AACCAGGTTCAATGCTCCGTTGCTGAATATCCGGAAGCGATCACCTATCTGCTGGAACAGTACGATCGTGTTGAAGCAGAAGAAGCGCGTCTGTCCGATCTGATCACCGGCTTTGTTGACCCGAACGCAGAAGAAGATCTGGCACCTACCGCCACTCACGTCGGTTCTGAGCTTTCCCAGGAAGATCTGGACGATGACGAAGATGAAGACGAAGAAGATGGCGATGACGACAGCGCCGATGATGACAACAGCATCGACCCGGAACTGGCTCGCGAAAAATTTG  >  NC_000913/3213428‑3213710

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: ‑