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

CCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTC  >  NC_000913/2867214‑2867495                                                                                                                                                 |                                                                                                                                          ccACGATTGCCATAACGGCGGGGAATTTTTACCAACAGACGCAAATTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGGGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAAACAAgggaa                                                                                                                                      >  2:503096‑M1/1‑145 (MQ=255) ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaa                                                                                                                                     <  2:181286‑M1/151‑7 (MQ=255)        tGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgag                                                                                                                               <  1:751988‑M1/150‑13 (MQ=255)         gCCATAACGGCGGGCAATTTTTAACACCAGACGCAAGTTACTCTCGATCCTCCGGCGGCGAGAGGCGACATCTCCCCGCAGTGCGCGACGCGCAAAATAAAATTCTTCTTCGGCCGTTAACAGGGGGGAATAAACaaggtaaagctgag                                                                                                                               >  2:770390/1‑137 (MQ=255)            aTAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcg                                                                                                                          <  1:86164‑M1/151‑18 (MQ=255)                      gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacac                                                                                                                >  1:648002‑M1/1‑124 (MQ=255)                       cAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacg                                                                                                               <  2:573634‑M1/151‑29 (MQ=255)                       cAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacg                                                                                                               >  1:283566‑M1/1‑123 (MQ=255)                         aTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgct                                                                                                             >  2:626796‑M1/1‑121 (MQ=255)                          tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctg                                                                                                            >  1:68235‑M1/1‑120 (MQ=255)                          tttttACCAACAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGGGAATAACCAAggtaaagctgagtcgcgtccaacacacgct                                                                                                             >  1:253833‑M1/1‑120 (MQ=255)                                    caGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccc                                                                                                  >  1:400435‑M1/1‑110 (MQ=255)                                     aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccct                                                                                                 <  2:652023‑M1/151‑43 (MQ=255)                                            aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataa                                                                                          <  1:71576‑M1/151‑50 (MQ=255)                                             gTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataac                                                                                         >  2:873350‑M1/1‑101 (MQ=255)                                               tACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataaca                                                                                        <  2:685863‑M1/150‑52 (MQ=255)                                                aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacag                                                                                       <  2:178088‑M1/150‑53 (MQ=255)                                                   ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagtt                                                                                     <  2:684333‑M1/149‑55 (MQ=255)                                                     cGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctc                                                                                 <  2:330078‑M1/151‑59 (MQ=255)                                                        tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttc                                                                              >  2:280556‑M1/1‑90 (MQ=255)                                                         catcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttc                                                                              <  1:800292‑M1/150‑62 (MQ=255)                                                           tcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggc                                                                           >  1:543008‑M1/1‑87 (MQ=255)                                                              ggcggcGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggcca                                                                         <  2:172875‑M1/150‑67 (MQ=255)                                                               gcggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCTTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaa                                                                        >  1:503021‑M1/1‑83 (MQ=255)                                                                 ggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatc                                                                     >  1:100404‑M1/1‑81 (MQ=255)                                                                    gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgtt                                                                  >  1:513092‑M1/1‑78 (MQ=255)                                                                     agagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgtt                                                                  >  2:211777‑M1/1‑77 (MQ=255)                                                                           gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactg                                                           <  2:631621‑M1/151‑81 (MQ=255)                                                                               tctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggtt                                                       >  1:374634‑M1/1‑67 (MQ=255)                                                                                  ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggtt                                                       <  1:126836‑M1/148‑85 (MQ=255)                                                                                  ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggtt                                                       <  1:443095‑M1/148‑85 (MQ=255)                                                                                      gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttc                                                >  1:605677‑M1/1‑60 (MQ=255)                                                                                      gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgtt                                                 >  1:852826‑M1/1‑60 (MQ=255)                                                                                       cAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttct                                               <  1:56226‑M1/151‑93 (MQ=255)                                                                                         gTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctac                                             <  2:584039‑M1/151‑95 (MQ=255)                                                                                            cgcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaa                                          <  2:297034‑M1/151‑98 (MQ=255)                                                                                               gACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggc                                       >  2:170905‑M1/1‑51 (MQ=255)                                                                                                aCGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggc                                       >  1:516660‑M1/1‑50 (MQ=255)                                                                                                  gcgcAAAATAAACTTCTACTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggcct                                     <  2:99242‑M1/150‑103 (MQ=255)                                                                                                     cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttc                                 <  2:797866‑M1/151‑107 (MQ=255)                                                                                                         aTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtc                              >  1:264747‑M1/1‑41 (MQ=255)                                                                                                             aCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtcaaaa                          >  1:893093‑M1/1‑37 (MQ=255)                                                                                                               ttcttcttcGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtcaaaaacc                       >  2:85684‑M1/1‑35 (MQ=255)                                                                                                                   tcttcGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtcaaaaacctca                    >  1:753320‑M1/1‑31 (MQ=255)                                                                                                                    cttcGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtcaaaaacctcaac                  >  1:914815‑M1/1‑30 (MQ=255)                                                                                                                      tcGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtcaaaaacctcaact                 <  1:512681‑M1/150‑123 (MQ=255)                                                                                                                       cGGCCGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtcaaaaacctcaactc                <  1:728915‑M1/150‑124 (MQ=255)                                                                                                                           cGTTAACAGTGGTGAATAACCAAggtaaagctgagtcgcttccaacacacgctgtgtggctccctgcgataacagttcctcgtcggccaaatcgttatcactgggttcctggtctactaatgccttttcgtcaaaaacctcaactccgtt            <  2:19343‑M1/150‑128 (MQ=255)                                                                                                                                  aGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtcaaaaacctcaactccgttctcatcaa    <  2:612740‑M1/151‑136 (MQ=255)                                                                                                                                  aGTGGTGAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtcaaaaacctcaactccgttctcatca     >  1:699959‑M1/1‑16 (MQ=255)                                                                                                                                    tggtgAATAACCAAggtaaagctgagtcgcgtccaacacacgctgtgtggctccctgcgataacagttcctcttcggccaaatcgttatcactgggttcctgttctactaaggccttttcgtcaaaaacctcaactccgttctcatcaaat  >  2:402211‑M1/1‑14 (MQ=255)                                                                                                                                                 |                                                                                                                                          CCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTC  >  NC_000913/2867214‑2867495

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