Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,131,341	Δ128 bp	IS5‑mediated	yghQ ←	putative transport protein YghQ

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3130384–3131345	3131468	124–1085	2 [1]	[0] 29	[insH9]–[yghQ]	[insH9],[yghQ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	29 (1.100)	29/488	0.1	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	3131469 =	0 (0.000)					coding (940/1068 nt)	yghQ	putative transport protein YghQ

GGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCG  >  NC_000913/3131229‑3131712                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    gctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCAc                                                                                                                                                                                                                                           <  1:210711‑M2/11‑1 (MQ=255)      aaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgctgattcgcaccttccCTAACAACCACGGTCt                                                                                                                                                                                                                                      <  1:84989‑M2/16‑1 (MQ=255)         acgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCtggt                                                                                                                                                                                                                                   <  2:138223‑M2/19‑1 (MQ=255)           gataaccaactggcgatgttattcacgctggccaacctgtttcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCtggtgg                                                                                                                                                                                                                                 >  1:54108‑M2/231‑251 (MQ=255)                                                                   gatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaagcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGAGTAAAAACTTTTACCCAGCAAACCCGCt                                                                                                                                                                         <  2:32932‑M2/77‑1 (MQ=255)                                                                     tacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTgg                                                                                                                                                                       >  1:253767‑M2/173‑251 (MQ=255)                                                                      acgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTgg                                                                                                                                                                       <  1:67719‑M2/79‑1 (MQ=255)                                                                                                    ggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCGAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATgg                                                                                                                                          <  2:10303‑M2/108‑1 (MQ=255)                                                                                                     gggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcgtattcgccccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTGTACCCAGCAAACCCGCTGGCg                                                                                                                                                                     <  2:222161‑M2/81‑1 (MQ=255)                                                                                                     gggataaagccttaaatggcgaagaaacggtctaaataggctgattccaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCg                                                                                                                                                                     >  1:222161‑M2/141‑221 (MQ=255)                                                                                                             gccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaaggaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAAt                                                                                                                               <  1:94171‑M2/119‑1 (MQ=255)                                                                                                                   aatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattctcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCcg                                                                                                                         <  1:255081‑M2/125‑1 (MQ=255)                                                                                                                            gaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttatgcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATACCCCGGCGGCAGCg                                                                                                                <  1:11282‑M2/134‑1 (MQ=255)                                                                                                                                 ggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTcc                                                                                                           >  1:267818‑M2/113‑251 (MQ=255)                                                                                                                                    ctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaagaa                                                                                                                                                  >  2:80998‑M2/110‑209 (MQ=255)                                                                                                                                    ctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaagaa                                                                                                                                                  <  1:80998‑M2/100‑1 (MQ=255)                                                                                                                                      aaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTGCCGGCGGCGTCGAAGAATGTCAGGGCGATTGTAAATAACCCGGCGGCAGCTGGTCCTaa                                                                                                        <  2:26140‑M2/142‑1 (MQ=255)                                                                                                                                                                                          gaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGAGCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTg                                                  <  2:75460‑M2/196‑1 (MQ=255)                                                                                                                                                                                               gaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGAGCCCCACCAACACTGTGCTGCACGAGTTACGCGCGGACCAGATGGAGTGGGCa                                              <  2:54108‑M2/200‑1 (MQ=255)                                                                                                                                                                                                aaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAAt                                            >  2:148403‑M2/50‑251 (MQ=255)                                                                                                                                                                                                 aatgacggagtcagccgagaagaagttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATg                                           <  2:253767‑M2/203‑1 (MQ=255)                                                                                                                                                                                                             agccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAg                               <  2:87149‑M2/215‑1 (MQ=255)                                                                                                                                                                                                               ccgagaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGAc                             >  1:36142‑M2/35‑251 (MQ=255)                                                                                                                                                                                                                   gaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGTTTGACCAGACaaaa                         <  2:101682‑M2/221‑1 (MQ=255)                                                                                                                                                                                                                   gaagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACaaa                          >  2:106600‑M2/31‑250 (MQ=255)                                                                                                                                                                                                                    aagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGATAAAAc                        >  2:161490‑M2/30‑251 (MQ=255)                                                                                                                                                                                                                    aagaatttccccgcttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAAc                        <  2:66885‑M2/222‑1 (MQ=255)                                                                                                                                                                                                                                  ttattcgcaccttccCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAACCTCCCCGCGCCttt           >  2:247665‑M2/16‑250 (MQ=255)                                                                                                                                                                                                                                          accttccCTAACAACCACGGTCTGGTGGGGCGAGGATATAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGCAGACTGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTACTAACACGATCCCCACCAACACTGTGCTGCACGACATACGCGCCGCCCAGATGGCGTGGGCAATGTTGGTTGACCAGACAAAACACCACGCGCCTTTAATGTATCg  >  2:236990‑M2/8‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    GGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCG  >  NC_000913/3131229‑3131712

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