Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NZ_CP009273	2,399,465	Δ7 bp	coding (650‑656/939 nt)	lrhA ←	transcriptional regulator LrhA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NZ_CP009273	2399465	2399471	7	38 [0]	[0] 38	lrhA	transcriptional regulator LrhA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 2399464	0 (0.000)	38 (0.710)	24/166	0.6	100%	coding (657/939 nt)	lrhA	transcriptional regulator LrhA
	?	NZ_CP009273	2399472 =	0 (0.000)					coding (649/939 nt)	lrhA	transcriptional regulator LrhA

CCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAAT  >  NZ_CP009273/2399385‑2399551                                                                                        |                                                                                cgcaggtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAg                                                                               >  1:1196127‑M2/88‑90 (MQ=255)      gtccgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGt                                                                          <  2:121934‑M2/8‑1 (MQ=255)         cgggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGc                                                                       <  1:82090‑M2/11‑1 (MQ=255)           ggctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCa                                                                     >  1:540294‑M2/78‑90 (MQ=255)            gctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCAt                                                                    <  2:1130603‑M2/14‑1 (MQ=255)             ctcatcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATg                                                                   <  2:1196127‑M2/15‑1 (MQ=255)               catcatctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGa                                                                 >  2:623013‑M2/74‑90 (MQ=255)                    tctcaaccggccttgccgtcacgccaagaccggctttcactgccgcccgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                            <  2:768179‑M2/22‑1 (MQ=255)                    tctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                            <  2:777677‑M2/22‑1 (MQ=255)                    tctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                            <  2:1528599‑M2/22‑1 (MQ=255)                    tctcaaccggccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                            <  1:1055016‑M2/22‑1 (MQ=255)                             gccttgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTttatt                                                   >  2:1209996‑M2/60‑90 (MQ=255)                                 tgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATg                                                                   <  2:1324016‑M2/15‑1 (MQ=255)                                 tgccgtcacgccaagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATg                                                                   >  1:1324016‑M2/56‑70 (MQ=255)                                            caagaccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAc                                    >  2:1028447‑M2/45‑90 (MQ=255)                                                 ccggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCatat                               >  1:1217667‑M2/40‑90 (MQ=255)                                                   ggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                            <  1:703378‑M2/22‑1 (MQ=255)                                                   ggctttcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATc                                                            >  2:703378‑M2/38‑59 (MQ=255)                                                        tcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                        <  2:1202990‑M2/58‑1 (MQ=255)                                                        tcactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGaaa                        <  1:894932‑M2/58‑1 (MQ=255)                                                         cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                       <  2:1215367‑M2/59‑1 (MQ=255)                                                         cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                       <  2:1217667‑M2/59‑1 (MQ=255)                                                         cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                       <  2:1263283‑M2/59‑1 (MQ=255)                                                         cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                       <  2:133020‑M2/59‑1 (MQ=255)                                                         cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                       <  1:1028447‑M2/59‑1 (MQ=255)                                                         cactgccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                       <  2:449330‑M2/59‑1 (MQ=255)                                                             gccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGc                   >  1:1093210‑M2/28‑90 (MQ=255)                                                             gccgcacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGc                   >  2:1067481‑M2/28‑90 (MQ=255)                                                                 cacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAgg               <  2:114666‑M2/67‑1 (MQ=255)                                                                 cacgaacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAgg               <  2:540294‑M2/67‑1 (MQ=255)                                                                     aacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCtt                                                       <  2:742065‑M2/27‑1 (MQ=255)                                                                     aacggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCtt                                                       >  1:742065‑M2/20‑46 (MQ=255)                                                                        ggccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATcc        <  2:205104‑M2/74‑1 (MQ=255)                                                                         gccggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCa       <  1:251019‑M2/75‑1 (MQ=255)                                                                           cggaagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGc     <  1:1441847‑M2/77‑1 (MQ=255)                                                                              aagcgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAAt  <  1:60787‑M2/80‑1 (MQ=255)                                                                                 cgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                       <  2:537850‑M2/59‑1 (MQ=255)                                                                                 cgtcgagAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAAc                       >  1:537850‑M2/8‑66 (MQ=255)                                                                                        |                                                                                CCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAAT  >  NZ_CP009273/2399385‑2399551

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