Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	CP009273	2,093,044	Δ2,186 bp	100%	IS5‑mediated	[ugd]–wbbL	[ugd], gnd, wbbL

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	CP009273	2093044	2096362–2095227	2184–3319	78 [1]	[26] 27	[ugd]–[insH1]	[ugd],gnd,wbbL,[insH1]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 2093043	1 (0.020)	54 (0.900)	51/276	NT	98.2%	coding (52/1167 nt)	ugd	UDP‑glucose 6‑dehydrogenase
	?	CP009273	2095230 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region

CAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTC  >  CP009273/2096282‑2096512                                                                                  |                                                                                                                                                      cAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGTTCCAGATGACAAATATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATcc                                                                                   >  1:156484/1‑150 (MQ=14)     aaaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGACCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTaa                                                                               <  1:135784/150‑1 (MQ=21)         tctcTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACAt                                                                            >  2:327515/1‑149 (MQ=37)           tctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATaaa                                                                         >  2:854688/1‑150 (MQ=39)                tCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGac                                                                     <  1:553728/149‑1 (MQ=255)                tCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGa                                                                      >  1:248854/1‑148 (MQ=255)                 cTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGac                                                                     <  2:248854/148‑1 (MQ=255)                       tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATaa                                                              >  1:891132/1‑149 (MQ=255)                       tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATaa                                                              >  1:112916/1‑149 (MQ=255)                         aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATaaa                                                             <  2:867109/148‑1 (MQ=255)                           ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCa                                                         <  1:750599/150‑1 (MQ=255)                            gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAAt                                                           >  1:742103/1‑147 (MQ=255)                            gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAg                                                        <  1:29290/150‑1 (MQ=255)                              gcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGaa                                                      <  2:714803/150‑1 (MQ=255)                               cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGa                                                       >  1:695871/1‑148 (MQ=255)                                  aCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGc                                                   <  1:302733/149‑1 (MQ=255)                                       tgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAgg                                              >  2:630768/1‑149 (MQ=255)                                                  tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTAc                                  <  1:866738/150‑1 (MQ=255)                                                  tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTAc                                  <  2:1035109/150‑1 (MQ=255)                                                       gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGaa                              >  2:106806/1‑149 (MQ=255)                                                              aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCAt                       <  2:990331/149‑1 (MQ=255)                                                                   atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTa                    >  2:396568/1‑147 (MQ=255)                                                                    tgaCTCATGATGACCCCTGTTCTATGGCCCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTa                    >  2:264043/1‑146 (MQ=255)                                                                       cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGAttt             >  2:456837/1‑150 (MQ=255)                                                                        tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTc            >  1:245264/1‑150 (MQ=255)                                                                                 aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCAtct   <  1:158275/150‑1 (MQ=255)                                                                                  cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTAGCACCTTCCTTAATCCCTAACATAAATAACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCAtctc  <  2:156484/150‑1 (MQ=255)                                                                                  |                                                                                                                                                      CAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTC  >  CP009273/2096282‑2096512

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