Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_007779	3,539,230	3461 bp→CATATGTCGACTCTAGAGGATCCGCATGC		rhaB–[rhaD]	rhaB, rhaA, [rhaD]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_007779	3539230	3542690	3461	39 [0]	[0] 42	rhaB–[rhaD]	rhaB,rhaA,[rhaD]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_007779	= 3539229	0 (0.000)	29 (0.880) +CATATGTCGACTCTAGAGGATCCGCATGC	27/142	0.3	100%	intergenic (‑284/‑4)	rhaS/rhaB	DNA‑binding transcriptional activator/rhamnulokinase
	?	NC_007779	3542691 =	0 (0.000)					coding (283/825 nt)	rhaD	rhamnulose‑1‑phosphate aldolase

CTTTACTGGCAAATACACCGTTTATTGTCACCGGCTCGGGCAAATTCTTCCGTAACGTCCAGCTTGATCCTGCGGCTAACTTAGGCATCGTAAAAGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTCCTTTCCCACTGCGA  >  NC_007779/3542596‑3542788                                                                                                |                                                                                                  agtaacgagaaggtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGAc                                                                                              >  1:243775‑M2/96‑101 (MQ=255)   taacgagaaggtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAg                                                                                            <  1:664120‑M2/8‑1 (MQ=255)   taacgagaaggtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAg                                                                                            >  1:340311‑M2/94‑101 (MQ=255)     acgagaaggtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCg                                                                                          <  3:846510‑M2/10‑1 (MQ=255)     acgagaaggtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCg                                                                                          >  3:1024488‑M2/92‑101 (MQ=255)      cgagaaggtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGa                                                                                         <  1:1290044‑M2/11‑1 (MQ=255)          aaggtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGgc                                                                                     <  3:311085‑M2/15‑1 (MQ=255)          aaggtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGgc                                                                                     <  3:862788‑M2/15‑1 (MQ=255)             gtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGcgcg                                                                                  >  1:333412‑M2/84‑101 (MQ=255)             gtcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGcgcg                                                                                  <  4:81404‑M2/18‑1 (MQ=255)              tcgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCgg                                                                                 <  3:1328643‑M2/19‑1 (MQ=255)               cgcgaattcaggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCggg                                                                                >  3:87859‑M2/82‑101 (MQ=255)                 cgaattccggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCt                                                                              >  4:148774‑M2/80‑101 (MQ=255)                        aggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACAt                                                                       <  3:646724‑M2/29‑1 (MQ=255)                         ggcgctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgcctctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACAtt                                                                      <  1:511797‑M2/30‑1 (MQ=255)                            gctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCtt                                                                   <  1:50870‑M2/33‑1 (MQ=255)                             ctttttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCttt                                                                  >  4:87528‑M2/68‑101 (MQ=255)                                tttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTggg                                                               >  3:1216190‑M2/65‑101 (MQ=255)                                tttagactggtcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTggg                                                               <  1:762662‑M2/37‑1 (MQ=255)                                        gggcgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAAcc                                                       <  1:969894‑M2/45‑1 (MQ=255)                                           cgtaatgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAAc                                                    <  1:872493‑M2/48‑1 (MQ=255)                                               atgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAg                                                <  3:1215943‑M2/52‑1 (MQ=255)                                                tgaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGc                                               <  1:979247‑M2/53‑1 (MQ=255)                                                 gaaattcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGcc                                              <  1:175080‑M2/54‑1 (MQ=255)                                                    attcagcaggatcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTc                                           <  1:1339680‑M2/57‑1 (MQ=255)                                                              atctccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCg                                 <  1:71305‑M2/67‑1 (MQ=255)                                                               tcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGa                                >  3:1189526‑M2/34‑101 (MQ=255)                                                               tcaccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGa                                >  1:682477‑M2/34‑101 (MQ=255)                                                                caccatatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGaa                               >  1:751247‑M2/33‑101 (MQ=255)                                                                     tatgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTcc                          >  1:434136‑M2/28‑101 (MQ=255)                                                                      atgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCg                         <  1:253412‑M2/75‑1 (MQ=255)                                                                      atgtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCg                         <  2:148902‑M2/75‑1 (MQ=255)                                                                        gtcgactctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGc                       >  1:1195973‑M2/25‑101 (MQ=255)                                                                               ctagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTc                >  1:406627‑M2/18‑101 (MQ=255)                                                                                tagaggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTcc               >  3:1415589‑M2/17‑101 (MQ=255)                                                                                   aggatccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTCCttt            <  1:1014741‑M2/88‑1 (MQ=255)                                                                                        ccgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTCCTTTCCCAc       >  3:1263486‑M2/9‑101 (MQ=255)                                                                                         cgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTCCTTTCCCACt      >  3:1251201‑M2/8‑101 (MQ=255)                                                                                         cgcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTCCTTTCCCACt      <  3:1136060‑M2/94‑1 (MQ=255)                                                                                          gcatgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTCCTTTCCCACTg     >  3:90438‑M2/7‑101 (MQ=255)                                                                                           catgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTCCTTTCCCACTGc    <  1:110819‑M2/96‑1 (MQ=255)                                                                                             tgcGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTCCTTTCCCACTGCGa  >  1:214304‑M2/4‑101 (MQ=255)                                                                                                |                                                                                                  CTTTACTGGCAAATACACCGTTTATTGTCACCGGCTCGGGCAAATTCTTCCGTAACGTCCAGCTTGATCCTGCGGCTAACTTAGGCATCGTAAAAGTCGACAGCGACGGCGCGGGCTACCACATTCTTTGGGGGTTAACCAACGAAGCCGTCCCCACTTCCGAACTTCCGGCTCACTTCCTTTCCCACTGCGA  >  NC_007779/3542596‑3542788

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